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Chapter 1 

Integrable dynamical systems 



The usually accepted definition of an integrable system in the sense of Liouville is a 
system with phase space of dimension 2n for which one knows n conserved quantities in 
involution. This a rather puzzling definition since by Darboux theorem one can always 
find locally a system of canonical coordinates on phase space (Pi, • • • , P n ; Qi, ■ • • , Q n ) 
with H = P\ hence fulfilling the hypothesis of Liouville theorem. However the generic 
dynamical system is certainely not what we mean by integrable system, so the hypothesis 
must be made more precise by requiring some global existence properties of the conserved 
quantities. A good starting point is to ask, following Moser, that the conserved quantities 
exist on an open domain of the phase space invariant under the dynamical flow, that is 
any trajectory starting in the domain stays in it. 

In all the examples that we will consider the conserved quantities are even analytic 
functions of canonical coordinates on some open domain and the known solutions are 
similarly analytic. 



1.1 The Liouville theorem 

We consider a dynamical hamiltonian system with phase space M, dim M = 2n. In- 
troduce canonical coordinates pi , qi such that the non degenerate Poisson bracket reads 
{.Pit Qj} = As usual a non degenerate Poisson bracket on M is equivalent to the data 
of a non-degenerate closed 2-form to, duj = 0, defined on M. In the canonical coordinates 
uj = J2j dpj A dqj. Let H be the hamiltonian of the system. For any function / on M, 
the equations of motion are Hamilton's equations: 

Here and in the following, a dot will refer to a time derivative. 

Definition 1 The system is Liouville integrable if it possesses n independent conserved 
quantities Fi, i = 1, ■ ■ ■ , n, {H, Fj} = 0, in involution, 

{Fi,Fj} = 
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There cannot be more than n independent quantities in involution otherwise the Poisson 
bracket would be degenerate. In particular, the hamiltonian H is a function of the Fj's. 

Theorem 1 (The Liouville theorem.) The solution of the equations of motion of an 
integrable system is obtained by quadrature. 

Proof . Let a = Yli Pidqi be the canonical 1-form and uj = da = J2i dpi A dqi be the 
symplectic 2-form on the phase space M. We will construct a canonical transformation 
(Pii Qi) ~^ (Fi, such that the Fj's are among the new coordinates, i.e., a transforma- 
tion such that 

lo = dpi A dqi = dFi A d$>i 

i i 

If we succeed to do that, the equations of motion become trivial: 

F j = {H,F j } = 

a it 

i>j = {HM = — = %(F). 

Their solutions are: 



Fj(t) 


Fj(0) 




= Vj(o) + t%. 



To construct this canonical transformation, we exhibit its generating function S. Let 
Mf be the level manifold Fi(p,q) = Suppose we can solve for pi, p, t = Pi(f,q), and 
consider the function 



rm pq 

S(F,q)= a= ^2pi(f,q)dqi, 
J ran J an 



' mo J go 

where the integration path is drawn on Mf and goes from the point of coordinate 
{p{f,qo),qo) to the point (p(f,q),q), where qo is some reference value. 

If this function exists, i.e. if it does not depend on the path from mo to m, it is the 
function we are looking for. Indeed, from the definition of S, pj = J^. Defining tpj by 

as 

we have 

dS = iPjdFj + Pjdqj, 

3 

Since d 2 S = we deduce that u> = J2j dpj A dqj = J2j dFj A dijjj. This shows that if S 
is a well defined function then the transformation is canonical. 
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Figure 1.1: Integration path on the level manifold Mf. 



To show that S exists, we must prove that it is independent of the integration path, ie. 
we have to prove that 

da\ Mf = w\mj = 

Let Xi be the hamiltonian vector field associated to Fj, defined by dFi = u(Xi, •), 

^ dq k dp k dp k dq k ' 
These vector fields are tangent to the manifold Mf because the Fj are in involution, 

X i (F j ) = {F i ,F j } = 

Since the Fj are assumed to be independent functions, the tangent space to the subman- 
ifold Mf is generated at each point m £ M by the vectors X{\ m (i = 1, . . . ,n). But then 
uj(Xi,Xj) = dFi(Xj) = and we have proved that u\mj = 0, and therefore S exists. ■ 

Remark. From the closedness of a on Mf, the function S is unchanged by contin- 
uous deformations of the path (mo,m). However, if Mf has non trivial cycles, which is 
generically the case, S is a multivalued function defined in a neighborhood of Mf. The 
variation over a cycle 



AcydeS = a 
J cycle 

is a function of F only. This induces a multivaluedness of the variables ipj-. A cyc i e if)j = 

d A c 
dFj cycled- 



1.2 Lax pairs 



The new concept which emerged from the modern studies of integrable systems is the 
notion of Lax pairs. A Lax pair L, M consists of two functions on the phase space of the 
system, with values in some Lie algebra Q, such that the hamiltonian evolution equations 
may be written as 



dL 

~dt 



[M,L] 



(1.1) 

Here, [ , ] denotes the bracket in the Lie algebra Q. We will denote by G the connected 
Lie group having Q as a Lie algebra. 

Although eq.(l.l) only requires a Lie algebra structure to be written down, we usually 
use finite dimensional representations of Q, so that L and M are matrices. However it 
will be important to keep in mind the more abstract formulation. 

The interest in the existence of such a pair originates in the fact that it allows for an 
easy construction of conserved quantities. Indeed the solution of eq(l.l) is of the form 

L(t)=g(t)L(0)g-\t) 

where g(t) G G is determined by the equation 

dg 
dt 



M 
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It follows that the eigenvalues of L are conserved. We say that the evolution equation 
(1.1) is isopsectral which means that the spectrum of L is preserved by the time evolution. 
Alternatively the quantities 



Hn 



Tr (L n ) 



are conserved. 



Integrability of the system in the sense of Liouville demands (i) that the system is 
Hamiltonian, (ii) that the number of independent conserved quantities equals the number 
of degree of freedom, and (iii) that these conserved quantities are in involution. 

The interest in the concept of Lax pairs relies on the existence of a tool allowing to 
produce such pairs fulfilling these constraints. 



1.3 The Zakharov-Shabat construction. 

Given an integrable system, there does not yet exists a useful algorithm to construct a 
Lax pair. There does exist however a general procedure, due to Zakharov and Shabat, to 
construct consistent Lax pairs giving rise to integrable systems. This is a general method 
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to construct matrices L(X) and M(A), depending on a spectral parameter A G C, such 
that the Lax equation 

d t L(X) = [M(X),L(X)} (1.2) 

is equivalent to the equations of motion of an integrable system. The main result is 
eq.(1.13) expressing the possible forms of the matrix M in the Lax pair. 

Let us consider matrices L(X) and M(A) of dimension N x N. We will assume that 
the matrices L(A) and M(X) are rational functions of the parameter A. Let {e^} be 
the set of their poles, namely the poles of L(A) and those of M(A). With the above 
notations, assuming no pole at infinity, we can write quite generally: 

-l 

L(X) = L + Y,L k (X), with L k (X)^ L kA*-tkY (1-3) 

k r=-n k 

and 

-i 

M(A) = M + Y,M k (X) with M k (X) = M ktr (X-e k ) r (1.4) 

fc r=—m k 

Here and m*; refer to the order of the poles at the corresponding point e k . The 
coefficients L kjT and Mfc jT , are matrices. We will assume that the positions of the poles 
e k are constants independent of time. 

We now want to impose that the Lax equation (1.2), with L(A) and M(X) given by 
eqs.(1.3,1.4), holds identically in A. It is important to realize that this is a very non 
trivial equation. Indeed looking at eqs.(1.2) we see that the pole at e k in the left hand 
side is a priori of order n k while in the right hand side it is potentially of order n k + m k . 
Hence we have two types of equations. The first type does not contain time derivatives 
and comes from setting to zero the coefficients of the poles of order greater than n k in 
the right hand side of the equation. This will be interpreted as m k constraint equations 
on M k . The equations of the second type are obtained by matching the coefficients of 
the poles of order less or equal to n k on both sides of the equation. These equations 
contain time derivatives and are thus the true dynamical equations. It turns out that 
one can solve the constraints equations. 

We introduce a notation. For any matrix valued rational function /(A) with poles of 
order n k at points e k at finite distance, we can decompose /(A) as 

-l 

/(A) = /o + ^/fc(A), with f k (X)= Y /fe,r(A-e fe ) r , 

k r=-n k 

with /o a constant. The quantity /fc(A) is called the polar part at e k . When there is no 
ambiguity about the pole we are considering, we will often use the alternative notation 
/_(A) = fk(X). Around one of the point e&, /(A) may be decomposed as follows: 

/(A) = /(A) + + /(A)_ (1.5) 
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with /(A)+ regular at the point and /(A)_ = /fc(A) is the polar part. 

Assuming that L(X) has distinct eigenvalues in a neighbourhood of e^, one can per- 
form a regular similarity transformation ff^(A) diagonalizing L(A) in a vicinity of e&. 



L(A)= 5 ( fc )(A) A«(A) ffW-^A) 



(1.6) 



where A^ k \\) is diagonal and has a pole of order nj, at e^. Obviously, we can write the 
polar decomposition of L(A) as 




(1.7) 

A first consequence of the Lax equation is that M(A) admits a similar polar decompo- 
sition 

Proposition 1 The decomposition of M(A) m po/ar parts reads 



(1.8) 



where B^ k \\) is diagonal and has a pole of order mk at e&. 
Proof . Defining B^ k \\) by 

M(A) = (A) B«(A) 5 ( fe )- 1 (A) + %«(A) ^(A) 



(1.9) 



the Lax equation becomes: 



i( fe )(A) = [B«(A),i (fc) (A)] 



This implies = as expected (because the commutator with a diagonal matrix has 

no element on the diagonal) , and moreover if we assume that the diagonal elements of 
A (k) 

are all distinct this equation implies that is also diagonal. Finally the term 
dtg^ g^^ 1 is regular and does not contribute to the singular part of M at e&. Hence 
Mk = (gW B^ g^^ 1 )- which only depends on . This simultaneous diagonalization 
of L(A) and M(A) works around any point where L(A) has distinct eigenvalues. 

This proposition clarifies the structure of the Lax pair. Only the singular parts of 

and B^ contribute to and M^. The independent parameters in L(X) are thus 
Lq , the singular diagonal matrices A^ of the form 



A - )= E ^,r(A-e fe y 



(1.10) 
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and jets of regular matrices of order n k — 1, denned up to right multiplication by a 
regular diagonal matrix dP*\\): 



n fc -l 
r=0 



(1.11) 



From these data, we can reconstruct the Lax matrix L(A) by defining L = Lq + 
with 

L k =(gV>AV>gW-<)_ (1.12) 

Then around each e&, one can diagonalize £(A) = gW.AMgW -1 . This yields an extension 
of the matrices and to complete series and </ fe ) in (A — e&). Finally to 
define M(A) = Mo + £^ fc -^fe> we cnoose a se t of diagonal polar matrices (E>( k \X))- and 
use the series g^) to define by eq.(1.8). 

In the vicinity of a singularity, L(X) and M(A) can be simultaneously diagonalized 
if the Lax equation holds true. In this diagonal gauge, the Lax equation simply states 
that the matrix ^4^ fe ^(A) is conserved and that B^ k \X) is diagonal. When we transform 
these results into the original gauge, we get the general solution of the non dynamical 
constraints on M(A): 

Proposition 2 Let L(\) be a Lax matrix of the form eq.(l.S). The general form of the 
matrix M(A) such that the orders of the poles match on both sides of the Lax equation 
is M = M + Y, k M k with 



M, 



(p {k) (L,X)) 



(1.13) 

where P^ k \L,X) is a polynomial in L(X) with coefficients rational in X and ( )_ denotes 
the singular part at X = e k . 

Proof . It is easy to show that this is indeed a solution. We have to check that the order 
of the poles is correct. Let us look at what happens around A = e k - Using a beautiful 
argument first introduced by Gelfand and Dickey we write: 



[M k ,L\. 



(pW(L,\j)_,L 
?W(L,A)-(P W (L,A)) ,L 



(pW(L,A)) ,L 



where we used that a polynomial in L commutes with L. From this we see that the order 
of the pole at e k is less than n k . To show that this is a general solution, recall eqs. (1.6, 
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1.8). Since A^ k \X) is a diagonal N x N matrix with all its elements distinct in a vicinity 
of ek , its powers up to N — 1 span the space of diagonal matrices and one can write 

B W = pW(#) )A ) (1.14) 

where P^(A^ k \X) is a polynomial of degree N — 1 in A^ k \ The coefficients of P^ 
are rational combinations of the matrix elements of A^ and hence admit Laurent 
expansions in A — in a vicinity of e^. Inserting eq. (1.14) into eq. (1.8) one gets 

Mk = (^p( k \L, X)j . Moreover in this formula the Laurent expansions of the coefficients 

of P< fc ) can be truncated at some positive power of A — since a high enough power 
cannot contribute to the singular part, yielding a polynomial with coefficients Laurent 
polynomials in A — e&. 

The above propositions give the general form of M(A) as far as the matrix structure 
and the A-dependence is concerned. One should keep in mind however that the coeffi- 
cients of the polynomials p( k \L, A) are a priori functions of the matrix elements of L 
and require further characterizations in order to get an integrable system. In the setting 
of the next section these coefficients will be constants. 

Remark. The Lax equation is invariant under similarity transformations, 

L - U = gLg-\ M - M 1 = gMg' 1 + d t gg- 1 (1.15) 

If this similarity transformation is independent of A, it will not spoil the analytic prop- 
erties of L(X) and M(A). We can use the gauge freedom eq.(1.15) to diagonalize Lq, 

Lq = Diag(ai, • • -,a N ) 

Consistency of eq.(1.2) then requires Mq to be also diagonal and thus Lq = [Mq, Lq] = 0. 
Hence Mq is a polynomial P of Lq, so that replacing M(A) — > M(A) — P(L(X)) gets rid 

Of Mq. 

1.4 Coadjoint Orbits. 

(k) 

In this section we show that the Zakharov-Shabat construction, when the matrices A_ 
are non dynamical, can be interpreted as coadjoint orbits. This introduces a natural 
symplectic structure in the problem and gives a Hamiltonian interpretation to the Lax 
equation. 

Let G be a connected Lie group with Lie algebra Q. The group G acts on Q by the 
adjoint action denoted Ad: 

X — (Ad g)(X)=gXg- 1 g G G, X G Q 

Similarly the coadjoint action of G on the dual Q* of the Lie algebra Q (i.e. the vector 
space of linear forms on the Lie algebra) is defined by: 

(Ad* g.E)(X) = E(Adg-\X)), g G G, S G Q\ X G Q 



14 



The infinitesimal version of these actions provides actions of the Lie algebra Q on Q and 
Q* denoted ad and ad* respectively and given by: 



adX{Y) = [X,Y] X,Y£G, 
ad*XS(Y) = -E([X,Y]) X,Y e G,E e G* 

On the space T{Q*) of functions on Q* there is a canonical Poisson bracket called 
the Kostant-Kirillov bracket. Let S G Q* and X, Y G Q, we define 



{E(X),E(Y)} = E([X,Y}) 



(1.16) 

If e a is a basis of Q and e* a is the dual basis of Q* , then we have 

x = Y J X a e a eg, E = Y J ^ a e* a eg* 

a a 

and 

E(X)=^ a Xa 

Setting X = e a , Y = ei, in eq.(1.16) we find 

{^aj^b} = fabc^c 

where we have introduced the structure constants of the Lie algebra Q 

[fid! Sfe] — fabc&c 

This formula define the Poisson bracket of the coordinates E a on Q* . We can extend it 
on the functions on Q* (i.e. functions of the E a ) in the usual way 

i F r\m- V dF dG ^ ~\- t dF dG 

{F,G}^) -^^a^i-a,^} -2.-/^^ 



Introducing the differentials dF, dG £ Q as 

dE 
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the above formula can be rewritten in the more invariant way 

{F,G}(E) = E([dF,dG}) 

The Kostant-Kirillov bracket is degenerate, meaning that there exists functions which 
Poisson commutes with everything. For instance if the basis e a is chosen so that the 
structure constants f a bc are totally antisymmetric the function 



- 2 = \ " - 2 
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is in the center of the Poisson bracket. Indeed 



{S 2 , E b } = 2 fabc^a^c = 

a,c 

In the context of Hamiltonian mechanics, it is important to identify all such functions 
and to set them to constants since they cannot contribute to the dynamics. This is 
where the notion of coadjoint orbit plays a very important role. The coadjoint orbit of 
an element A G Q* is the set of elements of Q* defined as 

Orbit(A) = {Ad* • A,Vg G G} 

The center of the Kost ant-Kir illov bracket consists of the functions which are Ad*- 
invariant, i.e. which are constant on coadjoint orbits. In fact for such a function we 
have 

F(E) = F(Ad* g ■ E) 
Taking an infinitesimal g = 1 + eX, this translates to 

F{E) = F(E + e ad*X ■ E) = F(E) + e ad*X ■ E(dF) + 0(e 2 ) 

hence, the Ad*-invariance of F(E) can be written as 

ad*X • E(dF) = E([dF, X]) = 0, \/X G Q 
On the other hand, if F(E) is in the kernel of the Kostant-Kirillov bracket we have 

{F(E),E(X)} = E([dF,X]) = 0, MX 

so that the Ad*-invariant functions are in the kernel and vice versa. On coadjoint orbits, 
the Kostant-Kirillov bracket becomes non degenerate. 

To see how these notions relate to our problem, let us consider first a Lax matrix 
with only one polar singularity at A = 0: 

L(A) = (g{\) A_(A) g-\\)) _ (1.17) 

with A_(A) = Y^r=-n ^rA 7 ', and g{\) has a regular expansion around A = 0. 

Let G be the loop group of invertible matrix valued power series expansion around 
A = 0. The elements of G are regular series g{\) = Y1T=q9^^ T - ^ ne product law is 
the pointwise product: (gh)(X) = g(X)h(X). Formally, the Lie algebra Q of G consists 
of elements of the form X = J2^Lo-^rX r . Its Lie bracket is given by the pointwise 
commutator. 
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The dual Q* of Q can be identified with the set of polar matrices £(A) = J2 r >i E r A 
where the sum contains a finite but arbitrary large number of terms, by the pairing: 



(E,X) = Tr Res A=0 (E(A)X(A)) = J^Tr (E r+l X r ) 



where ResA=o is defined to be the coefficient of A . 

The coadjoint action of G on Q* is defined by {{Ad*g) ■ E) (X) = E(g~ 1 Xg) for 
H 6 CP and any X G Using the above model for Q*, and since (E,g~ l Xg) = 
(gEg-\X) = ({gEg- l )_,X), we get 



(Ad* ff )-H(A) = ( 5 -S- 5 - 1 )_ 



This is precisely eq.(1.17). The Lax matrix can thus be interpreted as belonging to the 
coadjoint orbit of the element A_ (A) of Q* under the loop group G. 

L(A) = (g{\) A_{\) g-\\))_ = Orbit(A_(A)) 

If we take any element of the orbit and try to diagonalize it, the singular part of the 
matrix of eigenvalues is precisely A- (A) which is therefore an Ad*-invariant function and 
should be put to constants. This interpretation of L{\) as a coadjoint orbit therefore 
assumes that A- (A) is not a dynamical variable. The coadjoint orbit is then equipped 
with the Kostant-Kirillov symplectic structure. 

This construction can be extended to the multi-pole case. We consider the direct 
sum of loop algebras Qk-, around A = e&: 

Q = ®Q k 

k 

An element of this Lie algebra has the form of a multiplet 

X(A) = (A 1 (A),X 2 (A),---) 

where Xk(X), defined around e^, is of the form Xk(X) = J2 n >o^k,n (A — e^)™. The Lie 
bracket is such that [Xfc(A), X/(A)] = if k ^ I. The group G is the direct product of 
the groups G^ of regular invertible matrices at e^: 

G = (G 1 ,G 2 ,---) (1.18) 

The dual Q* of this Lie algebra consists of multiplets 

S=(S 1 (A), 32(A), ■■■) 
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where E k (X) around e k is of the form E k (X) = J2 r >i ^k,r( x ~ e fc) T • In this sum the 
number of terms is finite but arbitrary. The pairing is simply 

(E,X) = Y,(*k,X k ) = ]TTr Re Sek (E k (X)X k (X)) 
k k 

The coadjoint action of G on Q* is given by the usual formula: if g = (51, 52, ■ • •) € G 
and S = (Hi,H 2 ,-") G 

(Afs).E(A) = ((^H^r 1 )., (^jfe 1 )-, • • •) 

A coadjoint orbit consists of elements E k with a fixed maximal order of the pole. Then, 
we can interpret eq.(1.7) as the coadjoint orbit of the element ((Al)_, (A2)-, ■ ■ ■)■ 

Alternatively, we can consider the function on Q* 

L{\) = L G + Y J ^k (1.19) 

k 

with poles at the points e k . Given this function we can recover the E k by extracting the 
polar parts. The constant matrix Lq is added to match the formula for the Lax matrix 
eq.(1.7). By choice it is assumed to be invariant by coadjoint action. The pairing can 
be rewritten as 

(L,X) = ^TrRes efc L(A)A,(A) 

k 

Remark that only E k contributes to the residue at e k and the formula is compatible with 
the matrix Lq being invariant by coadjoint action. 

1.5 Classical r-matrix. 

We can now use this symplectic form to evaluate the Poisson brackets of the elements of 
the Lax matrix and show that they take the r-matrix form. Let us first introduce some 
notations. Let be the canonical basis of the N x N matrices, (Eij) k [ = 5i k Sji. We 
can write 

L(\) = Y,L ij (\)E ij 

ij 

Let 

Li(A) = L(A) ® 1 = J2 L ij( X )( E iJ ® 1), L 2^) = 1 ® L(n) = ^Ly(/x)(l E i:j ) 

ij ij 

The index 1 or 2 means that the matrix L sits in the first or second factor in the tensor 
product. More generally when we have tensor products with more copies, we denote by 
L a the embedding of L in the a position, e.g. L3 = 1<8>1<8>£<8>1<8>---. Finally, we define 
{Li(A), L2(fi)} as the matrix of Poisson brackets between the elements of L: 

{Li(A), L 2 (fi)} = ^2{Lij(X), LubifiEij ® E kl 

ij,kl 
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We assume that each L k (X) is a generic element of an orbit of the loop group 
GL(N)[X] that Lq and the are non-dynamical. Each orbit L k (X) is equipped with 
the Kostant-Kirillov Poisson bracket and 



{L k (X),L k >(Y)} = 0, k + k' 



(1.20) 



Proposition 3 With these assumptions, the Poisson brackets of the matrix elements of 
L(X) can be written as: 




{Li(X) , L 2 (fJ,)} 

A - I! 

(1.21) 

with C\2 = Ylij Eij®Eji where the E^ are the canonical basis matrices. The commutator 
in the right hand side of eq.(1.21) is the usual matrix commutator. 

Proof . Let us first assume that we have only one pole and L = (gA^g^ 1 )-. Because 
we are dealing with a Kostant-Kirillov bracket for the loop algebra of gl(N), we can 
immediately write the Poisson bracket of the Lax matrix using the defining relation 
{L(X),L(Y)} = L([X,Y]). Using, L(X) = Tr Res A =o(£(A)X(A)), this gives: 

{L(X), L(Y)} = Tr Res A=0 (L(X)[X(X),Y(X)}) (1.22) 

By definition of the notation {L\, L 2 }, we have: 

{L(X),L(Y)} = <{L!(A), L 2 (/x)} , X(X) ® F(/x)) 

where (, ) = Tri 2 ResARes^. We need to factorize X(X) <8> Y([i) in eq.(1.22). To this end, 
we introduce the operator, assuming |A| < |/x|,: 



A'' 



Ci 2 (A,M) = C 12 ^— t 



C\2 

A — /i' 



C 



12 



E 



3 « 



n=0 ' i,j 

This operator is such that for Y(X) = Yl^Lo ^"A™ we have 

Fi(A) =Tr 2 Res„ C 12 (A, n)Y 2 {n) 

We can now write: 



(L(A)pT(A), Y(A)]) = <[C 12 (A, »), L(A) ® 1] , X(A) >») 

Consider the rational function of A: <p(X) = {Li(A), L 2 (/i)} — [Ci 2 (A, /x), L(A) <8) 1]. By 
inspection (/? contains only negative powers of (i, and we have {<p,X(X) ® Y(/u)) = 0. 
Hence <^ contains only positive powers of A and is regular at A = 0. It has a pole at 
A = fx, due to the form of C(X, (i). We remove this pole by subtracting to (p the quantity 
[Ci 2 (A, fx), l<g)L(//)] which contains only positive powers of A and is therefore in the kernel 
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of (■, X(X) <8> Y(fi)}. The pole at A = fi disappears since [C12, L(fi) <g) 1 + 1 <g) L(fi)] = 0. 
The redefined 99 is regular everywhere and vanishes for A — > 00 hence vanishes identically. 
This proves eq.(1.21) in the one-pole case. 

We can now study the multi-pole situation occuring in eq.(1.7). Consider L = 
Lo + J2k Lk- Each Lk lives in a coadjoint orbit as above equipped with its own symplectic 
structure. From eq.(1.20) they have vanishing mutual Poisson brackets 

{(L fe (A))i,(L fc ,(M))2} = 0, k + k' 

We assume further that Lq does not contain dynamical variables 

{(£0)1, (£0)2} = 0, {(L )i, GMA)) 2 } = 

Then since Ci2/(A — fx) is independent of the pole e&, it is obvious that the r-matrix re- 
lations for each orbit combine by addition to give eq.(1.21) for the complete Lax matrix 
L(A). ■ 



We have obtained a very simple formula for the r-matrix specifying the Poisson 
bracket of L(X): 



ri2(X,n) = -r 2 i(/x,A) 



C 



12 



(A-/x) 



1.23) 



The Jacobi identity is satisfied because this r-matrix verifies the classical Yang-Baxter 
equation 



[ri2,ri 3 ] + [ri2,r 23 ] + [ri 3 ,r 23 ] = 



where r^ stands for rjj(Aj, Xj). Note that n 2 is antisymmetric: ri 2 (Ai, A 2 ) = — r 2 i(A 2 , Ai). 

These Poisson brackets for the Lax matrix ensure that one can define commuting 
quantities. The associated equations of motion take the Lax form. 

Proposition 4 The functions on phase space: 



#(")(A) = Tr (L n (A)) 



are in involution. The equations of motion associated to H^ n \/j,) can be written in the 
Lax form with M = M^: 




(1.24) 
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Proof . The quantities H( n \\) are in involution because 

{Tr L n {\) , Tr L m (/x)} = nmTi l2 {L 1 {\),L 2 { l i)}L n l -\\)L^- 1 ^) 

-Tr 12 ([C^LUm™- 1 ^) + [a2,L?(ri]Ll-\\)) = 



A — ji 

where we have used that the trace of a commutator vanishes. Similarly, we have: 



L(\) = {HW(ri,L(X)} = nTr 2 
Performing the trace and remembering that Tr2 (C12M2) = Ml, we get 



A — 



L(A) = [MW(A,/x),L(A)], MW(A,/,) = n^M (1 . 25 ) 

This M^(A, fi) has a pole at A = /i and is otherwise regular. According to the general 
procedure we can remove this pole by subtracting some polynomial in L(X) without 
changing the equations of motion. Obviously one can redefine: 

MH(A, M ) - M(")(A,/x) -n^-^ = -n L "' 1(A x } ~ L "' lM (1.26) 

This new M has poles at all e& and is regular at A = /x. Decomposing it into its polar 
parts, we write M = with 

M k (X) = -n < Ln ~ 1{X) 

This is of the form eq.(1.13) with 



p( fc )(L,A) = -— L n-1 (A) (1.27) 
A — /x 

Notice that the coefficients of the polynomial P^ k \h, A) are pure numerical constants. 

This proposition shows that the generic Zakharov-Shabat system, equipped with 
this symplectic structure, is an integrable Hamiltonian system (the precise counting of 
independent conserved quantities will be done in Chapter [2]). 

1.6 Examples. 

1.6.1 The Jaynes-Cummings-Gaudin model. 

We consider the following Hamiltonian 

n—l n—l 

H = ^ lerf +ubb + gY^ (bsj + (1.28) 

j=0 j=0 
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The Sj are spins variables, and b, b is a harmonic oscillator. The Poisson brackets read 

{s°,s b j } = -e abc s c j , {b,b} = i (1.29) 
The Sj brackets are degenerate. We fix the value of the Casimir functions 

% ' s j = s 

Phase space has dimension 2(n + 1). The equations of motion read 

n-l 

b = —iujb — ig~^2 sj (1.30) 

3=0 

s* = ig(bsj-bs+) (1.31) 

s+ = 2ie jS + - ligbs] (1.32) 

sj = -2ie jS ~ + 2z 5 6s| (1.33) 



We introduce the Lax matrices 



L(A) = Aa* + -(fcr + + &0- V + E flf ( L34 ) 

y y y J=0 ■? 

M(A) = -i\a z -ig(ba + + ba-) (1.35) 
where a a are the Pauli matrices. 

a ± = ^(a x ±ia y ), [a z , = ±2a ± , [a + , a~] = a z 
It is not difficult to check that the equations of motion are equivalent to the Lax equation 

L(\) = [M(\),L(\)\ (1.36) 



Let 



we have 



L(X) _(A(X) B(\) \ 
m -{c(X) -A(X)J 



w £-1 a? 



9 A — e,- 
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It is very simple to check that 

{A(\),A(ri} 
{B(\),B(ri} 
{C(\),C(l*)} 
{A{\),B(n)} 

{A(X),C(»)} 

{B(X),C(»)} 



A — n 



(B(X) - B(»)) 
-(C(X)-C^)) 



A — fi 

One can rewrite these equations in the usual form 

" Pl2 



X — 

2> (A(X) - A(fi)) 



{L 1 (X),L 2 ( l i)} = -i 



X — n 



,L 1 (X) + L 2 (fi) 



where 



12 



It follows that 



(\ 0\ 

10 

10 

\0 1/ 



;TV(L 2 (A)) =A 2 (A)+B(A)C(A) 



generates Poisson commuting quantities. One has 



2 ^ H; 



n—l 



W(A) = l(2A-.) 2 + ^ + ^^ + £ ( ^ 



2 • <T 9 

where the (n + 1) Hamiltonians read 

H n = bb + s , 

3 

and 



a* * — ' A — 6j 



3=0 



Hj = (26, - u,) S | + g(bsj + bsj) + 



Sj • Sfc 
• e 3 ~ £ fc 



(1.37) 



j = 0, • • ■ , n - 1 



The Hamiltonian eq.(1.28) is 

n— 1 

H = coH n + Y J Hj 
j=o 

Let us see how these formulae fit into our general scheme. The Lax matrix is a sum 
of simple poles at 6j . The loop group at each one of these points is 

g®(\) = gg ) + (\-e j )g? ) + - 
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where gjp are SU{2) matrices. In fact the pole being simple, only contributes to 
the coadjoint orbit. 

HX) = (>(A) ^ ^(A)) = ^ g?*^- 1 

At infinity we consider the loop group 

5 M (A) = 1 + A -i/-) + ... 

and the coadjoint orbit 

Loo(X) = \ \a z g^'\\)) = ^A^ + ^bSVl = + 2 {ba + + ba~) 

g z \ J+ g z g z 9 9 

Identifying 

g 2 

we do have 

71-1 

L(A) = L + L oo (X) + Y J L j W 

j=0 

To see how M(A) also fits into the scheme, we consider separately the evolution with 
respect to the Hamiltonians Hj, H n . Since Hj is the coefficient of g i^_ e ^ i n ^Tr L 2 (/j,) 
we just have to extract this coefficient in eq.(1.26), which for n = 2 reads (there is an 
extra factor i coming from the definition of the Poisson bracket) 

A — /i 

We find 

,2 



g s 3 • g 

jV ' v '2 A < 



Mj(A) = i- 



; y 



Similarly, for coH n we have to extract the term in /i u in the same expression. We find 

Moo(A) = -i^a z (1.38) 



Hence 



2 / n ~i -> -> 

M(A) = M 00 (A) + M,(A) = i^- "o* + X;^if 

V 5 i=o J 



2 

-i(Acj 2 + </(6<t+ + to")) + *y L (A) 



Since the L(A) term does not contribute to the Lax equation, we have recovered the 
expression eq.(1.35) for the matrix M(A). 
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1.6.2 The KdV hierarchy. 

In the one pole case, we have seen that the general structure of a Lax equation is 



L(A) 



(p(L,\))_,L(\) 



(1.39) 



The KdV hierarchies are constructed exactly on the same pattern but replacing loop 
algebras by the algebra of pseudo differential operators which we now describe. 



The algebra of pseudo-differential operators is the algebra of elements of the form 

N 

A= Yj ^ 

i=— oo 

with N finite but arbitrary. The coefficients a» are functions of x, dis the usual derivation 
with respect to x and the "integration" symbol, <9 _1 is defined by the following algebraic 
rules: 

d^d = dd- 1 = 1 

oo 

a _1 a = J2(-l) i (d i a)d- i - 1 (1.40) 

i=0 



We denote hy V = = S-oo j tne se * °f formal pseudo-differential operators. 
Let V + = [a = zf= 

aj3*| be the subalgebra of differential operators, and let V- - 

= ^I^ajC^j be the subalgebra of integral operators. We have the direct sum 
decomposition of V as a vector space: 

V = V+ © V- 

Notice that V is naturally a Lie algebra. P + and V~ are Lie subalgebras. but V+ and 
V- do not commute. 

The formal group G = exp(V ) is called the Volterra group. We have G ~ 1 + V- 
because powers of elements in V- are in V- . Let be an element of G: 

oo 

$ = 1 + Y^w i d~ i G (1+V-) (1.41) 
i=l 

The coefficients of its inverse < & _1 = 1 + Y1T can be computed recursively from 

the relation = 1. 
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For A G V we define its residue, denoted Resg A, as the coefficient of d in ^4: 



Resg ^4 = a-i(x) 

1 (1.42) 

On V there exists a natural linear form called the Adler trace, denoted ( ), defined by: 



{A) = J dxResoA = J dx a_i(x) 



(1.43) 

This linear form satisfies the fundamental trace property {AB) = {BA). 

Returning to integrable systems, we now let L be a differential operator of order n 



n-2 



L = d n ~Y J Uid i 



(1.44) 



In the algebra of pseudo differential operators its n-th root exists 

Q = L^, Q = d + q- 1 d- 1 + ••• 

The generalized KdV hierarchies are defined by the Lax equations (compare with 
eq.(1.39). In both cases the projection is on the dual of the Lie algebras Q[X] and V- 
respectively.) 



L 



(1.45) 



These equations are consistent for all k G ]N in the sense that we have a differential 
operator of order n — 2 on both side. To see it, we notice that Q k , V/c G N commutes 
with L since LQ k = Q n+k = Q k L. Then, we have: 



Q 



Q ,L 



Q 



L 



Q 



is of order less 



From the last equality, it follows that the differential operator (Q k ) + , L 

or equal to n — 2, so that the Lax equation eq.(1.45) is an equation on the coefficients 
of L. This is the original Gelfand-Dickey argument. 

The differential operator L is an element of V+. If we view V+ as the dual of the 
Lie algebra V- through the Adler trace, there is a natural Poisson bracket on V+: the 
Kostant-Kirillov bracket. For any functions / and g on V+, it is defined as usual by: 



{f,g}(L) = (L,[df,dg}} V L e V + 



(1.46) 
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where we understand that df, dg G V— In particular, for any X = Yl'jLo ® lx j e ^-j 
we define the linear function fx(L) by: 

f x (L) = {L,X) (1.47) 

we have df x =Xe V- Therefore {/*, M = <L, [X, Y]> = f [xy] for any IJe P_. 

Proposition 5 Let L G "P+ 6e i/te differential operator of order n as in eq.(1.44)- Define 
the functions of L by 




(i) The functions H k (L) are the Hamiltonians of the generalized KdV flows under the 
bracket eq. (1.4-6): 



d tk L = {H k ,L} 



L~> 



L 



(1.48) 



(ii) The functions H k {L) are in involution with respect to this bracket. 



Proof . We first need to compute the differential of the Hamiltonian H k . Let L and 5L 
be differential operators of the form eq.(1.44). One has, using the cyclicity of Adler's 
trace: 

((L + SLY) = (L u ) + v{L u ~ 1 5L) + • • • 

which implies d{L u ) = u(L u ^ 1 )^ n ^ where the notation ( )(_ n ) means projection on V- 
truncated at the first n — 1 terms. This projection appears because 5L = —5u n ~2d n ~ 2 — 
• • • — 5uo which is dual to elements of the form 8~ 1 xq + • • • + d~ n+1 x n -\ under the Adler 
trace. Hence: 



dH k (L) = [Ln 



(-n) 



Q 



(-n) 



G V- 



We call 6 



(k) 
» 



the terms left over in the truncation: 



k 

L n 



dH k + e 



-(n) 



(1.49) 



(1.50) 



We now prove eq.(1.48). Consider the function fx{L) = {LX), then 

fx = (L,X) = {H k ,f x }(L) = (L,[dH k ,dfx}} = ([L,dH k ],X) 

where we used the invariance of the Adler trace. Since X G V-, only [L,dH k ] + con- 
tributes to this expression. But 



[L,dH k \. 



L, (L« 



L. 



,(*) 



where we have used = 0, and the fact that [L,0^ n ^] + = 0. So [L,dH k ] + is 

differential operator of order at most n — 2, and this proves eq.(1.48). 
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Next we show that the Hamiltonians H k are in involution. We have: 
{H k , H k ,}{L) = (L, [dH k , dH k ,)) = ([L, dH k } + , dH k ,) 



< 



( L ") + ' 



, dH k i) 



Using again the fact that [L,dH k ] + is of order at most n — 2, we can replace dH k > by 



L 



and get: 

{H k , H k i}{L) = ( 



L~< 



L 



k' 



) = ( 



Li 



,L 



Ln) 



(fc' \ k' 
L~^\ by L^. Finally, 

from the invariance of the trace, we obtain: 

{H k ,H k ,}(L) = <(l£)_ 



+ L 



L, L" 



) = 



These systems are called the generalized KdV hierarchies. They are field equations 
or infinite dimensional systems. They are integrable in the sense that they possess an 
infinite number of Poisson commuting conserved quantities but we are already beyond 
the strict framework of the Liouville theorem. 

The KdV hierarchy corresponds to n = 2 and the generalized ones to n = 3,4, • • -. 
Let us consider the KdV case n = 2. The operator L is the second order differential 
operator 

L = d 2 — u 

We first find Q such that Q 2 = L. One has Q 2 = d 2 + 2g_i + (2g_ 2 + dq-^d' 1 + • • • so 
that q-\ = — \u, q~2 = \du, etc... 

Q = 0-i«0- 1 + i(0 u )0- 2 + ... 

We again check on this simple example that all the q-j are recursively determined in 
terms of u by requiring that no d~i terms occur in Q 2 . To obtain the KdV flows, we only 
have to compute (Q k ) + , k = 1,2,- ■ -. For k = 1, we have (Q)+ = d, and d±L = [d,L]. 
This reduces to the identification <9 tl = 5. For k = 2, we have (Q 2 )+ = L and we get 
the trivial equation dt 2 L = 0. The first non trivial case is k = 3. We have 

(Q 3 ) + = d 3 - 3 -ud- 3 -(du) 

so the Lax equation reads dt 3 u = [(Q 3 ) + , d 2 — u] which is the Korteweg-de Vries equation: 

4<9 tg u = d 3 u — 6u(du) 



This is the first of a hierarchy of equations obtained by taking k = 3, 5, 7, ■ ■ ■ called the 
KdV hierarchy (note that for k even we get trivial equations). 
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Chapter 2 

Solution by analytical methods 



We present the general ideas for the solution of Lax equations when a spectral parameter 
is present. The method uses the geometry of the spectral curve and its complex analysis. 



2.1 The spectral curve. 



Let us consider a N x N Lax matrix L(\), depending rationally on a spectral parameter 
A € C with poles at points 



L(A) =L Q + Y,L k {\) 



(2.1) 



As before Lq a constant matrix is independent of A and Lfc(A) is the polar part of L(A) 
at €k, i.e. 

-i 

£fc = ^2 L ktr (X - ek) r 

r=-n h 

The analytical method of solution of integrable systems is based on the study of the 
eigenvector equation: 

(L(A)-/xl)*(A,Ai) = 

(2.2) 

where ^(A, fi) is the eigenvector with eigenvalue fi. The characteristic equation for the 
eigenvalue problem (2.2) is: 



r : r(A,/i) = det(L(A)-/il) = 

(2.3) 

This defines an algebraic curve in C 2 which is called the spectral curve. A point on T 
is a pair (A, fi) satisfying eq.(2.3). Since the Lax equation L = [M,L] is isopsectral, this 
curve is independent of time. 
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If N is the dimension of the Lax matrix, the equation of the curve is of the form: 

JV-1 

T : T(X^)^(- f ,) N +J2r q (X)^ = (2.4) 

q=0 

The coefficients r q (X) are polynomials in the matrix elements of L(X) and therefore have 
poles at €k- The coefficients of these rational functions are independent of time. 

From eq.(2.4), we see that the spectral curve appears as an iV-sheeted covering of the 
Riemann sphere. To a given point A on the Riemann sphere there correspond N points 
on the curve whose coordinates are (A, fi±), • • • (A, /jln) where the \Xi are the solutions of 
the algebraic equation T(X,fi) = 0. By definition fa are the eigenvalues of L(X). 

Our goal is to determine the analytical properties of the eigenvector *(A, fi) and 
see how much of L(X) can be reconstructed from them. The result is that one can 
reconstruct L(X) up to global (independent of A) similarity transformations. This is not 
too surprising since the analytical properties of L(X) and the spectral curve are invariant 
under global gauge transformations consisting in similarity transformations by constant 
invertible matrices. So from analyticity we can only hope to recover the system where 
global gauge transformations have been factored away. 

In general, we may fix the gauge by diagonalizing L(X) for one value of A. To be 
specific, we choose to diagonalize at A = oo, i.e. we diagonalize the coefficient Lq. 

Lq = lim L(X) = diag(ai, ■ ■ ■ , ajv) (2-5) 

A^oo 

We assume for simplicity that all the Oj's are different. Then on the spectral curve, we 
have N points above A = oo: 

Qi = (A = oo,m = en) 

In the gauge (2.5) there remains a residual action which consists in conjugating the Lax 
matrix by constant diagonal matrices. Generically, these transformations form a group 
of dimension N — 1 and we will have to factor it out. 

2.2 Riemann surfaces. 
2.2.1 Desingularisation 

A Riemann surface is a compact analytic variety of dimension one. This means that 
there is a covering by open neighborhoods Ui and local homeomorphisms mapping them 
to the open disks \zi\ < 1. On the intersection Ui n Uj the local parameters z% and Zj 
are related by an analytic bijection. 

We will deal however with curves in C 2 as in eq.(2.4). To relate it to the abstract 
definition we have to find around any point of V a local analytic parameter z. 

A point Ao, hq is regular if none of the derivatives of T(A, /i) vanish at that point. 

d x T(X o ,HQ)^0, ^r(A ,/i )/0 
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Around such a point we can choose the local parameter z as either A — Ao or li — liq and 
use the equation T(A, /x) = to express the other one analytically in terms of z. 

A point Ao,/xo is a branch point if one of the derivatives of T(A,/i) vanish at that 
point, but not both. Let us assume 

d x T(X o ,n )^0, ^r(A o ,Mo) = 

Around such a point the curve looks like (for a branch point of order 2) 

r(A, $ ~ a A r(A , /i )(A - a ) + ^r(A , ^ - ^? + ■■■ 



Clearly, one cannot use A — Ao as a local parameter because then fi — liq ~ ay 7 A — Ao + - • • 
is not analytic at A = Ao- However choosing z = \i — liq is perfectly legal. 

A point Ao, fJ*o is a singular point if both derivatives of T(A, fx) vanish at that point 

d x T(X ,fi ) = 0, ^r(Ao,Mo) = 

An example of such a point is given by the curve 

^ = aX 2 + bX 3 

To give a meaning to the singular point we perform a birational transformation 

A = z , 11 = zy (2.6) 

which can be inverted rationaly as long as (A, fi) / (0, 0) 

z = X, y = - 

Under the transformation eq.(2.6) the curve becomes 

y 2 = a + bz (2.7) 

Now, instead of the singular point (0,0), we get two regular points (z = 0,y = ±\/a) 
which are mapped to the singuar point (0, 0) by eq.(2.6). We say that we have resolved 
(or blown up) the singularity by performing the birational transformation. It is always 
the desingularized curve like eq.(2.7) that we must consider to give a meaning to a 
singular point. We can use z as a local parameter on the desingularized curve. 

Finally a special care must be taken for the points at oo where A and (or) fi become 
infinite. Then we set 

X = 1/x, fi = l/y (2.8) 

this brings it to the origin (0, 0). In general we get a singular point and we have to blow 
it up with the above procedure. 

Important examples are the cases of hyperelliptic curves. They are defined by 

^ 2 = i 3 2n+ 2 (A), or ^ = P 2n+1 (A) (2.9) 
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where f*2n+2(A) and P2n+i(A) are polynomials of degree 2n + 2 and 2n + 1 respectively. 

i 3 2n+ 2 (A) = a\ 2n+2 + b\ 2n+1 + • • • , or P 2n+ i(A) = aA 2n+1 + b\ 2n + • • • 
To understand the point at oo, we perform the transformation eq.(2.8). We get 

y 2 (a + bx^ ) = x 2n+2 , or y 2 (a + bx^ ) 

The point (x, y) = (0, 0) is now a singular point. To blow it up we set 



„2n+l 



X = X 

y = x >n+l y > OT 



X = X 

V = x' n y' 



we obtain 



y 



12 



1 

(a + bx' + 



or y 



12 



X 

(a + bx' + 



In the first case the singularity has been resolved into two regular points, the local 
parameter can be taken to be z = x' . In the second case it has been resolved into a 
branch point, the local parameter can be taken to be z = y' . Summarizing, at infinity 
we have 



U _| 



or 



( X 



l l 



+ 



a n z 2n + l ' 



Notice that the number of branch points in both cases is 2n + 2. 



2.2.2 Riemann-Hurwitz theorem. 

This is the tool to compute the genus of a Riemann surface. Given a triangulation of 
the surface, the Euler-Poincare characteristic is defined by 



X = F - A + V 



where F is the number of faces of the triangulation, A the number of edges and V the 
number of vertices. The Euler-Poincare characteristic is a topological invariant. It is 
related to the genus by the formula 



X = 2 - 2g 



For instance in the case of the sphere, there is a triangulation with 8 faces, 12 edges and 
6 vertices. Hence 

Xo = 2, go = 

If r is a branched covering of Tq, we can lift a triangulation of To to T. Let us choose the 
triangulation of Tq such that the projections of the branch points are among its vertices. 
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Let Fq, Aq, Vq be the number of faces, edges and vertices of this triangulation. Let N 
be the number of sheets of the covering r — > To- When we lift the triangulation of To 
to T, we get a triangulation with F = NFq faces, A = NAq edges and V = NVq — B 
vertices where B is the total index of the branch points. At a branch point the index is 
the number of sheets that coalesce minus one. Hence we find 

X = N X o - B 

or 

2g - 2 = iV(2 50 - 2) + B 



This is the Riemann-Hurwitz formula. Let us apply it to compute the genus of the 
hyperelliptic curves given by eqs.(2.9). They are coverings of the Riemann sphere with 
two sheets. Each branch point is of index 1. In both cases we have seen that there are 
exactly 2n + 2 branch points. Therefore 2g — 2 = 2(2 x — 2) + 2n + 2, that is g = n. 

2.2.3 Riemann-Roch theorem. 

This is the tool to count the number of meromorphic functions on a Riemann surface. 
A Divisor is a formal sum of points with multipli cites. 

D = n\P\ + 112P2 H h n r P r , eZ 

The degree of the divisor is 

deg D = ^2 n i 

i 

If / is a meromorphic function, we denote by (/) its divisor of zeroes and poles. Let 
M. (D) the space of meromorphic functions whose is such that 

(f)>D 

that is M(D) is the det of meromorphic funtions whose order of the poles are at most 
the one specified by D and the order of zeroes are at least the one specified by D. 
The Riemann-Roch theorem asserts that 



dim M(-D) = i(D) + &egD-g + l 



where i{D) is the dimension of the space of meromorphic differentials to such that 

M > D 

There are two cases where the theorem leads to simple answer. If D = 0, then M.{D) 
is the set of holomorphic functions on the Riemann surface. We know that the only such 
function is the constant. Hence dim M{D) = 1. Similarly i(D) is the dimension of the 
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space of holomorphic differentials. The theorem then give i(D) = g. We therefore get 
the very important result 



dim (holomorphic differentials) = g 



If D < 0, then dim M(D) = and this allows to count the meromorphic differentials 



i(D) = -deg D + g - 1, deg D < 



Notice that it is the degree of the divisor which is relevant. The freedom gained by 
adding a pole is compensated by the restriction of adding a zero. 

The next simple case is when deg D > g where generically i(D) = 0. For instance if 

D = Pi + P 2 + ■ • ■ + Pg 

then i{D) is the dimension of the space of holomorphic differentials with zeroes at the 
point of D. To construct such differentials we expand them on a basis uj-i of the g 
holomorphic differentials and try to impose the conditions 

^2ui(Pj)ci = 0, j = l,---,g. 

i 

This system in general has no solution because for a generic set of g points Pj we have 
detu)i(Pj) / 0. Hence i(D) = 0. Then the theorem gives 



dim M{— D) = deg D — g + 1, deg D > g 



The difficult case is when < deg D < g and a careful investigation is needed. 
2.2.4 Jacobi variety and Theta functions. 

Consider a Riemann surface T of genus g. Let ai,bi be a basis of cycles on T with 
canonical intersection matrix (a^ • a,j) = (bi ■ bj) = 0, (ai ■ bj) = 5ij. One can continuously 
deform these loops without changing the intersection index which is the sum of signs ±1 
at each intersection according to the orientation of the tangent vectors. In particular 
one can deform the loops a« and bi so that they have a common base point and then 
cut the Riemann surface along them. We get a polygon with some edges identified. The 
boundary of this polygon can be described as a± ■ b\ ■ a^ 1 ■ b^ 1 ■ ■ ■ a g ■ b g ■ a" 1 • b~ x where 
the identifications are obvious. The common base point becomes all the vertices of the 
polygon. 

The globally defined analytic one-forms on T are called Abelian differentials of first 
kind. They form a space of dimension g over the complex numbers. There is a natural 
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pairing between these forms and loops obtained by integrating the form along the loop. It 
can be shown that the pairing between a-cycles and differentials is non degenerate (note 
they have the same dimension g). We choose a basis of first kind Abelian differentials, 
which we denote by LOj, j = 1, ■ ■ ■ , g, normalized with respect to the a-cycles : 



Wi = 8i3. (2.10) 

The matrix of 6-periods is then defined as the matrix B with matrix elements : 

Bij = <b ujj (2.11) 



Taking the example of an hyperelliptic surface y 2 = P2 g +i(x) where P (x) is a poly- 
nomial of degree 2g + 1 , a basis of regular Abelian differentials is provided by the forms 

x J_1 dx 

uj = for j = 1, ■ ■ • , g 

V 

These forms are regular except perhaps at the branch points and at oo. At a branch 
point the local parameter is y and we have y 2 = a(x — b) + • • • hence x 3 ~ 1 dx/y = 
(2b J ~ 1 /a)(l + • • -)dy which is regular. At oo we take x' = 1/x and y' = y/x^ 9+1 ^ so that 
y' 2 = ax' + ■ ■ ■ and x^ 1 dx/y = by' 2<y9 ~^ dy' which is regular for 1 < j < g since y' is the 
local parameter. Of course these forms are unnormalized. 

Similarly Abelian differentials of second kind are meromorphic differentials with poles 
of order greater than 2. Given a point p on T, there exists a unique normalized (all a- 
periods vanish) Abelian differential of second kind whose only singularity is a pole of 
second order at p. Indeed applying the Riemann-Roch theorem with deg D = — 2, we 
find i(D) = o + l. The g comes from the first kind differentials which are included in this 
counting. Adding a proper combination of differentials of first kind one can always insure 
that all a-periods of the second kind differential vanish and the differential becomes 
uniquely determined. 

We define Abelian differentials of third kind as general meromorphic differentials with 
first order poles whose sum of residues vanish (this condition results from the Cauchy 
theorem). Given two points p and q there exists a unique normalized (all a-periods 
vanish) third kind differential whose only singularities are a pole of order 1 at p with 
residue 1, and a pole of order 1 at q with residue -1. 

On a Riemann surface on which we have chosen canonical cycles, there is a pairing 
between meromorphic differentials. Namely, let fii and H2 be two meromorphic differ- 
entials on S. The pairing (fli • f^) is defined by integrating them along the canonical 
cycles as follows: 

j> 0,i j> SI2 — j> ^2 <j> ^1 

V". bj dj bj J 



E 

i=i 



37 



The Riemann bilinear identity expresses this quantity in terms of residues: 

Proposition 6 Let g\ be a function defined on the Riemann surface, cut along the 
canonical cycles, and such that dg\ = fii. We have: 

(fti • n 2 ) = 2m res(sifi 2 ) (2.12) 

poles 

Corollary 2 The matrix of b-periods B is symmetric. 

Corollary 3 Let Cl 2 be a normalized differential of second kind with a pole of order 
n, with principal part z~ n dz at z = for some local parameter z. Let £l\ = iO k be a 
normalized holomorphic differential expanded as 

CO 

i=0 

around z = 0. One has: 

Cl 2 = 27ri- n 



By linearity, if is a normalized second kind differential with principal part dP(z) 
where P{z) = J2n=iPn z ~ n > then we have 

— I ft (p) = -Res (to k P) (2.13) 
2^7^ J bk 



Consider a divisor of degree which can always be written T> = Y^iiPi — q%), with 
non necessarily distinct points. Choose paths 7$ from to pi and associate to V the 
point in C 9 of coordinates: 

PkiP) = ^ / u k, k = g 

i T> 

where the tOi are holomorphic differentials. Such sums are called Abel sums. If the paths 
are homotopically deformed these integrals remain constant by the Cauchy theorem. If 
one makes a loop around a k then pi — > p\ + 5 k i- If one makes a loop around b k , then 
Pi - ► Pi + &kl- Hence the maps p k give a well-defined point on the torus: 

J(r) = C 9 / (Z 9 + B7L 9 ) (2.14) 

where £> is the matrix of the 6-periods. If one permutes independently the points pi and 
qi the point in the torus does not change. To see it, let the paths j[ connect q\ to p 2 , 
72 connect q 2 to p\ and a connect q\ to q 2 . One has u = J^, u — u up to periods 

and J^ 2 to = f^,uj + J a lo up to periods, so w cancels in the sum. 

The theorems of Abel and Jacobi state that the point on the torus J(T) characterizes 
the divisor V up to equivalence. 
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Theorem 4 (Abel) A divisor V = J2i(Pi~Qi) ^ s the divisor of a meromorphic function 
if and only if, for any first kind Abelian differential uj, the Abel sum J2i f T w vanishes 
modulo Z + BTL for any choice of paths 7, from qi to pi. 

Theorem 5 (Jacobi) For any point u G J(T) and a fixed reference divisor Vq = 
X)f=i Qi> one can fi n d a divisor of g points V = p\ + • • • + p g on £ such that Pk{L) — A)) 
maps to u. Moreover for generic u the divisor V is unique. 

One can embed the Riemann surface V into its Jacobian J(T) by the Abel map. 

Namely, choose a point qo £ T and define the vector A(p) with coordinates Ak(p) 
modulo the lattice of periods: 

A:T 1 — ► J(r) (2.15) 

A k (p) = [ P io k (2.16) 

Clearly, the Abel map depends on the point qo. But changing this point just amounts 
to a translation in J(T). 

One can show using Riemann bilinear type identities that the imaginary part of the 
period matrix B is a positive definite quadratic form. This allows to define the Riemann 
theta-function: 

9(zi, ...,z g )= e 27ri (™> 2 )+™( Sm > m ). (2.17) 

Since the series is convergent, it defines an analytic function on C 9 . 

The theta function has simple automorphy properties with respect to the period 
lattice of the Riemann surface: for any / G Z 9 and z € C 9 

0(z + l) = 9{z) 
9(z + Bl) = exp[-iir(Bl,l)-2iir(l,z)]9(z) (2.18) 

The divisor of the theta function is the set of points in the Jacobian torus where 9{z) = 0. 
Note that this is an analytic subvariety of dimension g — 1 of the torus, well-defined due 
to the automorphy property. 

The fundamental theorem of Riemann expresses the intersection of the image of the 
embedding of T into J(T) with the divisor of the theta function. 

Theorem 6 (Riemann) Let w = (w±, ■ ■ ■ , w g ) G C 9 arbitrary. Either the function 
9(A(p) — w) vanishes identically for p G T or it has exactly g zeroes p±, ■ ■ ■ ,p g such that: 

A(pi) + ■ ■ ■ + A(p g ) = w - K (2.19) 

where K, is the so-called vector of Riemann's constants, depending on the curve V and 
the point qo but independent of w. 
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2.3 Genus of the spectral curve. 



Before doing complex analysis on T, one has to determine its genus. A general strategy 
is as follows. As we have seen, T is a iV-sheeted covering of the Riemann sphere. There is 
a general formula expressing the genus g of an iV-sheeted covering of a Riemann surface 
of genus go (in our case go = 0). It is the Riemann-Hurwitz formula: 

2g - 2 = N (2g - 2) + B (2.20) 

where B is the branching index of the covering. Let us assume for simplicity that the 
branch points are all of order 2. To compute B we observe that this is the number of 
values of A where T(A, fi) has a double root in /x. This is also the number of zeroes of 
dfjT(X, n) on the surface T(A, ji) = 0. But d^T(X, [/,) is a meromorphic function on T, and 
therefore the number of its zeroes is equal to the number of its poles and it is enough 
to count the poles. These poles can only be located where the matrix L(A) itself has a 
pole. So we are down to a local analysis around the points of V such that L(A) has a 
pole. Around such a point the curve reads 

(^"(A^ + ---)---(^(A^)^ + ---) =0 

where lj are the eigenvalues of nfc that are assumed all distinct. When A tends to 
efc, n tends to infinity. We bring this point to the origin by setting 

fi=-, A - e k = z 

y 

Around the point (0,0) the curve reads 

{hy - z n * + ■■■)■■ ■ (l N y ~ z n " + ■ ■ •) = 

Clearly the point (0, 0) is a singular point. To desingularise the curve, we set y = z nk y' 
and we find 

(hy' - 1 + •••)•• • (W - 1 + •••)= 

The singular point has blown up to the N points z = 0,y' = lj 1 . Hence, above a pole 
efc , we have N branches of the form 

h 

H = ^ + ■ ■ ■ > X- ek = z 

On such a branch we have dfjT(X, /u)|(a,^-(a)) = Tli^j (VjW ~ which thus has a 

pole of order (N — l)rifc. Summing on all branches the total order of the poles over is 
Af(Af — l)nfc. Summing on all poles of L(X) we see that the total branching index is 
B = N(N - 1) J2 k n k- Th is gives: 



= — ^ J>fc-iV + l 

k 
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2.4 Dimension of the reduced phase space. 



For consistency of the method it is important to observe that the genus is related to 
the dimension of the phase space and to the number of action variables occuring as 
independent parameters in eq.(2.4), which should also be half the dimension of phase 
space. The original phase space M is the coadjoint orbit 

k 

(k) 

Let us compute its dimension. The matrices A_ characterize the orbit and are non- 
dynamical. The dynamical variables are the jets of order {n k — 1) of the 's which gives 
N 2 n k parameters. But L k is invariant under — ► g( k )$® with d^ a jet of diagonal 
matrices of the same order. Hence the dimension of the L k orbit is (N 2 — N)n k , and the 
dimension of the orbit is the even number: 



dimM = (A 2 -N)^n k 

k 



The reduced phase space A4 rc duced is obtained by performing the quotient by the residual 
global diagonal gauge transformations. 

Proposition 7 The reduced phase space A4 rc d U ccd has dimension 2g and there are g 
proper action variables in eq.(2.4)- 

Proof . The residual global gauge transformations act by diagonal matrices as g k — ► dg k , 
or L(A) — ► dL(\)d~ 1 . This preserves the diagonal form of Lq. The orbits of this action 
are of dimension (N — 1), since the identity does not act. The action of this diagonal 
group is Hamiltonian and its generators are given just below. The phase space -M re duced 
is obtained by Hamiltonian reduction by this action. First one fixes the momentum, 
yielding (A — 1) conditions, and then one takes the quotient by the stabilizer of the 
momentum which is here the whole group since it is Abelian. As a result, the dimension 
of the phase space is reduced by 2(A — 1), yielding: 



dim .M rc d UC cd = (N 2 - N) ^ - 2(N - 1) = 2g 

k 



Let us now count the number of independent coefficients in eq.(2.4). It is clear that 
7j(A) is a rational function of A. The value of r?(A) at oo is known since fij — > dj. Note 
that rj(X) is the symmetrical function aj(fi±,- ■ ■ , hn) where fii are the eigenvalues of 
L(X). Above A = e k , they can be written as 

Hk ciP 

N = Yj (\ — \n + regular (2.21) 
n=i 1 k> 
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where all the coefficients c\ , ■ ■ ■ , Cnl are fixed and non-dynamical because they are the 
matrix elements of the diagonal matrices (A^)-, while the regular part is dynamical. 
We see that rj(X) has a pole of order jn^ at A = e&, and so can be expressed on j J^k n k 
parameters, namely the coefficients of all these poles. Summing over j we have altogether 
a pool of ^N(N + 1) rik parameters. They are not all independent however, because 

in eq.(2.21) the coefficients Sn are non dynamical. This implies that the highest 
order terms in rj(X) are fixed and yields Nn k constraints on the coefficients of rj(X) . 
We are left with ^N(N — 1) J2k n k parameters, that is g + N — 1 parameters. 

It remains to take the symplectic quotient by the action of constant diagonal matrices. 
We assume that the system is equipped with the Poisson bracket (1.21). Consider the 
Hamiltonians H n = (1/n) Res A=00 Tr (L n (A)) dX, i.e. the term in 1/A in Tr (L n (A)). 
These are functions of the rj(X) in eq.(2.4). We show that they are the generators of 
the diagonal action. First we have: 

Res A=00 Tr(L"(A))dA = nRes A=00 Tr (L^ 1 ^ L k (X))dX 

k 

= nRes A=00 Tr(L^ 1 L(A))dA (2.22) 
since all -Lfc(A) are of order 1/A at oo. Using the Poisson bracket we get 

{H n , L(n)} = - Res A=00 Tn L™" 1 ® 1 

The term L(X) ® 1 in the commutator does not contribute because the Lq part produces 
a vanishing contribution by cyclicity of the trace and all other terms are of order at least 
1/A 2 . The term 1(g) L(/i) yields — [Lq" 1 , L{y)\ which is the coadjoint action of a diagonal 
matrix on L(/j,). Since Lq is generic, the Lq generate the space of all diagonal matrices, so 
we get exactly N — 1 generators Hi, ■ ■ ■ , i?jv-i- I n the Hamiltonian reduction procedure, 
the H n are the moments of the group action and are to be set to fixed (non-dynamical) 
values. Setting 

Mi(A) = a j + ^ + --- (2.23) 

around the point Qj = (oo, aj), we have H n = Y2j a1 j~ 1 ^j- So, both (by definition) 
and bi are non dynamical. On the functions rj(X) this implies that their expansion at 

infinity starts as rj(X) = + ^ h • • -, with and rj -1 -* non dynamical. Hence 

when the system is properly reduced we are left with exactly g action variables. 

The constraints eqs.(2.21,2.23) can be summarized in an elegant way. Introduce the 
differential 5 with respect to the dynamical moduli. Then our constraints mean that the 
differential 5[idX is regular everywhere on the spectral curve because the coefficients of 
the various poles being non dynamical, they are killed by 5: 



5fi dX = holomorphic 



Cl2 
X — jJL 



,L(X) 1 + 1® L(fi) 



dX 
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Since the space of holomorphic differentials is of dimension g, the right hand side of the 
above equation is spanned by g parameters which shows that the space of dynamical 
moduli is g dimensional. Notice that these action variables are coefficients in the pole 
expansions of the functions rj(X), and thus appear linearly in the equation of T. Hence 
eq.(2.4) can be written in the form 



r : R(X, n) = R (\, + R i( X > ») H 3 = 



(2.24) 



2.5 The eigenvector bundle. 

Let P be a point on the spectral curve. We assume that P = (A, fi) is not a branch point 
so that all eigenvalues of L(A) are distinct and the eigenspace at P is one-dimensional. 
Let ^(P) be an eigenvector, and ipj(P) its N components: 



*(P) 



MP) \ 

^n(P)J 



Since the normalization of the eigenvector ^(A, fi) is arbitrary, one has to make a choice 
before making a statement about its analytical properties. We choose to normalize it 
such that its first component is equal to one, i.e. 

M P ) = !> at any point P G T. 
It is then clear that the ipj(P) depend locally analytically on P. As a matter of fact: 

Proposition 8 With the above normalization, the components of the eigenvectors ^f(P) 
at the point P = (X,fi) are meromorphic functions on the spectral curve T. 

Proof . Let A(A,/z) be the matrix of cofactors of (L(X) — f/,1), which, by definition, is 
such that (L(X) — /ul)A = r(A, fi)l. Therefore at P = (A, fi) £ T, each column of the 
matrix A is proportional to the eigenvector ^f(P). Hence we have 



A{P) 



Aij(A,/i) 



Aij(A,Ai) 

which is a meromorphic function on T. 

In fact the matrix A(P) is a matrix of rank one, since the kernel of (L(X) — /xl) is of 
dimension one. Hence, for P £ T the matrix elements of A(P) are of the form a.i{P)f3j{P) 
and the components of the normalised eigenvector are ipi (P) = ^(p)^[p) = ^(p) ■ We 
thus expect cancellations to occur when we take the ratio of the minors and we cannot 
deduce the number of poles of the normalized eigenvector by simply counting the number 
of zeroes of the first minor. 
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Proposition 9 We say that the vector ^f(P) possesses a pole if one of its components 
has a pole. The number of poles of the normalized vector *(P) is: 



m = g + N - 1 



(2.25) 



Proof . Let us introduce the function W(X) of the complex variable A defined by: 

W{X) = (det§(A)) 2 
where *(A) is the matrix of eigenvectors of L(X) defined as follows: 



*(A) 



/ 1 1 

MPi) MP2) 



1 \ 

MPn) 



(2.26) 



\Vjv(Pi) iPn{P2) ■■■ ^Pn(Pn)' 



where the points Pi are the N points above A. The function W(A) is well-defined as a 
rational function of A on the Riemann sphere since the square of the determinant does 
not depend on the order of the Pj's. It has a double pole where ^f(P) has a simple pole. 
To count its poles, we count its zeroes. First notice that W(X) only vanishes on branch 
points where there are at least two identical columns. Indeed, let Pi = (/Xj,A) be the 
N points above A. Then the *(P?) are the eigenvectors of L(A) corresponding to the 
eigenvalues fii are thus linearly independent when all the fii's are different. Therefore 
W(X) cannot vanish at such a point. The other possibility for the vanishing of W(X) 
would be that the vector ^f(P) itself vanish at some point (all components have a common 
zero at this point), but this is impossible because the first component is always 1. Let 
us assume now that Ao corresponds to a branch point, which is generically of order 2. 
At such a point W(A) has a simple zero. Indeed let z be an analytical parameter on 
the curve around the branch point. The covering projection P — ► A gets expressed as 
A = Ao + X\z 2 + 0(z 3 ). The determinant vanishes to order z, hence W vanishes to order 
z 2 . This is precisely proportional to A — Ao- Hence W(A) has a simple zero for values 
of A corresponding to a branch-point of the covering, therefore m = B/2. Recall that 
from eq.(2.20) the number of branch points is B = 2{N + g — 1). 



We now need to examine the behavior of the eigenvector around A = 00. At the 
points Qi above A = 00, the eigenvectors are proportional to the canonical vectors e^, 
(ei)k = Sik, since L(X = 00) is diagonal, cf eq.(2.5). While this is compatible with the 
normalization ipi(P) = 1 at the point Qi, it is not compatible at the points Qi, i > 2, 
if the proportionality factor remains finite. The situation is described more precisely by 
the following: 

Proposition 10 The k th component i^k{P) of^(P) has a simple pole at Qk and van- 
ishes at Qi for k = 2, 3, ■ ■ • , N. 
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Proof . Around Qk(X = oo, fi = a^), k = 1, • • • , N, the eigenspace of L(X) is spanned by 
a vector of the form Vfc( A) = efc+0(l/A). The first component of Vk is V k x = 5\k+0{\/ X). 
To get the normalized \I> one has to divide Vk by V^. So we get: 



/ 1 \ 

0(1/A) 



*{P)\p~Qi 



*(P)| 



/ 1 \ 

0(1) 



0(A) 
0(1) 



k > 2 



(2.27) 



\0(1/A)/ \0(1)/ 
where O(A) is the announced pole of the k th component of ty(P)\p^Q k . 

The previous proposition shows that fixing the gauge by imposing that L(X) is diag- 
onal at A = oo introduces N — 1 poles at the positions Qi, i = 2, • • • , N. The location of 
these poles is independent of time, and is really part of the choice of the gauge condition. 
These poles do not contain any dynamical information. Only the positions of the other 
g poles have a dynamical significance. Let D be the divisor of these dynamical poles. 
We call it the dynamical divisor. Recall that the vector ^(P) possesses a pole if one 
of its components has a pole. Therefore the two previous propositions tell us that the 
divisor of the k th components of the eigenvector *(P) is bigger than (-D + Qi — Qk)- 
This information is enough to reconstruct the eigenvectors and the Lax matrix. 

Proposition 11 Let D be a generic divisor on V of degree g. Up to normalization, 
there is a unique meromorphic function ipk{P) with divisor (ip^) > —D + Q\ — Qk- 

Proof . This is a direct application of the Riemann-Roch theorem, since tpk is required 
to have g + 1 poles and one prescribed zero. Hence it is generically unique apart from 
multiplication by a constant ipk ~^ dkipk- 

Equipped with these functions ipk(P) for k = 2, ■ ■ ■ , N we construct a vector function 
with values in C^: 

' MP) 



*(p) 



Wtv(P)/ 



Consider the matrix \£(A) whose columns are the vectors ip(Pi) with Pj = (A,//j) 
are the points above A, cf. eq.(2.26). This matrix depends on the ordering of the 
columns, i.e., on the ordering of the points Pj. However, the matrix 



L(A) = *(A) -ju-^-^A) 



(2.28) 



does not depend on this ordering and is a well defined function on the base curve. Here 
fi is the diagonal matrix /2 = diag (/ii, • • • , /ijv). 
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It should be emphasized however that ^(P) is defined up to normalizations ip^ — > 
dkipk- In fact the normalization ambiguity of the i^k translates into left multiplication 
of the vector ^(P) by a constant diagonal matrix d = diag (1, cfo, ■ ■ ■ , djv)- On the Lax 
matrix L(A) this amounts to a conjugation by a constant diagonal matrix. Hence the 
object we reconstruct is actually the Hamiltonian reduction of the dynamical system by 
this group of diagonal matrices as emphasized at the beginning of this chapter. 



2.6 Separated variables. 

We show here in the example of the Jaynes-Cummings-Gaudin model that the coordi- 
nates of the dynamical poles of the eigenvectors form a set of separated variables. A 
general proof exists but it is cumbersome. 

Let us write the Lax matrix as 

(A{\) B{\)\ 
L{X) -{C(X) -A(X)) 

The spectral curve reads det(L(A) — jj) = or 

r : n 2 - A 2 {\) - B{X)C(X) =0 

Clearly we have 

where Q2n+2(A) is a polynomial of degree 2n + 2. Defining y = fJ-Ylj(X-ej), the equation 
of the curve becomes 

V 2 = <22n+2(A) 

which is an hyperelliptic curve of genus n. The dimension of the phase space of the 
model is 2(n + 1). However, we have to reduce that model by the action of the group of 
conjugation by diagonal matrices which in our case is of dimension 1. Hence we confirm 
that 

g = -dim -M rc duccd = n 

The generator of the group is the Hamiltonian H n as shown by eq.(1.38). We can 
compute easily 

1 , / 4 „„ 2 ^ SHj 



VV2n+2(A) V g 9 a 



3=0 ■ ■ €j 



The coefficients of 5Hj are polynomials of degree n — 1 = g — 1, hence give rise to 
holomorphic differentials. The coefficient of 5H n is a polynomial of degree n = g and 
lead to a singular differential. It is absent in the reduced model where 5H n = 0. 
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At each point of the spectral curve, we can solve the equation 

(L(A) = 

Normalizing the second component (instead of the first for later convenience) of to be 
1, we find 



1 J' T C(A) 

Hence the poles of are located above the zeroes of C(A). Note that if C(Afc) = 0, then 
the points on T above A& have coordinates (ik = ±A(Afc). The pole of ^ is at the point 
Hk = ^4(Afc), since at the other point fj,^ = — A(X/ t ) the numerator of tpi has a zero. 

Recalling that 

C (A)-g + g^— gI3i#^ (2.29) 

This shows that indeed the eigenvector has n = g dynamical poles. 
At infinity, we have two points 



1 

9 

Remembering that 



Q± : n = ±^ (2A - w)(l + 0(A~ 2 )) 



± 



A(A) = ^(2A-w) + 0(A- 1 ) 

C(A) = j + OiX- 1 ) 

we find the following behavior of the function ipi at the two points Q 

Q+ : = i(2A-^) + 0(A- 1 ) 
Q_ : ^i = 0(A _1 ) 

showing that the eigenvector has a pole at Q+ and a zero at Q_ in agreement with the 
general result. 

We can now compute the symplectic form in the coordinates \k,fik- From the con- 
straint {Sj) 2 + s^sj = s 2 we we can eliminate sj. Remembering the Poisson bracket 
{s?, s'f} = is'j ', we can write the symplectic form as 



n = -i5b f\8b + i2_^ -qr- A 6s Z j 
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From eq.(2.29), we see that s+ is the residue of C(A) at A = 6j, so that 



it follows that 



therefore 



But 



Therefore 



n = -iSb A (56 + i j- A 2^ <5s| + 1 2^ 2^ — ^ 

6 i fc j Afc_ej 



6A(X k ) = %5X k + V — ^- - gj a fc 



O = _iJbA«56 + «^- A ^ + i ^ 5A fc ASA(X k ) 



.8b 



5{bb) + Ssf\ + i SX k A M(A fc ) 



Finally 



f2 = z<5 log 6 A SH n + i ^ 5Xk /\5jji k 



This shows that the variables X k ,fi k are canonically conjugate. Remark that the sepa- 
rated variables are invariant under the diagonal group action 



{H n ,X k } = 0, {H n ,fi k } = 
so that they are really coordinates on the reduced phase space. 



(2.30) 



Let us explain why the variables (X k , fi k ) form a set of separated variables. Consider 
the function 



S({F j },{X k })= / a = J2 

J mo k J Aq 



fi k dX k 



The integration contour is done on the level manifold Fj = fj, and jj, is obtained from 
the spectral curve. Just as in the proof of the Liouville theorem, this function does not 
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depend on local variations of the integration path. It is explicitly separated since it can 
be written as a sum of functions each one depending on only one variable X k : 

S({F j },{X k }) = J2Sk({F j },X k ) 

k 

Since Ylk^kdXk is the canonical form, this function is a reduced action and satisfies 
the Hamilton-Jacobi equation. It depends on n arbitrary constants Hj and is thus the 
complete integral of the Hamilton-Jacobi equation. The variables have been explicitly 
separated. 

2.7 Riemann surfaces and integrability. 

Consider a curve in C 2 



T : R(X, fi) = Ro(X, ri + Y, R j( X > ^) H 3 = 

j'=i 



(2.31) 



where the Hi are the only dynamical moduli, so that Ro(X, fj,) and Ri(X, /x) do not contain 
any dynamical variables. If things are set up so that T is of genus g and there are exactly 
g Hamiltonian Hj, then the curve is completely determined by requiring that it passes 
through g points (Aj,/Uj), i = 1, ■ ■ ■ ,g. Indeed, the moduli Hj are determined by solving 
the linear system 



Y j Rj{Xi,Hi)Hj + R (Xi,fii) = 0, i = l,---,g 



whose solution is 
where 



H = -B~ L V 



(2.32) 
(2.33) 



H 



H 



/i?i(Ai,/xi) ••• R g {\i,ni)\ 



B 



Ri(Xi, m 



Rg{Xi, Hi 



V 



\Rl(Xg, Hg) ■■■ Rg(Xg,Hg)/ 



/i?o(Ai,Mi)\ 



Ro{Xi,fa) 



\R0(Xg,fJ,g) / 



Here, of course, we assume that generically det B / 0. 



Theorem 7 Suppose that the variables (Aj,/ij) are separated i.e. they Poisson commute 
for i^j: 

{Xi,Xj} = 0, {m, fj,j} = 0, {Xi,(ij} =p(Xi,Hi)5ij (2.34) 
Then the Hamiltonians Hi, i = 1 ■ ■ ■ g, defined by eq.(2.33) Poisson commute 



{H i ,H :i } = 
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Proof . Let us compute 



BiS 2 {(b- 1 k)i, (£rV) 2 } = { J B 1 ,5 2 }(s-V) 1 (5- 1 y) 2 

-{Bi, F 2 }(£rV)i - {V u B 2 }(B- 1 V) 2 + {V lt V 2 } 
Taking the matrix element i,j of this expression, we get 

(B 1 B 2 {(B- l V) l ,(B- 1 V) 2 }) = SijT^BuKB-WUB^Vh 

-6ij Y,{ B ik, Vi}{B- l V) k - Sij BaUB-'V)! + 5 i3 {V u V{\ = 



where 5^ occurs because the variables are separated. ■ 

It can hardly be simpler. The only thing we use is that the Poisson bracket vanishes 
between different lines of the matrices, and then the antisymmetry. We did not even need 
to specify the Poisson bracket between Aj and fa. The Hamiltonian are in involution 
whatever this Poisson bracket is. This is the root of the multihamiltonian structure of 
integrable systems. 

Lax matrices built with the help of coadjoint orbits of loop groups lead to spectral 
curves of the very special form eq.(2.31) where the Hj are the Poisson commuting Hamil- 
tonians. The coefficients Rj(X,fi) have a simple geometrical meaning. As we have seen 
varying the moduli Hi at A constant one has 



Sfi dX = holomorphic 



(2.35) 



Since 



we see that the coefficients Rj(X, /j,) are in fact the numerators of a basis of holomorphic 
differentials on V: 



Define the angles as the images of the divisor (Xk,Hk) by the Abel map: 



where aj(X,(i)dX is any basis of holomorphic differentials. This maps the dynamical 
divisor {Xk, fJ-k} to a point on the Jacobian of T . 
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Theorem 8 Under the above map, the flows generated by the Hamiltonians Hi are linear 
on the Jacobian. 

Proof . We want to show that the velocities d^Oj are constants, or 



One has (no summation over k) 
dti^k = {Hi, Afc} 



= -{B u l V h \ k } 

= B^i^rs, Afcl-ET^VJ - B'^m, X k } 



= -B 



ik 



{Bks, Xk}H s + {Vk, Xk} 



where in the last line, we used the separated structure of the matrix B and the vector 
V. Explicitly (p(X, ji) = 1 in eq.(2.34) for coadjoint orbits) 

duXk = B~f} d^R s (Xk, Hk)H s + d^RoiXk, Hk) 

or 



duXk — B ik 1 d fJ ,R(X k , Hk) 



which we rewrite as (remember that Bkj = Rj(Xk, fJ>k)) 

r, > Rj(X k ,fl k ) R -1 R 

Recalling eq.(2.36), we have shown 



(2.37) 



(2.38) 



(2.39) 



There are cases where the condition eq.(2.35) is modified. This happens for instance 
when coadjoint orbits are non generic, or as we will see in the next Chapter when the 
Lax matrix belongs to the group G* instead of the Lie algebra Q*. In those cases the 
generalized condition reads 



5fi 



f(X, n) 



dX = holomorphic 



(2.40) 



Obviously the counting argument still works in this case. Moreover by adapting the 
function p(A,/x) entering the Poisson bracket, eq.(2.34), we can preserve the condition 
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that the flows linearize on the Jacobian. The condition eq.(2.40) defines the holomorphic 
differentials as 

* j (X,»)dX = fn R { { ^l , d\ 

Then eq.(2.39) becomes 

T^duXk ^'^ (Afe,/^) = Sij (2.41) 
which produces a linear flow on the Jacobian when f(X,fi) = p(X,fi). 

2.8 Solution of the Jaynes-Cummings-Gaudin model. 

In the Jaynes-Cummings-Gaudin model, the spectral curve reads 

i? = 1(2A- ^) 2 + - 2 H n + - 2 Y^^- + E7V^A2 

so that 

Rj{\n) = M\ M) = V + 4(2A - + + £ 

3 * - ej g g '—r [a — ej) 

Imposing that the points (A&, belong to the spectral curve, we get the set of equations 



\^ n 3 _ y „2 1 m ,,n2 9 rr V ST S 3 ' * J 



\2 



The matrix i?^ is the Cauchy matrix 

= - J— (2.42) 

Since the Hamiltonian is H = ujH n + J2j Hj the equation of motion takes the form, 
using eq.(2.37) 

a* = E ^ A * = -^ 2 ^ E V ( 2 - 43 ) 

or 

.• . 2 rij(Afe-ej) . 2 \/ ( 52n+2(Afc) 

Afc = ig Mfc rr = «ff 



lli^fc(Afe - Aj) llz^(Afc - A ; ) 
(2.44) 



Let us check that eq.(2.44) is of the form eq.(2.39). In fact, multiplying eq.(2.43) by 
A l jfc 1 -Bfci with Bfa given by eq.(2.42), and summing over k, we find 

EBki . 2 
Afc = 
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or equivalently 

The coefficients of A& are precisely a basis of holomorphic differentials. In general such 
an equation is solved by the Abel transformation and ^-functions (see e.g. [5]). 

2.8.1 Degenerate case. 

Let us consider the degenerate case where 

4£L _. 2 



Q 2 n + 2(X) = ^H(X-E l 



y 1=1 

This happens for instance when we start from an initial condition where 

6 = 6 = 0, sf = 0, s] = ejs, ej = ±1 (2.45) 
The energy of this configuration is 



H = 2s e 



j 



At time U we have B(\)\ t=ti = C{\)\t=t l = and 



The zeroes of Q2ra+2(A) are located at the zeroes of ^(A)^^ and are thus the roots of 
the equation 

3=0 3 



This equation has also a remarkable interpretation. Its solutions E\ are the eigen 
frequencies of the small fluctuations around the configuration eq.(2.45). To perform the 
analysis of the small fluctuations around this configuration, we assume that 6, 6, are 
first order and s z - = sej + 5s j. Then 5s z - is determined by saying that the spin is of 
length s and is of second order. 

5s z = -^-s~s + 

This is compatible with eq.(1.31). The linearized equations of motion are 

6 = -icob-ig^sj (2.47) 
3 

sj = -2i€jSj + 2isgejb (2.48) 
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and their complex conjugate. We look for eigenmodes of the form 



b(t) = 6(0)e 



-2iEt 



s; = sj(0)e 



-2iEt 



We get from eq.(2.48) 




Inserting into eq.(2.47), we obtain the self-consistency equation for E 



lo sg 2 sr^ £j 



(2.49) 



2 2 ^ E - ej 

A J 



3 



which is exactly eq.(2.46). 

Let us assume first that ej > 0. The minimal energy state among the configurations 
eq.(2.45) corresponds to all the spins down: ej = —1, j = 0, ■ ■ ■ ,n — 1. The graph of 
the curves y = A — to/ 2 and y = Yl]=o when all spins are down is presented in 
Fig. [2.1]. We see that we have n + 1 real roots, meaning that the state is localy stable. 



Let us assume next that €j < 0. The minimal energy state among the configurations 
eq.(2.45) now corresponds to all the spins up: ej = 1, j = 0, ■ ■ ■ , n — 1. The graphs now 
look like Fig[2.2]. We see that we have n + 1 real roots if lo > u) sup or lo < uJi n f. In 
between we have n — 1 real roots and a pair of complex conjugate roots, which means 
that an instability develops. The question then arises to determine the time evolution 
of the non linear system. 

The equations of motion of the separated variables are in that case 




5 



Figure 2.1: The solutions of eq.(2.49) when ej > 0, ej = — 1. 



Xi = 2i 



(2.50) 
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Figure 2.2: The solutions of eq.(2.49) when ej < 0, ej = 1. 



One can solve explicitly these equations, hence finding the solution of the reduced 
model on the unstable surface. From eq.(2.50), we have 



A, 



2i- 



Uk^li^i ~ E k) 



Therefore 



EA, 



A,; — Ei 



2i 



Aj - Ei Ylj^ii^i ~ Aj 



2i(Ei- <7i + 



where is a big circle at infinity surrounding all the Xj. Hence 

rt 



log ]J(Ei - Xj) = 2iE l t + 2iJ di(Si(t) - <7i(S)) 



where Ei(t) = ^ • Aj(i). Define 



so that 



Introducing 



log 



7 = 2i(Ei(t)-ffi(£))7 

7 - 1 (i)J](^-A J ) = 2iE,t + C, 
i 

P(A,t) = n(A-A,) 

j 
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the above equation becomes (we parametrize e Cl = A\ Ylk^ii^i ~ ^k) f° r convenience) 
V{E U t) = -r(t)Ai l[(Ei - E k ) e 2iEl \ I = 1, • ■ ■ n + 1. 

These are n + 1 conditions on the polynomial V(X, t) of degree n in A. It can then be 
reconstructed by the Lagrange interpolation formula 



P(X,t) = 7 (t)^V iEit I](^ E k ) 

i kf-l 



(2.51) 



Imposing that the term of degree n should be normalized to A n , we get 



lit) 



E t Aie 2lElt 



(2.52) 



From 



V(X,t) = 1 {t)i^A l e 2iE ^\ n - 1 {t)i^{a{E)-E l )A^ 



we deduce that 



EiAiE ie ^ _ 17 



which shows the consistency of the construction. The original variables are given by 
b{t) 



b(t) = 7 (*) 

2 -f(t)r(ej,t) 



9llk^j( e j- e k) 

4 V(€j,t) 



s + 



9 2 Uk^j( e j- e k) 



€j + Si(t)- <7i(e)- 



This is an exact solution of the Jaynes-Cummings-Gaudin model. It can be viewed as a 
non linear superposition of the eigenmodes of the system. 

2.8.2 Non degenerate case 

In the non degenerate case, we can find the generalization of eq.(2.51). The result is 
eq.(2.66). Consider the hyperelliptic curve 



2g+2 
i=0 



V 2 = P2g+2(x) 



(2.53) 
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We choose a partition of the branch points into g + 1 disjoint pairs. The branch points at 
which the cuts start are called ai and the end points of the cuts are called bi, i = 0, 1, ■ ■ ■ g. 
A basis of holomorphic differentials is 



x % l dx 


i = l,-- 






•,9 


y 





The dual basis of second kind differentials is 



o + 1 - j)pk+i+j 



x dx x 3 d 



(P 2 g+ 2 ( x ) x 2j ) + , j = l,---,9 



k=j 



4y Ay dx 



where ()+ means the polynomial part in the expansion at x = oo. The first thing to 
check is that rjj is a second kind differential, i.e. has no residue at infinity. For this, we 
remark that 



TyTx {P2 ^ x)x ~ 23)+ = TyJx\ 



(P 2g+2 {x)x- 2 i - (P 2g+2 (x)x- 2 i)^ 



(2.54) 



The first term is a total derivative and cannot have a residue. The second term is regular 
at infinity. To compute (rjj • u>i), let us apply eq.(2.12). Using eq.(2.54) we see that we 
can take g\ = \ {yx~i) + regular. Hence 

(rjj • uJi) = ^res 00 (x l ~ j ~ 1 dx + regular) = ^5ij 
The periods of these differentials are denoted as follows (i, j = 1, • • • , g) 

UJi 



a 3 

( 2 v)ij = / m, ( 2 v')ij = Vi 

J cij J hi 



The Riemann bilinear identities imply 



J 1 uj - uj t J = 0, r]' t u - r) t J 
Let us introduce the symmetric function 



ITT 



Id, ri^rj-r] V = 



F(XI,X 2 ) = P 2 g+ 2 {^/XlX 2 ) + P 2g+2 (-^/XlX 2 ) + {P 2 g+ 2 {^/XlX 2 ) - P 2g+2 {-^XlX 2 )) 

l^JX\X 2 

9+1 9 



2 S ^ i V2iX\x\ + (Xl + X 2 ) y^p 2 i + ix\x l 2 



i=0 



i=0 
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With it, we construct the fundamental symmetrical bi-differential 



£l(xi,x 2 ) = 



2yiy 2 + F(xi,x 2 ) dx\ dx 2 
A(x\ - x 2 ) 2 yi y 2 



If x 2 — > X\ + e, we have 

F(x u x 2 ) = 2P 2g+2 { Xl ) + eP^ g+2 ( Xl ) + 0(e 2 ) = 2y{x 1 )y{x 2 ) + 0{e 2 ) 
from which it follows that 

tt(xi,x 2 ) = ( — + O(l) J dziete 2 

X 2 ) J 



Alternatively, we can write 



£l(xi,x 2 ) 



Of y\+V2 



dx 2 \2yi(x\ — x 2 
9 ( y\+V2 



dx x dx 2 + ^2uji(xi)r)i(x 2 ) 
i=i 

)dxidx 2 + V] 7]i(xi)LOi(x 2 ) 

oxi \2y 2 {xi - x 2 ) J ^ 



(2.55) 

(2.56) 
(2.57) 



i=i 



In general the fundamental Kleinian cr-function is defined by 



a(z) 



1 \ 4 / tt9 i 



D(F) J V det(2u/ 



ezWw l \*)e[e R ](z,u,tJ) 



where [cr] is the characteristic of the vector of Riemann constants. The function D(T) 
is the discriminant of the equation R(\, /x) defining the curve T. Its main property is 
that it is invariant under modular transformations 

a(z,LO,co') = cr(z,u,(jj') 

where (a),u/) is related to (o>,u/) by a Sp(2g,Z) transformation. From this we define 
the Kleinian £ and p-functions by 

dloga(z) . 
Ci(«) = ~ . *=1," -5 



and 



a 2 iog ( j(z) 

dzidzj 



The functions pij(z) are automorphic functions with respect to the Sp(2g,Z) transfor- 
mations. 
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The functions cr(z), Q( z ) an d Pij( z ) have the following periodicity properties: 

<t(z + \Q(m,m'))) = e (. E ( m > m ')\ z + n ( m > m ')) e -in(m\m)+2in((m\q')-(m'\q)) (T ^ 

(i(z + n(m, m')) = Q(z) + E^m, m!) 
pij(z + Q(m, m')) = pij(z) 
where we have introduced the vectors of periods 

\E(m, m')) = 2r]\m) + 2rj'\m'), \Q(m, m!)) = 2uj\m) + 2oj'\m!) 

For an hyperelliptic curve, things can be made more explicit. The function cr(z) 
reads 

a ^ = e ( 2 |(2^)- 1 r,| Z >+4i7r( g n(2^)- 1 | Z >+j7r( (/ '|r|'?'>-2ivr( g '|'7> 6 )((2 w )-l z _ # ) 



where K ao is the vector of Riemann constants 



k=l Ja ° 



Because ao is a branch point, K ao is a half period and can be written as K a = q + r</ 
with half integers g and g'. The characteristic [cr] is 



v 



Notice that if 



3 

fc=l l/a ° 



we have (2a;) 1 z — K ao = (2u) 1 u where 



9 




Ui = ^_ 


^ j Xk x % l dx 


k= 


-l Ja k V 



(2.58) 



Theorem 9 Let (ao,y(ao)), (x,y) and be arbitrary points on T. Let 

{(xi,yi), (x g ,y g )}, {(fiuvi),-- ■ , (Vg^g)} 
be arbitrary distinct points. Then the following relation is valid 



rx 9 rx 

/ E / 



Xi 



log < 



K/>-5Xi -IXiJ>) 
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Proof . Using eq.(2.55), the left hand side in this formula behaves like log(x — Xi) when 
x — ► xi and — log (a; — m) when x — > m and these are the only singularities. Moreover 
it vanishes when x — > /j,. This behaviour is clearly reproduced by the right hand side. 
Then using eq.(2.56) one can check that the periods on both sides are the same. Hence 
the two expressions are identical. ■ 



Theorem 10 (Bolza) 



Ja k=Q Ja k * h=n 



9 \ z—Xk 
k=0 



z=x k 



R'(x k ) 



Proof . Setting /ij = a«, we get 

9 



I (7 I I '' UJ // ) <T f f UJ 

«(*,**) = log ) ao ; Kao 

We now take the derivative of both sides with respect to Uj. We obtain 

rx 9 



rx a f j' X \ ( \ 

I e«<*. **> 4=-<> {J m "-J+ <> u. - - "J 



With the help of eq.(2.57) we find 

(j ( I u - u J - Q ( ! uj- 

\J an / \J an 



1 ^ 1 dx k 

2 ^ y k duj L 



fc=l 



i k 



V + Vk _ v + Vk 
diij 



From eq.(2.58) we have 



Hence 

0(/ w - u j - Ci ( / w - « 

V ■/ ao / V ■/ ao 



y dx rjj(x) + - ^ 



(2.59) 



(2.60) 



1 dx h 
Vk duj 



V + Vk _ v + Vk 



Applying the hyperelliptic involution a(x, y) = (x, —y), cr(fi, v) = (/i, —v), we get as well 



Cj[ I u + u j - Cj ( [ u + u 

ao J \J an 



f X 1 9 

/ dxr} j (x) + -^ 
J n 1 k=i 



y -yk v -yk 



1 dxk 

y k duj [x - x k n - x k 
(2.61 
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Let us introduce the Vandermonde matrix 



i-i 



then eq.(2.60) gives 



Vki = 4- 



1 dx k _ ! = (P(z)z J )+\ z=Xk 



where 



Vkduj * P'{x k ) 
a 

p(z)=n(z- Xk ) 

The last formula is easy to prove 



(2.62) 



k=i 



z )z ^)^.\ z=Xk xj 



but 



D^-W" 1 = E-r 1 E^ = E^ E -r 1 ^' = P( 1 _ f 



1 <?'<»■ 



If z — ► Xfe 7^ a;/, we get zero because P(x k ) = P(xi) = 0, while if 2 — > xz, we get P'(xi). 
Hence we have shown that Y2j VijVjk 1 = <W 



Let 

R{z) = {z- x)P(z) 
By simple polar decomposition, we have 

(P(z)z-*) + \ g=x _ = ^ (P(z)z-*) + \ z=Xk 1 



1 



R'{x) 



Next we have the identity 



= Yvr k 1 — 
P'(x k ) x - x k 3k x - x k 



(2.63) 



z—x k 



z=x k 



z-x k 



z=x k 



(P(z)z-i) + \ z=Xk 1 



R'(x k ) 
This is because 
i?(z) 



P'(x fc ) x - Xfc 



(2.64) 



y 9~3-r(_ 



Z — X\ 

(P(z)z-i), = 



(-l) r [xa r _i(x 2 , ■ • ■ ,x g ) + 0>(x 2 , ■■■ , Xg)\ 



£ 

r<g-j 

(-l) r [XlCT r „l(X2, ■ • • , Xg) + (T r (X2, • • • , Xg)\ 

r<g-j 
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so that 



Z — X\ 



It follows that 



A -{P(z)z-*) + = (x-xi) ^- J '- r (-l)V r _i( i r 2 ,---,x ff ) 



^— x k 



z=x k 



Z — Xk 



z=x k 



R'(x k ) 



P'(x k ) 



-1 

jk 



X X k 



(2.65) 



We can now evaluate the sums appearing in the right hand side of eq.(2.61). We 
have 



A 1 dx k y - y k Vy-i 1 V v' 1 - 

^ y k duj x-x k V i k ^ yfe J' fe r - 



fe=i 



k x Xk k 



X Xk 

( ( fl« „- 



y- 



(P(z)z^) 



+ \z=x 



R!{x) 



+E 



Z-Zfc 



+ _ (S— ) 



z=xjt 



P'(x k ) 



\ 



Noticing that P(z) = R(z)/(z — x) and defining (xq, yo) = (x, y), we get 



E 



1 dx k y- y k 
y k duj x-x k 



(m. z -i\ 

9 \z-x k * ) 
k=0 



z=x k 



R'(x k ) 



E 



+ 



Vk~ 



z=x k 



P'(x k ) 



k=l J k=0 v ' k=l 

Defining (x g+ i,y g+ i) = (fi, v), R(z) = (z - n)P(z), a g+1 = a Q , we obtain 



E 

k=i 



1 dx k 



V-Vk v -y k 



y k duj [x - x k fi- x k _ 



2®. r-3 



9 \ z—x k 
k=0 



R>{x k ) 



3+1 \z-x k 
z=x k \ 

— --l^Vk 

k=l 



z=x k 



R'{x k ) 



Hence 



z—x k 



/ rx \ 9 i-x k \ ® ( 

\L u+u ) + £ L mix) ~ t £* — 



z=:r fe 



9+1 



1^ ( 



/ ft 1 \ /"^fc 1 

o / +E / #)-oE 

V-'ao / k=1 Ja k z fc=1 



2=X fe 



i2'(x fc ) 
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The left hand side is a symmetric function of x, x\, ■ ■ ■ , x g , but the right hand side 
does not depend on x. Therefore the whole expression is a constant independent of 
x, xi, • • ■ ,x g . Applying the hyperelliptic involution, we see that this constant vanishes. 



The main interest in the Kleinian functions is that they give a very explicit solution 
to the Jacobi problem. 

Theorem 11 Let T as in eq.(2.53). Then the preimage of the point u € Jac(r) is given 
by the set of points (P±, P2, • • • , P g ) where (x±, X2, • • • , x g ) are the zeroes of the polynomial 



V(x, u) = x 9 - V g {u)x 9 - 1 - V g -i{u)x 9 - 2 Vi(u) 



(2.66) 



where 




The coordinates (yi, y2, • • • , y g ) are given by 

dP(x,u) 



Vk = - 



The constants Ci are given by eq.(2.69). 



du„ 



(2.67) 



Proof . We now take the limit (x,y) — > P± in eq.(2.59). When x — > 00 we have, by 
eq.(2.65) 



+ 



\z-x k 

lim — 



z-x k 



z=x k 



z=x k 



P'(x k ) 



On the other hand 



m z -j 

z-x * 



2/0*0 - 



R'(x) 



= y(x) 



= y{x) 



P(x) 



P(x)x' j - (P(x)z- j )^ 
Pjx) 



When x — > 00, the first term behaves as y(x)x 3 ~ (y(x)x ■ ? ) + , while in the second 
term, (P{x)z~i) _ ~ (— l) s_ - 7+1 cr g _j + i(xi, X2, ■ ■ ■ , x g )x~ 1 + 0(x~ 2 ) and we find 



y(x) 



z-x 



R'(x) 



~ {y(x)x J ) + - ^/ P 2g + 2(-l) 9 J+1 cr g - j+1 (x 1 ,X2,' 



, Xg) 
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Since the limits of the left hand side of eq.(2.59) are finite, so are the limits 

cf>= lim (- T m + h yx ~r )+ \ (2 . 68) 

(x,j/)->P± \ J ao 2 J 

and they are obviously independent of x±, ■ ■ ■ , x g . Therefore we arrive at 

(j ( UJ + U j-Cj[ Uj + Uj =Cj + ^g^i-lY'^g-j+l (x!,x 2 , ■ ■ ■ , x g ) 

\J a () / V J ao / 



'ao / \J ao 

where we have set 



Cj = cf ] - cj. } (2.69) 



Multiplying by x- 7 1 and summing over j = 1, ■ ■ ■ , g, we get 

fP+ \ / rP- 



x3 1 (o (/ w + u ) - (/ w + u ) - c i) = - v^+2 \ Jl( x ~ Xfc ) 



Hence we have shown 



k l j=l 



which is eq.(2.66). From eq.(2.62), we have 

dx k y k 



du 9 Ul^k( x k ~ xi) 

and 



so that 

9 - X/) 



du 9 

This proves eq.(2.67). 



= -Vk 

x=x k 



For applications to the Jaynes-Cummings-Gaudin model, u = u(t) is chosen to be a 
linear function of time. Once the (xk(t),yk(t)) are found by solving the inversion Jacobi 
problem, it is easy to reconstruct the original variables of the model. 
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Chapter 3 

Infinite dimensional systems. 



3.1 Integrable field theories and monodromy matrix. 



For a system with a finite number of degrees of freedom, we have seen that a Lax matrix 
could be interpreted as a coadjoint orbit. In field theory we have seen that one possibility 
is to consider the algebra of pseudo differential operators. We obtained in this way the 
KdV hierarchies. A more general approach consists is starting from the zero curvature 
representation which assumes that the field equations can be recast in the form 



d t u - d x v - [v, u] = o 



(3.1) 



As before the matrices U and V will in general depend on an extra parameter A. The 
zero curvature (3.1) condition expresses the compatibility condition of the associated 
linear system 



(d x -U)* = 0, (dt-V)* = 



(3.2) 



The matrices U and V can be thought of as the x and t components of a connection. 
This connection will be called the Lax connection. Given U and V, the linear system 
(3.2) determines the matrix Vl/ up to multiplication on the right by a constant matrix, 
which we can fix by requiring Vl/(A, 0, 0) = 1. This ^ will be called the wave function. 

Choosing a path 7 from the origin to the point (x, t) the wave function can be written 
symbolically as 




^(x,t) =exp 

J -V 

(3.3) 

where exp denotes the path-ordered exponential. This is just the parallel transport 
along the curve 7 with the connection (U,V). Since the Lax connection satisfies the zero 
curvature relation (3.1) the value of the path-ordered exponential is independent of the 
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choice of this path. In particular if 7 is the path x G [0, 2ir], with fixed time t we call 
^(27r,t) the monodromy matrix T(X,t): 



T(X, t) =exp 



2tt 



U (A, x, i)cfcc 



L/0 



(3.4) 

where we assume that U(X,x,t) and V(X,x,t) depend on a spectral parameter A. 

It is the monodromy matrix which plays the role of the Lax matrix in the field 
theoretical context as the following proposition shows 

Proposition 12 Assume that all fields are periodic in x with period 2tt. Let T(A, t) be 
the monodromy matrix. Its time evolution is given by the Lax equation 



d t T(X,t) = [V(X,0,t),T(X,t)\ 



(3.5) 



As a consequence the quantities 



fl»(A) = Tr{T n {X,t)) 

(3.6) 

are independent of time. Hence traces of powers of the monodromy matrix generate 
conserved quantities. 

Proof . Thinking of the path-ordered exponential on [0, 2n] as 

f27T 



exp 



U (x)dx 



{l + SxU(x n ))---(l + SxU( Xl )) 



with a subdivision x\ = < X2 < • • • < x n = 2ir such that Xi + \ — Xi = Sx — > 0, we get 
(all exponentials are path-ordered exponentials): 

r 2w 



d t T(t) 



r 2 

= / dxe^ Udx dtU(x) efo Udx 
Jo 

= / dxef*" Udx (d x V + [V, U])e% Udx 
Jo 

= dxd x [e^ Udx Vefo Udx ^j 

Performing the integral, 

d t T(X, t) = V(X, 2vr, i)T(X, t) - T(A, t)V(X, 0, t) 



(3.7) 



So, if the fields are periodic, we have V(X, 27r,i) = V(X,0,t) and we obtain eq.(3.5). 
This is a Lax equation. It implies that H^ n \X) is time independent. Expanding in A we 
obtain an infinite set of conserved quantities. 
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3.2 Abelianization. 

We consider the linear system eq.(3.2) where U(X,x,t) and V(X,x,t) are matrices de- 
pending in a rational way on a parameter A having poles at constant values X k . 

-i 

U = U + J2Uk with U k = U k , r {X-\kY (3.8) 

k r=-n k 

-1 

V = Vo + Y j V k with V k = £ ^(A-A fe ) r (3.9) 



r=-m k 



The compatibility condition of the linear system (3.2) is the zero curvature condition 
(3.1). We demand that it holds identically in A. 

As for finite dimensional systems we first make a local analysis around each pole 
Afc in order to understand solutions of eq.(3.1). Around each singularity X k , one can 
perform a gauge transformation bringing simultaneously U(X) and V(X) to a diagonal 
form. The important new feature that must be to emphasized, as compared to the finite 
dimensional case, is that this construction is local in x. We have 

Proposition 13 There exists a local, periodic, gauge transformation 

d x -U = g^(d x -A^)g^-\ d t -V = g^ k \d t -B^)g^- 1 (3.10) 

where g^ k \X), A^ k \X) and B^ k \X) are formal series in X — X k 

oo oo oo 

r=0 r=—n r=—m 

such that the matrices A^ k \X) and B( k \X) are diagonal. Moreover dtA ( - k \X)—d x B ( - k \X) = 
0. 

As for finite dimensional systems, we can reconstruct all the matrices U k and V k , and 
therefore the Lax connection, from simple data. 

U = U + Y,U k , with U k =(g^A^g^- 1 )_ (3.11) 
k 

V = Vo + J^Vk, with V k =(gWB«! ) gW- 1 )_ (3.12) 

In this diagonal gauge it is easy to compute the conserved quantities: 

Proposition 14 The quantities Q( fc ) (A) = A^ k \X,x,t) dx are local conserved quan- 
tities of the field theory. They are related to eq.(3.6) by 

(*)/ 



fl»(A) = Tr exp n / A^ k \X,x,t)dx 

L Jo 



Tr exp 



nQ w {X) 
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Notice that there is no problem of ordering in the exponential since the matrices are 
diagonal. 



We now give some examples of 2 dimensional field theories having a zero curvature 
representation. 

Example 1. The first example is the non-linear a model. For simplicity, we look 
for a Lax connection in which U and V have only one simple pole at two different points 
and Uo = Vo = 0. Choosing these points to be at A = ±1, we can thus parametrize U 
and V as: 

"-5^1*. v = -xh J < < 3 - 13 > 

with J x and J t taking values in some Lie algebra. Decomposing the zero curvature 
condition [d x — U, dt — V] = over its simple poles gives two equations: 

dtJx ~ ^[Jx,Jt] = 0, 

d x J t + \[J x ,Jt] = 0. 

Taking the difference implies that [dt + Jt, d x + J x ] = 0. Thus J is a pure gauge and there 
exists g such J t = g~ 1 dtg and J x = g~ 1 d x g. Taking now the sum of the two equations 
implies dt J x + d x Jt = 0, or equivalently, 

dt(g- 1 d x g) + d x (g- 1 d t g) = 

This is the field equation of the so-called non-linear sigma model, with x, t as light-cone 
coordinates. 

Example 2. Another important example is the sinh-Gordon model. It also has a 
two poles Lax connection, one pole at A = 0, the other at A = oo. Moreover, we require 
that in the light cone coordinates, x± = x ± t, U(X, x±) has a simple pole at A = and 
V(A, x±) a simple pole at A = oo. The most general 2x2 system of this form is: 

(d x+ -U)v = o, u = u + \- 1 u 1 

(d x _-V)V = 0, V = V + XV 1 

The matrices Ui, Vi are taken to be traceless matrices, so contain 12 parameters. One can 
reduce this number by imposing a symmetry condition under a discrete group, Namely, 
we consider the group Z2 acting by: 



vfr(A) _^ ^(-AK- 1 , a z =(^ (3.14) 

and we demand that be invariant by this action. This restriction means that the wave 
function belongs to the twisted loop group. It follows that: a z U(—X)a z = U(X) and 
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a z V{— \)(T Z = V(X). We still have the possibility to perform a gauge transformation by 
an element g, independent of A, in order to preserve the pole structure of the connection, 
and commuting with the action of Z2, i.e. g diagonal. This gauge freedom can be used 
to set (Vo)u = 0. The symmetry condition then gives: 

tj = ( u o A _1 ui \ = / Xvi 

VA- 1 ^ -u J ' \Xv 2 

In this gauge, the zero curvature equation reduces to : 

d X -Uo - U1V2 + v\U2 = (3.15) 
d X -Ui = 0, d x _U2 = (3.16) 
d x+ v\ - 2viu = 0, d x+ v 2 + 2v 2 u = (3.17) 

From eq.(3.16) we have u\ = a(x+), U2 = f3(x+). We set uq = d x+ (p. Then, from 
eq(3.17) we have v\ = j(x-) exp2ip and V2 = S(x-) exp — 2</?. Finally eq(3.15) becomes : 

d x+ d x _<p + /3(x+) 7 (x_) e 2 ^ - a(x + )5(x-)e- 2< ? = 

This is the sinh-Gordon equation. The arbitrary functions a(x + ), @(x + ) and j(x-), S(x-) 
are irrelevant: they can be absorbed into a redefinition of the field ip and a change of 
the coordinates x + ,x_. Taking them as constants, equal to m, we finally get 

Hence the Lax connection of the sinh-Gordon model is naturally recovered from two- 
poles systems with Z2 symmetry. This construction generalizes to other Lie algebras, 
the reduction group being generated by the Coxeter automorphism and yields the Toda 
field theories. 

Example 3. The KdV equation reads: 

4d t u = -Qud x u + dlu (3.19) 
The KdV equation can be written as the zero curvature condition 

F xt = d x v - d t u - [u,v] = 

with the connection U, V, depending on a spectral parameter A: 

U-( ° M V-U dxU 4A-2u\ 

U ~\\ + u 0j' V ~ 4 \4X 2 + 2Xu + d 2 x u-2u 2 -d x u J 

Alternatively one can recast the KdV equation in the Lax form dfL = [M, L], where 
L and M are the following differential operators: 

L = d 2 -u (3.21) 

M = ^(4d 3 -6ud-3{d x u)) = (L§)+ 
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The operator d acts as d x , and the notation (L 2 )+ refers to the pseudo differential 
operator formalism. 

Of course these two descriptions are not independent. To relate them, consider the 
linear system: 

(dx-U)(Z)=°> (dt-V)(*)=0 (3.22) 



X J \x. 
The x-equation yields x = d x ^ and 

(L-A)* = with L = d 2 x -u (3.23) 

The time evolution of ^ is given by 4dt^ = d x u ■ ^ + (4A — 2u) d x ^f. Using eq(3.23), 
this may be rewritten as: 

(d t - M)* = with M = ^(4d 3 - 3ud - 3du) (3.24) 

The compatibility condition of eqs. (3.23,3.24) is the Lax equation dtL = [M,L], which 
is equivalent to the KdV equation. 

Eq.(3.23) is the Schroedinger equation with potential u. The parameter A gets an 
interpretation as a point of the spectrum of this operator. This is the origin of the 
terminology "spectral parameter" . 



3.3 Poisson brackets of the monodromy matrix. 

As we just saw, the zero curvature equation leads to the construction of infinite set of 
local conserved quantities. We want to compute their Poisson brackets. For this we will 
compute the Poisson brackets of the matrix elements of the monodromy matrix. 
In order to do it we assume the existence of a r-matrix relation such that: 

{£MA, x), U 2 (fi, y)} = [ri 2 (A - »), U^X, x) + U 2 (n, y)]S(x - y) (3.25) 

We assume that ri2(A — fx) is a r-matrix as in eq.(1.23). We say that the Poisson bracket 
eq.(3.25) is ultralocal due to the presence of 5(x—y) only. This hypothesis actually covers 
a large class of interesting integrable field theories, but certainly not all of them. 

Since we are computing Poisson brackets, let us fix the time t, and consider the 
transport matrix from x to y 

T(X;y,x) =exp (^j U(X, z)dz^J 

In particular the monodromy matrix is T(A) = T(X;2ir,0). The matrix elements [T]ij 
of T(A; y, x) are functions on phase space. We use the usual tensor notation to arrange 
the table of their Poisson brackets. 
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Proposition 15 // eq.(3.25) holds, we have the fundamental Sklyanin relation for the 
transport matrix: 

{Ti(X;y,x),T 2 (fi;y,x)} = [n 2 (A, //), Ti(A; y, x)T 2 (n; y, x)] (3.26) 

As a consequence, the traces of powers of the monodromy matrix H^ n \X) = Tr (T n (A)), 
generate Poisson commuting quantities: 

{tf (n) (A),# (m) (/i)} = (3.27) 

Proof . Let us first prove the relation (3.26) for the Poisson brackets of the transport 
matrices. Notice that A is attached to T\ and /x to T 2 , so that there is no ambiguity if we 
do not write explicitly the A and fx dependence. The transport matrix T(y, x) verifies 
the differential equations 

d x T(y,x) + T(y,x)U(x) = (3.28) 
d y T(y,x)-U(y)T(y,x) = 

Since Poisson brackets satisfy the Leibnitz rules, we have 

{T 1 (y,x),T 2 (y,x)} = (3.29) 

duduTi(y,«)T 2 (y,t;){i7i(u),?72(i;)}ri(«,ar)r2(i;,x) 



Replacing {U\{u), U 2 (v)} by eq.(3.25), and using the differential equation satisfied by 
T(y, x) this yields: 

{Tx(y,x),T 2 (y,x)} = 

[ [ dudv5(u-v).(T 1 (y,u)T 2 (y,v)r 12 (d u + d v )T 1 (u,x)T 2 (v,x) 

+{d u + d v ){T 1 {y,u)T 2 {y,v)) r\ 2 Ti(u,x)T 2 (v,x)^ 

dzd z (Ti(y, z)T 2 (y, z).ri 2 .Ti(z, x)T 2 (z, xfj 



Integrating this exact derivative gives the relation (3.26). Let us now show that the trace 
of the monodromy matrix H^ n \\) generates Poisson commuting quantities. Eq.(3.26) 
implies 

{Tf(A),T 2 ™(M)} = [r 12 (A, M ),T 1 "(A)T 2 m (^)] 

We take the trace of this relation. In the left hand side we use the fact that TV12 (A®B) = 
Tr (A)Tr (B) and get {# (n) (A), if( m )( / u)}. The right hand side gives zero because it is 
the trace of a commutator. 

Let us emphasis that it is the integration process involved in the transport matrix 
which leads from the linear Poisson bracket eq.(3.25) to the quadratic Sklyanin Poisson 
bracket eq.(3.26). 

The proposition shows that we may take as Hamiltonian any element of the family 
generated by H^ n \ii). We show that the corresponding equations of motion take the 
form of a zero curvature condition. 
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Proposition 16 Taking H^ n \fi) as Hamiltonian, we have 

£7(A,x) = {HW(n),U(\,x)} = d x V<ri(\,n,x) + [V^(X,fi,x),U(X,x)} (3.30) 
where 

V( n \\, fi; x) = nTn {t, (/x; 2tt, x)r 12 {fi, A)7i (/*; x, O)^ 1 ( M , 2vr, 0)) 

T/iis provides the equations of motion for a hierarchy of times, when we expand in [i. 

Proof . To simplify the notation, we do not explicitly write the A, [i dependence as 
above, noting that fi is attached to the tensorial index 1 and A to the tensorial index 2. 
We have: 

/•27T 

{T 1 (27r,0),f/ 2 (x)} = / dy ^(271^) {U 1 (y),U 2 (x)}T 1 (y,0) 

Jo 

= Ti(27r,x) [ri 2 ,C/ 1 (x) + C/ 2 (x)] Tx^O) 

Expanding the commutator we get four terms 

{T 1 (2n,0),U 2 (x)} = 

ri(27r,x) • na • C/i(x)T!(x,0) +Ti(2ir,s) • r 12 • U 2 {x) Ti(s,0) 

N v ' V v ' 

use fifi//. eq. commute 

- ri(27r, x) C/i(x) -na • Ti(s, 0) - Ti(27r, x) C/ 2 (x) -r 12 • T^x, 0) 

N v ' " v ' 

use diff. eq. commute 

Using the differential equations (3.28) and commuting factors as indicated gives 

{r 1 (27r,o),c/ 2 (x)} = d x v 12 (x) + [v 12 (x),u 2 (x)} 

where we have introduced V± 2 (x) = Ti(2ir, x)-ri 2 -T±(x, 0). From this we get {T±(2ir, 0), ?7 2 (x)} = 
d x V${x) + [V^\x),U 2 {x)\ with V${x) = Y,i T! r i ~ 1 yi2{x)T{. Taking the trace over 
the first space, remembering that H^ n \fj,) = TrT n (fi), and setting V^ n \\, fi, x) = 
TriV^x), we find eq.(3.30). 

3.4 Dressing transformations. 

We now introduce a very important notion, the group of dressing transformations, which 
is related to the Zakharov-Shabat construction. These transformations provide a way 
to construct new solutions of the field equations of motion from old ones. It defines a 
group action on the space of classical solutions of the model, and therefore on the phase 
space of the model. 

Dressing transformations are special non-local gauge transformations preserving the 
analytical structure of the Lax connection. These transformations are intimately related 
to the Riemann-Hilbert problem which we have discussed in the section on factorization. 
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We choose a contour T in the A-plane such that none of the poles A& of the Lax 
connection are on T. We will take for T the sum of contours r( fe ) each one surrounding 
a pole Afc as in the factorization problem. 

To define the dressing transformation, we pick a group valued function <?(A) £ G on 
r. From the Riemann-Hilbert problem, g{X) can be factorized as: 

g{X)=g-_ 1 {X)g + {\) (3.31) 

where g + (A) and <7_(A) are analytic inside and ouside the contour V respectively. In the 
following discussion we assume that g{\) is close enough to the identity so that there 
are no indices. 

Let U, V be a solution of the zero curvature equation eq.(3.1) with the prescribed 
singularities specified in eqs.(3.8,3.9). Let ^ = ^(X;x,t) be the solution of the linear 
system (3.2) normalized by ^(A;0,0) = 1. We set: 

0(A; x, t) = *(A; x, t) ■ g(X) ■ *(A; x, t)~ l (3.32) 

At each space-time point (x,t), we perform a A decomposition of 6(X,x,t) according to 
the Riemann-Hilbert problem as: 

0(A; x, t) = 9Z 1 (X; x, t) • #+(A; x, t) (3.33) 

with 6+ and #_ analytic inside and outside the contour V respectively. Then, 

Proposition 17 The following function, defined for X on the contour V , 

* 9 (A; x, t) = 9±(X; x, t) ■ *(A; x, t) ■ g±\X) (3.34) 

extends to a function ^+ defined inside T except at the points Xk where it has essential 
singularities and a function ^> 9 _ defined outside T. On T we have ^ 9 _~ 1 ^ 9 + \t = 1- So 
Vl/rj. define a unique function ^ 9 which is normalized by ^ 9 (A,0) = 1 and is solution of 
the linear system (3.2) with Lax connection U 9 and V 9 given by 

U 9 (X;x,t) = 0± • U • eg 1 + d x 6± ■ (3.35) 
V 9 {X;x,t) = 6± ■ V ■ 6±* + d t 9± ■ 6±* (3.36) 

The matrices U 9 and V 9 , which satisfy the zero curvature equation (3.1), are meromor- 
phic functions on the whole complex X plane with the same analytic structure as the 
components U(X) and V(X) of the original Lax connection. 

Proof . First it follows directly from the definitions of g± and 9± that for A on T, 

MA; x, t) ■ \I/(A; x, t) ■ g-\X) = 9_ (A; x, t) ■ *(A; x, t) ■ glHX) 

so that, the two expressions of the right hand side of eq. (3.34) with the + and — signs are 
equal, and effectively define a unique function ^ 9 on T. It is clear that this function can 
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be extended into two functions ^ respectively denned inside and outside this contour 
by: 



\I4 = 0± ■ * ■ gg 1 

These functions have the same essential singularities as * at the points A&. By construc- 
tion, they are such that ^ 9 _~ 1 ^ 9 + \r = 1- 

We may use to define the Lax connection U±, V± inside and ouside the contour T. 
Explicitly: 

u 9 ± = d x v 9 ± ■ v 9 ^ 1 = d x o±ei l + o±uoi l 

V 9 = d t ^f 9 ± ■ = d t 0±0± 1 + 0±V0± 1 

Since V 9 _ V 9 + \r = 1 we see that U + coincides with U- on the contour T and similarly 
V+ = V- for A € T and hence the pairs U±, V± define a conection U 9 , V 9 on the whole 
A-plane. Since 6± are regular in their respective domains of definition, we see that U 9 , 
V 9 have the same singularities as U, V. 

This proposition effectively states that the dressing transformations (3.34) map so- 
lutions of the equations of motion into new solutions. Given a solution U, V of the zero 
curvature equation with the prescribed pole structure and an element of the loop group 
G, we produce a new solution of the zero curvature equation with same analytical struc- 
ture. But since this analytic structure is the main information which specifies the model 
we have produced a new solution of the equations of motion. 

3.5 Soliton solutions. 

In general, a matrix Riemann-Hilbert problem like eq.(3.31) cannot be solved explicitly 
by analytical methods. This statement applies to the fundamental solution of the Rie- 
mann Hilbert problem i.e., the one satisfying the conditions det#± 7^ 0. However, once 
the fundamental solution is known, new solutions "with zeroes" can easily be constructed 
from it. This can be used to produce new solutions to the equations of motion. Starting 
from a trivial vacuum solution, we obtain in this way the so called soliton solutions. 

Let #±(A) be the fundamental solution of a Riemann-Hilbert problem. A solution of 
the Riemann-Hilbert problem with zeroes at ji\ , ■ ■ ■ , /ijV; ^1 ; * * * j is. 









C'(A) = fc'(A)(l-^_£*)x.- 


1 . Pn)xn 
V A — [in J 
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where Pi are projectors Pf = P. We assume that the Pi and the Xi are independent 
of A. When A = m, which we assume inside T, then 0+(A) contains as a factor the 
projector (1 — Pi) so that det 6(fii) = 0. Similarly, if A = Aj, which we assume outside 
r, then 9Z 1 (X) contains as a factor the projector (1 — Pi) so that det#(Aj) = 0. Hence 
the name Riemann-Hilbert problem with zeroes. We now extend the method of dressing 
transformations to the case of a Riemann-Hilbert problem with zeroes. 
Let 



u + — Xn-l 



1 



e (n)-l = ^-1 (A) 1 _ 



A — A n _i 

Ai - Mi 
A - Mi 



A- Ai 



A=//„ 



-L : Xn-l 



A — fJL. 



n—1 



A — A n 



Proposition 18 Given a Lax connection satisfying the zero curvature condition and 
the associated wave function ^(X,x,t), and given vector spaces V n (0), W n (0) we define 



uniquely the projectors 


P n by 










Ker P n (x,t) 


= e W V(X n ,x,t)V n (0) 




(3.37) 












Im P n (x,t) = 






(3.38) 










Then for any g(X) = g 


Z 1 (X)g + (X) on 


T, the transformation *I> — > 







#s = ©i^flfi 1 , ez 1 9+ = y^g-* 

is a dressing transformation, i.e. preserves the analytic structure of the Lax connection. 

Proof . We start with the linear system 

(d x - f/(A, x, t))tf = 0, (d t - V(X, x, t))* = 

and dress it with a solution with zeroes of the Riemann-Hilbert problem, according to 
eqs.(3.35, 3.36): 

U 9 = 6±-U-9g 1 + d x 9± ■ 9±\ V 9 = 9± ■ V ■ eg 1 + d t 9± ■ 9±* 

In general, the components of the dressed connection will have simple poles at the points 
fi n ,X n . We must require that the residues of these poles vanish. At A = /x n , isolating 
the terms containing P n , we have 

MA) ^ M n (i - p B )e£\ 9+\x) ~ ^^e^-'PnM- 1 
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so that 

e + (d x - u)6- + l ~ / ^ An M n (i - p B ) e^\d x - c/| M jef - 1 p n m- 1 

A — [l n 

To kill the pole at A = fj, n , we choose 

ImP n = ef ) ^( /[ /„,x,t)>V(0) 

A vector in Im P n is a linear combination, with coefficients possibly depending on x, t, of 
the column of the matrix in the right hand side. Note that the factor 1 — P n is important 
to kill the terms coming from this x, t dependence of the coefficients. Clearly the same 
argument works for dt — V. Similarly, when A ~ A n we have 

ez^x) ~ e^-^i - P n )N n , e-(\) * ^^iv- 1 P n e { - ] 

A — A n 

so that 

Mas, - u)ez l * ^^n- 1 p n e^\d x - u\ x je^-\i - p n )N n 

A — A n 

To kill the pole we choose 

Ker P n = Im (1 - P n ) = eL n) *(A„, x, t)V(O) 
Clearly the same argument works with dt — V. 

The interest of this procedure is that it yields non trivial results even if the Riemann- 
Hilbert problem is trivial i.e. g(X) = Id. Then its fundamental solution is also trivial 
0±(^) = Id, and the solutions with zeroes are constructed by purely algebraic means. 
The resulting 0±(\) are rational functions of A. 

To make this method effective, we need a simple solution of the zero curvature 
condition dtU — d x V — [V, U] = to start with. Simple solutions can be found in the 
form 

U = A(X,x), V = B(X,t), [A,B] = 

Then ^ = exp ^ Adx + J Q * Bdt^j . The solutions obtained by dressing this simple type 
of solutions by the trivial Riemann-Hilbert problem with zeroes are called soliton solu- 
tions. 



3.5.1 KdV solitons. 

Let us illustrate this construction for one soliton solution of the KdV equation. Recall 
the famous KdV soliton: 



^soliton ( X, t) — 



2k 2 



cosh(fcr + k 3 t) 



(3.39) 
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We obtain it by dressing the vacuum solution u(x, t) = 0. The solution of the vacuum 
linear system reads 



*( M ) = e(^W), 17(A) =(° J) 



We consider the trivial Riemann-Hilbert problem (<?(A) = 1). 
Let 



x _1 (i + 



P 



Hence 6*1 1 (A) has a zero at A = k 2 and #+(A) at A = 0. Since 

, s _ / x\ t/j,2 A — f cos h(A;a: + /c 3 t) /c _1 sinh(/cx + fc 3 i) 
\l/(U,x,ij= I Q Q l, M/(fc ,x,t) = y ks - mh ^ kx + k 3 t - ) ^(kx + k 3 t) 



we can choose 



so that 



Then, we find 



Im P 



, Ker P 



k 1 tanh(fcr + k 3 t) 



p _ . 1 -fc~ 1 tanh(fcx + fc 3 t)\ 
1 ) 



i f — fc tanh(fcr + k 3 t) 



C/soliton_X 1 X- k 2 ktanh(k 



lx + fc 3 t) ) X " 



We see at this stage that the A dependence is already essentially correct. Choosing 

X -- 

we find 



f 

kt<m\\(kx + kH) I 



U soliton 







A - u sohton (x,t) 

with u so \n on (x, t) given by eq.(3.39). We can repeat the procedure, and reach finally the 
N-solitons solution 



d 2 

^N-solitons(^, t) = ~ 2 g^2 lo § T N 



where 



ttv = det(l + W), W t 



V XiX i X . _ 2(kiX+kft) 

h +h- ' 1 ~ 1 

™i ~r rvi 
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3.6 Finite zones solutions. 



In the field theory case, we can use the previous constructions to find particular classes 
of solutions to the field equations, called finite zone solutions. The equations we have to 
solve are the first order differential system: 

{d x -U(\))V = (3.40) 
{dt-V(\))V = (3.41) 

whose compatibility conditions are equivalent to the field equations. The situation is 
very different as compared to the finite dimensional case. As we saw, the analog of the 
spectral curve is 

det(T(A) - n) = (3.42) 

where T(A) is the monodromy matrix of the linear system (3.40,3.41). This equation 
defines an algebraic curve of infinite genus and the analytical tools must be carefully 
adapted. 

If however, we restrict our goal to find only particular solutions to eqs.(3.40, 3.41), 
then we can look at situations where the curve eq.(3.42) is infinitely degenerate leaving 
only a finite genus curve. 

One common way to do that is as follows. The two equations (3.40, 3.41) are exactly 
of type of eq.(??) whose solution was built using the usual analytical tools. 

To interprete the two equations (3.40,3.41) as evolution equations with respect to 
two different "times" for a system with finite number of degrees of freedom we need a 
Lax matrix L(X) satisfying 

[d x - U(X),L(X)} = (3.43) 
[dt — V(A),L(A)] = 

To exhibit such Lax matrices, we consider the higher order flows associated to higher 
Hamiltonians. They provide a family of compatible linear equations (dt { — Vi)^f = for 
i = 1, 2, 3, ■ ■ ■ where we have identified t\ = x, V± = U and ti = t, V2 = V. Since these 
equations are compatible they satisfy a zero-curvature condition: 

Fa = d ti Vj - d tj v - [Vi, Vj] = 0, vi, j = 1, • • • , 00 

We now look for particular solutions which are stationnary for some given time t n , i.e. 
9t n Vi = for all i. The zero curvature conditions F n i = reduce to a system of Lax 
equations: 

^ = [Mi, L], i = 1, • • • , 00 with L = V n , Mi = Vi 

This is an integrable hierarchy for a finite-dimensional dynamical system described by 
the Lax matrix L(X). Taking n larger and larger, the genus of the corresponding spec- 
tral curve usually increases and we get families of solutions involving more and more 
parameters. 
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3.7 The Its-Matveev formula. 



Let us apply these ideas to construct solutions of the KdV equation. As explained, one 
way to get a Lax matrix compatible with the equations of the KdV hierarchy is to seek 
for stationary solutions with respect to some higher time tj. 

A very simple example of this situation occurs when u is stationary with respect to 
the first time t% = t. In that case the Lax matrix is V given in eq.(3.20). The associated 
spectral curve is: 

r : det(V - f i) = i_ l 2 - * A 3 + J(3u 2 - u")\ + i(W - u' 2 - 4u 3 ) = 

The zero-curvature condition becomes the Lax equation d x V = [U, V] and reduces to 
the stationary KdV equation 

6uu - u" = 

Integrating, one gets 3« 2 — u" = C\ and 2u 3 — u' 2 = 2C\u + C2 for some constants C±, 
C<i- So the spectral curve reads 

^ 2 = A 3 /4 - CiA/4 - C2/I6 

It is independent of x as it should be. This is a genus 1 curve and u is given by the 
Weierstrass function. 

u(x) = 2p(x + C) 



For higher times the matrices Vt 2j _ 1 are defined by 

2J-1 , 



(L—)+^> 



d x * 



The matrices Vt 2j _ 1 are not hard to compute but it is quite clear that they are 2 x 
2 traceless matrices. Stationarity with respect to any higher time always lead to an 
hyperelliptic spectral curve and we will just retain this feature. 



Hence we consider an hyperelliptic curve of the generic form: 

2g+l 

r : /x 2 = R{\) = \{ (A - A,) 
i=i 



(3.44) 



The point at co is a branch point, and a local parameter around that point is z = \f\. 
Our goal is to construct a function ^ on T satisfying the equations 



M/ = (3 2 -u)* = A*, d t2 .^ 



(3.45) 
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for some potential u. If we succeed to do it, u(x, £3, £5, • • •) will have to satisfy the KdV 
hierarchy equations by consistency 

Baker-Akhiezer functions are special functions with essential singularities on Rie- 
mann surfaces. It is a fact that there exists a unique function $onT with the following 
analytic properties 

• It has an essential singularity at the point P at infinity: 

V(t, z) = (l + ^ + 0(V* 2 )) (3-46) 

where z = VX and £(i, z) = ^>i z 2l-1 *2i-i- 

• It has g simple poles, independent of all times. The divisor of these poles is denoted 
D = (71, •••,7s)- 

This function is called a Baker-Akhiezer function. Even though Baker-Akhiezer 
functions are not meromorphic functions, they have the same number of poles and zeroes. 
In fact the differential form d(log *) is a meromorphic form. The sum of its residues is 
the number of zeroes minus the number of poles of and this has to vanish. The essential 
singularity does not contribute because around P we have d(log^) = d£(t, z) + regular 
and d£(t, z) has no residue (remember that 1/z is the local parameter around P). 

The uniqueness is then clear. If we have two such functions with the same singularity 
structure and divisor D, their ratio is a meromorphic function on T with g poles which 
can only be a constant. 

The existence will be proved by giving an explicit formula for ^ in terms of 6 func- 
tions. 

But before that, we will prove that this Baker-Akhiezer function solves the KdV 
hierarchy equations. Let ^> be the Baker-Akhiezer function as above with the following 
behaviour at infinity: 



00 

* = e^ z \l + 0(z- 1 )), £(t,z) = 2> 2j ' -1 *2i-i 



Then we have 

Proposition 19 There exists a function u(x, t) such that 

(d 2 x - u) * = 
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(3.47) 



Proof . Consider on T the function d 2 ^f — A*. To define this object as a function on the 
curve, A is viewed as a meromorphic function on T. Remark that A has only a double 
pole at oo and such a function exists only if T is hyperelliptic. We see that d 2 ^> — A^ 
has the same analytical properties as \F itself at finite distance on T. At infinity we have 
by eq.(3.46): 

efci - A^ = e^ t ' z \2d x a + 0(l/z)), z = v 7 ! 

So it is a Baker-Akhiezer function, but with a normalization 2d x a instead of 1 at infinity. 
By the uniqueness theorem of such functions, we have: 

d 2 x <S> - A* = u*, u = 2d x a (3.48) 

Having found the potential u, we construct the differential operator L = d 2 — u and 
show that the Baker-Akhiezer function ^ obeys all the equations of the associated KdV 
hierarchy. 

Proposition 20 The evolution of * is given by: 

d t2i _^ = (L^) + * 
where L = d 2 — u is the KdV operator constructed above. 

2i — l 

Proof . Consider the function dt 2i _ 1 i & — {L 2 )+'$?■ It has the same analytical properties 
as ^ at finite distance on T. At infinity we have d t2i _^ = z 2 ^ 1 ^ + e^O(l/z) and 

{L^r) + -$ = L^*- (L^r k )_<I' = z 2 ^ 1 ^ + e^O(l/z), where we have used = z 2 ^. 
Hence we get: 

9 i2 ,_ 1 $-(L 2! r) + $ = e* 2 )0( 2 - 1 ) z -► oo 
By unicity, this Baker-Akhiezer function which vanishes at oo vanishes identically. 

We now give a fundamental formula expressing the Baker-Akhiezer functions in terms 
of Riemann theta functions. Recall that a differential of second kind is a meromorphic 
differential with poles of order > 2. Let f^ 2 -' -1 *' be the unique normalized second kind 
differential (all the a-periods vanish) with a pole of order 2j at infinity, such that: 

q(2j-i) _ ^ ( z 2 i-! _|_ regular) , for z — ► oo 
Let U^~^ be its 6-periods: 

Define 

3 

and denote by 2iirU the vector of 6-periods of Q. 
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Proposition 21 If D = Ylf=i 7« ^ s a generic divisor of degree g, the following expression 
defines a Baker- Akhiezer function with D as divisor of poles: 



*(P) = const, exp ft 



6{A{P) + U-Q 



Po ) 9(A(P)-Q 

(3.49) 



Here ( = A(D) +/C with K. is the vector of Riemann's constants and A denotes the Abel 
map with based point Pq. 

Proof . It is enough to check that the function defined by the formula (3.49) is well- 
defined (i.e., it does not depend on the path of integration between Po and P) and has 
the desired analytical properties. Indeed when P describes some o-cycle, nothing hap- 
pens because the theta functions are a-periodic and f2 is normalized. If P describes the 
fy- cycle the quotient of theta functions is multiplied by exp(—2iirUj while the exponen- 
tial factor changes by exp(2i7rf/ ? ), so that ^ is well-defined. Clearly it has the right 
poles if deg D = g. ■ 



As a consequence, the normalized Baker- Akhiezer function with the divisor of poles 
D = (71, • • • , 7 S ) can be expressed as: 



d (A(P)-oe(u-o 



(3.50) 

where A(P) is the Abel map on T with base point 00, and £ = A(D) + K. with fC the 
Riemann's constant vector. 



The KdV potential, u, is given by the Its-Matveev formula: 



u(x,t) = -2d 2 x \oge(^2t 2j -iU {2j ~ 1) - + const. 

(3.51) 

In fact, in eq.(3.50) the integral has to be understood in the following sense: for z in a 
vicinity of 00, one defines f^ 2 - 7-1 ) as the unique primitive of r^ 2 -? -1 ) which behaves as 
z 2 ^ 1 + 0(l/z) (no constant term). Of course, when this is analytically continued on the 
Riemann surface, 6-periods will appear. However they will cancel out in eq.(3.50) due 
to the monodromy properties of ^-functions, leaving us with a well-defined normalized 
Baker-Akhiezer function. The formula for the KdV field is found by using: 



X + u = (d 2 x log*) + (8 x log^ 
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Setting fit 1 ) (2) = d(z + § + 0(z 2 )), where /? does not depend on times, we have: 

a x log* = z + ^iog^(^(P) + ^t 2j _ 1 ^- 1 )-c) - 
-a, log ^X;^--!^'- 15 - c) + - z + o(^ 2 ) 



We evaluate this expression when z — > 00. Using Riemann's bilinear identities, we can 
expand the Abel map A(P) around 00, and we have: 

- « ((*->" + («»- a?) <"■>+ 



-c 

Keeping the l/z terms, we obtain: 

d x log * = z - ^ log ( t 2j -iUto-V " C) + f + O(^) 

j 

Differentiating once more with respect to x, we also get c^log* = 0(1/ z). It follows 
that z 2 + u = z 2 - 2d 2 log 6 + 2j3 + 0(1/ z) proving the result. 
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Chapter 4 

The Jaynes-Cummings-Gaudin 
model. 



4.1 Physical context. 
4.1.1 Rabi oscillations. 

This is the interaction of a two levels atom and photons in a single cavity mode. When 
the electromagnetic field is classical, the system is described by the Hamiltonian 

H = /kj y + hft[a + b + cj-tf] 




We assume 

b(t) = e- iu) % 

The resonnance condition is A = ujq — uj = but we keep A / for a while. Let us 
denote 

|^ = e -uoblb t ( ^1 



i>1 



The Schroedinger equation becomes 
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From this we get a second order equation for tp\ 



V>i - iuipi + ^ 2 6q&o + |(^o - 2cj cj)J ^i = 
If the atom is in its fundamental state at time t = 0, i.e. ipi(0) = 0, then the solution is 

^ = ae -*f V 7 * - e" i7i ), Y>2 = -7^7-e i t*((2 7 + A)e** + (2 7 - A)e~^) 

2S Zoq 



1 



where we introduced the Rabi frequency 



fiLbi = nV + 6j6o), « 



A 
2ft 



The norm of the state is = + ^2^2 • The probability to find the atom in the 

excited state at time t is 

p m W(t)Mt) 2A h sin 2 Rabi t 

Pe(t) = ^#r = ^° 6o ^r 

These are the famous Rabi oscillations, the amplitude is maximal at the resonnance 
A = 0. Notice that 

{ip\a z \ip) _ y^i - ^2^2 = 2 i^h _ 1 



so that 



(4.1) 



What happens if the electromagnetic field is quantum? and in particular if the 
number of photons is small (5 < n < 40)? The system now is described by the Jaynes- 
Cummings Hamiltonian [1] 

H = hioo^- + hwtfb + Ml[a + b + a~b^] 

Where we recall the usual commutation relations 

[a z , a 1 * 1 ] = ±2a ± , [a + , a~] = a z , [b, tf] = 1 

It turns out that the model is still exactly solvable. The key is the existence of an 
extra conserved quantity. Let 

^ = tfb + X -a\ H = ko z + &V- + a + b, k = — 



We have 



H = toHi + QH , [H x ,H o ] = 
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The Heisenberg equations of motion are 

is z = [H, s z ] =n(H - 2ks z - 2bs~ 

But for the spin 1/2 representation, we have 

k b 



H 



so that 



Hence, we get 



s z = in 



(k - 2H s z ^j 



Since Hq is conserved the solutions is extremely simple (keep the order of the operators) 



s *(t) = ^H' 1 + e- 2mt H ° 



(«*(0) " ^ 



We now introduce the states 

|n,T) = (6 t r|0)®|T), M) = (6t)"|0)®||) 

We have 

H \n, t) = n\n, |) + |n + 1, j), i?oK I) = "|n - 1, |) - re|n, |) 
which implies 



Let us define 
It is then simple to show that 



Hl\nA) = {K 2 + n + l)M> 

n 2 n = n 2 ( K 2 + n+i) 



(m,T |s*Kt) = 2<WiH™) 



1 -2ft 2 (n + l) 



sin 2 O n i 
fi2 



We introduce the coherent states 6(0) | a) = a|a) 



I|a| 2 ^ « 



|n) — ( - " — —\n 



n=0 
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For such states, the average number of photons in the cavity is n = \a\ 2 . We obtain 



n=0 



(a, I |<r z (0)|a,|> 



n! 



Drawing this quantity as a function of time, we observe a collapse and resurgence phe- 
nomenon of the Rabi oscillations. 




Figure 4.1: The collapses and revivals of Rabi oscillations, n = 30, A = 2\/2, fi = 1. 



4.1.2 Cold atoms condensates. 

Consider alcali atoms like Li, K, Na, etc... Let / denotes the magnetic moment of the 
nucleus and S the spin of the electron. By choosing the isotope, we can arrange that the 
atom is a fermion or a boson. We will consider the case of a fermion. The Hamiltonian 
of a single atom in a magnetic field is 



H = gl ■ S + g B S ■ B 

For two atoms far apart, the Hamiltonian is simply the sum of the Hamiltonians of the 
idividual atoms. 

H = H 1 + H 2 

When they come closer together however they start to interact. We consider a situation 
where we can have formation bound states or "molecules" between two atoms. Atoms 
are fermions, molecules are bosons. 
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atom atom 



molecule 




Let jc an d Cj a be creation and annihilation operators for fermions in the state a =| or 
I in an orbital of energy ej. Let and b be creation annihilation operators of a molecule 
at zero momentum. The Hamiltonian of the boson-fermion condensate is 

This can be rewritten in terms of pseudo spins 

CT 

we get 

ra— 1 rt— 1 

3=0 3=0 

This is the Jaynes-Cummings-Gaudin Hamiltonian. 

In the Born-Openheimer approximation the energy levels become a function of the 
distance E — ► i£(r) as shown in Fig. [4.2]. A Feshbach resonnance occurs when a bound 
state becomes degenerate with a scattering state. By tuning the magnetic field, one can 
adjust the molecule state to be just above or below the atomic state. We can thus induce 
a transition in the system at will. A particularly interesting situation is the case of a 
soudain perturbation. At t = the system is in an atomic state, and at t > molecules 
start forming in the fundamental state (at zero temperature). What is the dynamical 
evolution of the system? What is the rate of formation of "Cooper" pairs? 



(4.4) 



4.2 Settings. 

Let b, and Sj be quantum operators 

[b, &t] = h, [*+ sj] = 2hs z , [s z j, sf] = ±hsf 
We assume that Sj acts on a spin s representation 

Sj\mj) = hmj\mj) , s^\rrij) = hyj 's(s + 1) — mj(rrij ± 1) |mj±l), rrij = — s, — ■ ■ • , s— 1, . 
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V(r 





Figure 4.2: Feshbach resonnance. 



where s is integer or half integer. Notice that 



so that 



(- + ) P l-> = J^0il- + r) 



(8+Y I - s)\\ 2 = 2slrl = r(25 + 1)r(r + 1) (4 5) 

15 j 1 8)11 (2s-r)\ r(2a-r + l) [ ' 



Hence, if r > 2s the norm is automatically zero because the T function in the denomi- 
nator has a pole. Similarly, we assume that b, act on the Fock space tf n \0). 

Instead of the representations above, we will work with the Bargman spaces. For the 
oscillator b, 6' this is the space 



Bb = entire function of z J |/( 



(z)\ 2 e ' h dzdz < oo 



On this space we have 

b = h^-, tf = z (4.6) 
dz 

For the spin operators, following Sklyanin [5], we set 

s z = h [ w— s 



dw 

s + = hw 

d 2 d 
s = n[-w-—^ + 2s— 
dw z dw 
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Notice that the value of the Casimir is 

(s z ) 2 + \( s+s ~ + s ~ s+ ) = h2s ( s + 1) 

The lowest weight vector corresponds to the constant function 1. Other states are 
obtained by applying s + ~ w. Let 1 be the ideal in the set of polynomials C[w] generated 
by w 2s+1 . The above operators are well defined on C[w]/I. This amounts to showing 
that s a I C X. We have only to check the dangerous case 

s~w 2s+1 = -(2s + l)2sw 2s - 2s(2s + l)w 2s = 

The representation space is C[w\/X and is of dimension 2s + 1. 



4.3 Bethe Ansatz. 



On this Hilbert space acts the Hamiltonian 

n— 1 n— 1 



Let 



where 



3=0 3=0 



L(X) _(A(X) B(X)\ 



2A w ^ s 2 



C(A) = - + V^- 

are now quantum operators. It is very simple to check that 

[A(A),A( M )] = 

[B(A),B( M )] = o 

[C(A),C( M )] = o 

[A(A),B( M )] = ^-(BCAJ-B^)) 

A — /i 

[A(A),C( M )] = -^(C(A)-C(m)) 
[B(A),C( M )] = ^(A(A)-A(a*)) 
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One can rewrite these equations in the usual form 



[L 1 (\),L 2 (Li)} = -h 



P12 
A — fi 



,Li(A) + L 2 (/x) 



where 



We now show that 



/l 0\ 

D _ 10 

12 ~ 10 

\0 1/ 



Tr (L 2 (A)) = 2A 2 (A) + B(A)C(A) + C(A)B(A) 



generate commuting quantities. 
Proposition 22 we have 



[TrL 2 (A),TYL 2 Gu)]=0 



Proof. First one has 



with 



[IV(L 2 (A)),L(/x)] q =[M(A, / x),L(/x)] £ 



M(A,/z) = -2h 



L(A)-L( M ) 



A — fi 

where we distinguished the commutators in the quantum space and the commutator in 
the auxiliary space. Alternatively, we have 



[Tr(L 2 (A)),L 2 (^)] q = 

It follows that 

[Tr(L 2 (A)),L 2 (^)] q 

and therefore 

[Tr(L 2 (A)),Tr(L 2 ( M ))] q 
2h 2 



2h 

A — fi 



[L 2 (A),L 2 (//)] £ 



2h 
A — fi 



Tri [P12L1 ( A) , L 2 (//)] (4.7) 



2h 
A — fj, 



[L 2 (A),L 2 (^)] £ 



2h 
A — 11 



TVipx.L^A)^ 2 ^)] 



2h 
A — /x 



Tr^P^L^A),^)] 



— ^Tr 12 P 12 ([P 12 , Li(A) + L 2 (fi)]L 2 (fi) + L 2 (/x)[P 12 , L : (A) + L 2 (/x)]) 
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(A-m) 



Expanding the four terms and using P 2 2 = 1 and the cyclicity of the trace (for P± 2 only) 
we arrive at 

[Tr (L 2 (A)), Tr (L 2 (/x))] = ^— -^Tr 12 {[U(X), L 2 (/x)] - P 12 MA), L 2 (/,)]) 

The first term in the right hand side vanishes because Tr L(X) = 0. The second term 
vanishes too because 

Ti 12 P 12 [L 1 (X),L 2 ( f x)} = — ^TrP 12 [P 12 ,L 1 (A)+L 2 (^)] 

A — \X 



h 



A — \i 



Tr(L 1 (A) + L 2 (/i))[P 12 ,P 12 ] =0 



We are now in a position to write the Bethe Ansatz. Let 

|0> = |0> <8) | — si> <8) ■ ■ ■ «8) | — s„), 6|0) = 0, sj\0,- Sj ) = 

We have 



B(A)|0> = 



and 



A(A)|0> = a(A)|0> 



, . . 2A uj v-^ hsj 
g 2 g 2 Y A_ei 



Since 



we also have 



With all this we have 



[5(A), C(A)] = 2hA'(X) 



B(X)C(X)\0) = 2ha'(X)\0) 



2 TrL 2 (A)|0) = (a 2 (A) + fe'(A))|0) 



Let 



^2, • • • , mm) = c(m)C(n 2 ) ■ ■ ■ c([i M )\o) 



97 



For the following, we need the 



[TrL 2 (A),C(^)] = - 4 ^(C(/x)A(A) - C(A)A(/x)) 

which is obtained from eq.(4.7) by pushing the to the right. We recall also that 

L4(A),C( M )] = -^-(C(X) - C{p)) (4.8) 

Proposition 23 One has 

^TrL 2 (A)f2(/ii, A*2,- • • , Mm) = A(A, //i, ■ ■ ■ , mm)^(Mi> M2, ■ ■ ■ , Mm) (4-9) 
+ X] A ^ A ' ^' " ' ' Mm)^(A, m, fi M ) (4.10) 

i 

T7ie /irsi term is called the wanted term, and the other ones are called the unwanted 
terms. Their coefficients are respectively 

A(A, mi, ■ ■ • , Mm) = « 2 (A) + W(A) + £ Y ^a(A) + 2 £ (4.11) 

Ai(A, mi, ■ ■ ■ , Mm) = t— ~~ «(Mi) + J" — " — 1 (4-12) 
A ~ Mi \ r^T Mi - Mj ' 

Proof . We have 

1, lrL"(A)S2(/xi,/x 2 , ■ ■ ■ ,Mm) = 

1 



2 TrL 2 (A)fi( /[ /i, A i 2 ,---,MM) = (a 2 (A) + ^a'(A))^(^i,^ 2 ,---,MM) (4.13) 



+ -[TrL 2 (A),C( m )C(^)---C( / u M )]|0) 
But 

^[TrL 2 (A),C( Ml )C( / u 2 )---C( m )]|0) = (4.14) 
= E T^-CVi) • • • (C(M0^(A) - C(A)i( tt ))c( W+1 ) • • • C(/iAf) |0) 

I 

We now push -A(A) and ^4(Mi) to the right, using eq.(4.8). Clearly, when we do so we 
will generate terms only of the form (remember that |0) is an eigenvector of both .A(A) 
and A(m)) 

0(mi,M2,---,Mm) = C(/xi)C(// 2 )---C(^Af)|0) 

n(A,//i, ■ ■ ■ ■ ■ ■ ,hm) = c{\)C{ni)---c{m)---c{nM)\Q) 
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The wanted term cannot come from the term C(X)A((j,j) in the above formula because 
one of the C(/j,i) has already replaced its argument by A and there is no way to recover 
it. Hence one has to use the first term C(/j,i)A(X) and push A(X) to the right using only 
the second term in eq.(4.8). We get in this unique way 



h 



fi(//l,/i 2 , • • ■ ,/iAf) 



Adding the wanted contribution, eq.(4.13), we obtain eq.(4.11). Let us see now how to 
get the first unwanted term, the one where fi\ has been replaced by A. Clearly this term 
has to come from the term C{X)A{m) in eq.(4.14). Then one has to push A(m) to the 
right using only the second term in eq.(4.8). We get in this unique way 



2h 



h 



fi(A,// 2 , • • • ,Mm) 



Since fi2, • • • , Atjw) is completely symmetrical in the //j, we have proved eq.(4.12). 
The unwanted terms vanish if the Bethe equations are satisfied 




(4.15) 



Taking into account these conditions, the eigenvalue can be rewritten as 




(4.16) 



4.4 Riccati equation. 

We now analyse the Bethe equations eqs.(6.8). We introduce the function 




Proposition 24 The Bethe equations (6.8) imply the following Riccati equation on S(z) 

2 

%~ 2 



S'{z) + S 2 {z) + ^({2z-uj)S{z)-2M) =Y,2 Sj S{z) ^ 



Z — €j 



(4.17) 
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Proof . The Bethe equations read 



2fii lo t—^ hsj s—^ h 
g 2 g 2 ym-tj j^fK-Vj 

we multiply by l/(z — Hi) to get 

2 fli LO 1 y. hSj 1 H 1 q 

g 2 z- fii g 2 z- ^ ^ m - €j z - m ^f. Hi- Hj z- Hi 
We now sum over i. We have 

M 

jr = y + = — m + ^(z) 

^ z - Hi t-r 1 z - Hi z- Hi 

t=i i 

y^ 1 _ 1 y( 1 , 1 \_ S(z)-S(e j ) 

^ Hi - €j z - Hi z - ej ^ \z - Hi Hi~ e jJ z~ e j 

yy^ L_ = ly*y^ 1 / 1 1 

if y 1 = I(y 1 

= \{S 2 {z) + S'{z)) 



In equation (4.17) the S(ej) appear as parameters. They can be determined as 
follows. Suppose first that Sj = 1/2. We let z — > e« into eq.(4.17) getting 



S'(ei) + S 2 (e t ) + * ((26, - w)^) - 2Af) = S'(e.) + £ S(ej) 



The remarkable thing is that S'(ei) cancel in this equation and we get a set of closed 
algebraic equations determining the S(ej). 

Proposition 25 Let Sj = 1/2. In that case the constants S(ej) are determined by the 
set of closed algebraic equations: 



S 2 ( £l ) + A ((2* - u)S(ei) -2M)=J2 S(et) S(e '\ • = 1, ■ ■ ■ , n 
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Suppose next that s = 1. We expand the Riccati equation around z = e^: 



(z-ei)° : 5'(e l ) + 5 2 (e J ) + -^((2e l -u;)5(e l )-2M) =2 < S , (e,)+2^ — 

5" „•_*.,.■ e * e i 



(z-e;) 1 : 5"( ei )+25(e i )5 , (e i ) + ^((2e l - W ) < S , (e J ) + 25(e l )) 

= S"{e l )-2Y J 



S(€j) - 5(6,-) S'(€j) 

(e; - e,) 2 e< - 



We see that in the second equation 5"'(ej) cancel. The first equation allows to compute 
S'(ei) and the second equation then gives a set of closed equations for the S(ei). We 
have shown the 

Proposition 26 Let Sj = 1. In that case the constants S'(ej) and S(ej) are determined 
by the of closed algebraic equations: 



5 2 (e i ) + F ^((2e i -o;)5(e i )-2M) = 5 , (e l )+2^ 

5 „•_£.• £ * £ i 



25(e l )5'(e l ) + ^((2e,-^)5'(e l ) + 25(e,)) = _ 2 £ _ 

The general mechanism is clear. For a spin s, we expand 

_ 2 /(^)~%) = y rn^2s s{m){e){z _ e)m -i 
z-e ^ ml 

m 

and we see that the coefficient of S^ 2s \e) vanishes in the term m = 2s. The equations 
coming from (z — e) m_1 for m = 1, ■ ■ ■ , 2s — 1 allow to compute 

S'(e),---,S^- 1 \e) 

by solving at each stage a linear equation. Plugging into the equation for m = 2s , we 
obtain a closed equation of degree 2s + 1 for 5(e) . 

P 2s +i(5(e)) = (4.18) 

Notice that if M < 2s, the system will truncate at level M because there always 
exists a relation of the form = P(S, S', ■ ■ ■ S^^). 

The S(ej) also determine the eigenvalues as well. Recall that 

A(A, Ml , • • • , M m) = a 2 (A) + ha'(X) + 27^ ° (A) ~ ^ 
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so that 



hence 



and 



g(A) - a(fij) 
A- m 



2 ^ — ^ Hs j 



9" ^\-{J>i\\-ej fii-ej 
2 ^ hsj 1 

tr j x ~ e j I' 1 ~ 'j 



a(A) - a(ni) _ 2M x - hsjS(ej) 
X- Hi g 2 X-ej 



A(A) = a 2 (A) + fia'(A) + 2ft ^ - £ j 



Expanding a(A) we find 



ft 2 Si ( Sj + 1) 



A(A) = i(2A- W ) 2 + V« + ^EA + E , 



with 



i? n = KM - hy^ sj + 



and 



9 9 



■ e j ~ e i 



4.5 Baxter Equation. 

We can linearize this Riccati equation by setting 



S(z) 



Obviously 



M 



4>(z) = Y[(z- im) 



The linearized equation reads 



2 / 2M v-^ HsjSUj) . 



ip"(z) + -a(z)if}'(z) + , , ,. , , 

n \ g z *-r* z — e 3 
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Here, we should understand that S(ej) are determined by the procedure explained in 
the previous section. For such values of the parameters, the equation has the following 
remarkable property. 

Proposition 27 For the special values of the parameters S(ej) coming from the Bethe 
equations, the solutions of eq. (4-22) have trivial monodromy. 

Proof . This is clear for the solution 4>\{z) defined by eq.(4.21) since it is a polynomial. 
A second solution can be constructed as usual 



ilo/ 



(z) J" exp (-1 j V a(t)dt-2]ogMv))dV = Mz) f 



\2.s 



e~^ 2 - 



V2(z)- V { ^ , . , ,. V . VM „. v LV .., , |r(y _^ 2 

The monodromy will be trivial if the pole at y = fii has no residue preventing the 
apparition of logarithms. Expanding around /Zj, we have 

rv \ e (-sf'«(<Mi-2E # ,(ft-«)) 

/ a(t)dt-2]Qgii> 1 (v)] = 



exp 



2 

h 



(y - 



exp ( -jXy-in) 



+ 0{z-mf 



but the coefficient of the dangerous {y — Hi) term vanishes by virtue of the Bethe equa- 
tions. ■ 



Next we set 



4>(z) = exp ( / a{y)dy ) Q(z) 



We obtain for Q(z) the equation 



h 2 Q"(z) - [ a\z) + ha'(z) + ^ - 2ft 2 £ | Q{z 



Comparing with eq.(4.19), we obtain Baxter's equation 



h 2 Q"{z) - A(z)Q(z) 



Notice that 



so that 



exp 



a{y)dy)=e-^ {z2 - wz) \{{z-e j y 



-Kj(z 2 -luz) -K i (z 2 -ojz) M 

Q[z) = m-^ iz) = wr^y " 



(4.23) 



(4.24) 



f[(z-tH) (4.25) 
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4.6 Bethe eigenvectors and separated variables. 

We recall that 

fi(/zi,/i 2 , • • ■ ,hm) = c(m) ■ ■ ■ c(fi M )\o) 

By definition of the separated variables, we have 



C(X) = 



2^ nr =1 (A-A.) 

9 m=i(A-Ci) 



Inserting into the Bethe state and remembering that 

il>(z) = Yl( z ~ A*i) 

we get 



n(Mi,/i2,---,/iji#)= (n^y (7 j n^)i°) 



Proposition 28 In the separated variables, the Hamiltonians read 

rr _ Ylk^-^k) V n ¥ A-^) ( d 2 2 d 

1 U k ^i - e k ) 2^ n M ,(A, - A fe ) [dXj + h { ]) d\ 3 h 
Proof . Write eq.(4.22) for each separated variable as 



V- 1 TT I / \ \ ( ^ 2 d M \ 

^A^^ ) = -(rfA| + ^ (A ^"^J 



(4.26) 



Since this formula holds for a basis of eigenvectors, we can "factor" by 4>{\i). Inverting 
the Cauchy matrix = l/(Aj — 6j), and taking care of the order of operators we obtain 
the Hamiltonians Hi in terms of the separated variables 



1 _ 2 ^_ _ M 

« VJ ' " %3 "A, - £j ' Kj " dX* + h a[Xj) d\ 3 h 



Hi - -B^Vj, B 



explicitly, they are just eqs.(4.26). 



Proposition 29 The Hamiltonians eq.(4-26) commute. 
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This is a general result and we will prove it in section 4.8. 
We now introduce the scalar product 

\M\ 2 = JdXidXi WW p(x 1 ,x 2 ,---,x n )\^(X i )\ 2 (4.27) 

where 

W = H(x t ~ Xj) 

and 

x _ i Ujh-Xj] 2 
h Ylk^i( e i- e k) 2 

The rationale for this is that we want to change variables from the spin variables s+ = 
hwj to the separated variables Xj, that is 



2z n.(gj-Ai) _ _zz 



The factor WW just comes from the Jacobian of this transformation. In the reduced 
model however, the variables z,z are integrated out. The measure p(x±,X2, ■ ■ ■ ,x n ) is 
determined by requiring that the Hamiltonian Hj are Hermitian. 

Proposition 30 The Hamiltonians Hj are Hermitian with respect to the scalar product 
eq.(4.27) with 

P(xi, x 2 , = J" dye-yy^-ZM+W J"J J - 2s ^0 ) (4.28) 

where J 2si+ i(x) is the Bessel function. For n = 1, the formula for p{x) can be simplified 
giving 

p(x) = d™ +1 [e-*x M - 2s ] = e-P M -2 S (x) 
where Pm-2s{x) is a Laguerre polynomial of degree M — 2s. 

Proof. We have to show that 



/ 



dX k dX k mX k )\ 2 W-il + yL a(Xj) + M j B^\ W \ 2 p( Xl , ...,x n ) (4.29) 



is real. Now 



where Ajj is the minor of the element Bji. It is clearly independent of Xj. Hence 

d -B~ 1 = B- 1 (—- A^ 1 — A 



dXj 13 13 \dXj dXi 
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We have 

a = rw^rw^) = n _ tj) . lw ^ 

1M A< e i> j+i i= i 

where we introduced 

rij( 6 i - ^j) _ zjZj 

Remark the important formula 

d 

-JyZk = Bj k Z k 

hence 

d . „, i d 



so that 



Next, we have 



-^-p(x U ■ ■ ■ ,X n ) = ^2 B jkX k -7^p{xi, • • • ,X n ) 



d 2 d 2 
-^p(x lr --,x n ) = B jkBjix k xi Qx Qx p(xi, • • • , x n ) 

3 k,l 

32 



E n 2 ^ / \ v — > 1 5 



fc k k,l 

Putting everything together eq.(4.29) becomes 



J d\ k d\ k \^{\ k )\ 2 Y^B^\W\ 2 \-Y,B 2 jk x k V k + Y,B jk O k + \ P(xi, 



where 



V k = x k d Xk + 2(s k + l)d Xk 
V ° = (j2 Xkdx * + M + n + 1 j 
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and 

O k = \( e k-^)-^Y^ — ~ {( x k d x k + s k + l)(xid Xl + s l + l)-a fc sj-l) 

l^k €k el 

The conditions on p(x±, ■ ■ ■ ,x n ) are that the quantity in the curly bracket in eq.(4.30) 
should be equal to its complex conjugate. When we do the sum over j, we have 

which is real and gives no condition. Next we have the identities 

j 



j 

The conditions on p(xi, ■ ■ ■ , x n ) then read 



kj=-i 



-(zt - Zi) [V iP + V p] + £ Zk 2 _ *** [V iP - V kP ) = 
Finally find the n conditions 



Vip-V k p = 0, k^i (4.31) 

Vip + Vop = (4.32) 

Notice that eq.(4.32) is independent of % if the conditions eq.(4.31) are satisfied. A 
solution of eq.(4.31) is 



p{x lr -- ,x n ) = J2 C p n 



xf 



P=0 <JlH h<2n=pi = l V ' 

Then eq.(4.32) gives 

Cp+i + (M + n + p + 1)C P = 

the solution of which is 



c p =(-i)»( m+ ; +p )p! 



Hence we have found 



P (x 1 ,x 2 ,..,x„ ) =B-i)'( M+ ; +J ')f! e n „, (2 ..; 1+ ,.), (« 3 > 

p=0 vp/ «i+-+9 n = P <=i ft -^ s * + 1 + ft ''- 



107 



This is just eq.(4.28) 



This important formula should be further studied. In particular, for ip{z) being a 
Bethe state, one should be able to compute it exactly because we know that by Gaudin 
formula 

||n(/xi,Ai2,---,MM)|| 2 = M!detA 
where A is the Jacobian matrix of Bethe's equations. This is still very mysterious. 



4.7 Quasi-Classical limit. 

The exact formula relating Q(z) and ip(z), eq.(4.25), allows to study the properties of 
the solutions of Bethe roots \i\ in the quasi-classical limit h — > 0. 
Let us set 

v{z) = h W) 

Then Baxter's equation, eq.(4.36), becomes 

hy'(z)+y 2 (z)=A(z) (4.34) 

where we recall that 

In the semi-classical limit eq.(4.34) becomes the equation of the spectral curve of the 
model (in that limit hsj = O(h )): 

y\z)=K{z) 

This also means that we can write in the semi-classical limit, as expected, 

Q(z) = e K / z vWdX ^ e \ f z VW)dX+O(h0) 

From eq.(4.23) we deduce that 

h 



y{z) = a{z) + ^ 



z — \i 



so that we expect in the semi-classical limit 

- \fW) ~ a ( z ) 

Z f-ii 

i 

This is a remarkable formula. It gives us the distribution of Bethe roots m in the 
semi-classical limit, as we now show. 
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Let y/ A(z) be represented as a meromorphic function in the cut z- plane. Let us put 
the cuts so that 

y/Mzj = 4 - 4 +0(z- 1 ), \z\^oo 
g 2 9 2 

and (we neglect terms of order h which do not contribute in the leading h approximation) . 

V / AM= f^ + 0(l), z^e j 

Z 6j 

By the Cauchy theorem, we have 

dz' y/AjFj 



C 2m Z 1 — Z 



where C is a big circle Co at infinity, minus small circles Cj around z = €j , minus contours 
Ai around the cuts of y/A(z). Hence 



r— f dz' y/Xjz 7 ) [ dz' y/Ajjj f dz' y[U/) 

V {Z> J Co 2m z'-z L> J c . 2m z> - z ^ J A . 2m z> - z 



J 

But 



dz' y/Ajz 7 ) _ { /T7 -^ _2z uj 

Jc 

and 



c 2i7T z> - z v """"" u 1 g 1 

dz' y/A(z') f dz' —hsj hsj 



, c . 2m z'-z J c . 2m (z 1 - ej)(z' - z) 
so that we arrive at 



/ 2z uj sr-^ frsj sr-^ f dz' y/ A(z') 

V{Z) = 7~7~^ ~ 4- J Ai 2m~^^7 



hence 



and therefore 



y/Ajz) ~ <Z) = Y, J 



dz' t/AJz 1 ) 
a. 2m z' — z 



v — = v f ^y^l + o(H) 

4^ z-m ^J Ai 2m z-z' 

Comparing both members of this formula suggests that the Bethe roots m accumulate 
in the semiclassical limit on curves Ai along which the singularities of both side should 
match. To determine these curves we assume that the Bethe roots \ii tend to a continuous 
function fj,(t) when h^0(t = hi and i = Oift^ 1 )). 



^ ^ f [ d ( 



dt\ 1 



dfj, J z — [i 
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Here A = J2^i- Hence, comparing with the semi-classical result, we conclude that the 
function fi(t) should satisfy the differential equation (we rename the function fi(t) to 

*(*)■) 

? = (4-35) 

The boundary condition is that the integral curve z(t) should start (and end !) at a 
branch point of the spectral curve y 2 = A(z). 

This result can be checked by numerical calculation. We consider the one spin-s 
system. A typical situation is shown in Fig. (4.3). The agreement is spectacular. 

We can say a word on how the Bethe equations were solved. We first determine 5(e) 
by solving the polynomial equation eq.(4.18) and then determine ip(z), eq.(4.21), by 
solving eq.(4.22). The Bethe roots are then obtained by solving the polynomial equation 
il)(z) = 0. 



4.8 Riemann surfaces and quantum integrability. 

Let us consider a set of separated variables 

[Aj, Xj] = 0, \pi, fij] = 0, [Xi,fij] = p(Xi, m)5ij 
We want Baxter's equation, so we start from the linear system 

^2 Rj(*i, Vi)Hj + Ro(X t , m) = (4.36) 

3 

Here the Hj are on the right, and in Rj(Xi, Hi), Ro(Xi, Hi), we assume some order between 
Xi,fii, but the coefficients in these functions are non dynamical. Hence we start from 
the linear system 

BH = -V (4.37) 

We notice that we can define unambiguously the left inverse of B. First, the determinant 
D of B is well defined because it never involves a product of elements on the same line. 
The same is true for the cofactor A^- of the element By (we include the sign (— \) l+] in 
the definition of Ajj). Define 

= {B-% = D-'Aji 

We have 

(B~ 1 B) i j = ^ D~ 1 A ki B k j 
k 

But Afcj does not contain any element Bki, hence the product A^B^ is commutative, 
and the usual construction of the inverse of B is still valid. If right inverse of B exists, 
and it exists at least classically, it coincides with the left inverse in an associative algebra 
with unit. So we have the identities 

(BB-% = B tk B^ 1 = B lk D~ 1 A jk = ^ (4.38) 
k k 
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We write the solution of eq.(4.37) as 

H = -B~ 1 V (4.39) 

Theorem 12 The quantities Hi defined by eq.(4-39), which solve Baxter's equations 
eqs.(4-36), are all commuting 

[H i ,H j ] = Q 

Proof . Using that V k and V\ commute, [V k , V{\ = 0, we compute 

[Hi, Hj] = ^[S-Vfe,^- 1 ^] (4.40) 

k,l 

= £ [B^'B^VM - B^[B^, V k ]V t + B-\B~ k \ V l \V k 

k,l 

Using 

[A- 1 ^- 1 ] = A- 1 B- 1 [A,B)B- 1 A- 1 = B' 1 A'^A, B]A~ 1 B' 1 

so that 

\ B ik' B Jl'\ = £ B ir lB jJ[ B rs,B r , sl }B-}B-£ 
rs,r's' 

= £ Bj^B^iBrsTBr^^B'^B-,} 

rs,r's' 

the first term can be written 

Y.^ B Tk\Bj l 1 ]V k V l = ^B^Bp[B ra ,B rl ^B^B^ + B^ 1 k B- l 1 )v k V l 

k,l 

= £ 2 B ^ B ^ ^ rS ' Br ' S '} { B sk B s'l + B sl lB s'k) V k V l 

Using that [B rs , B r > s i] = 5 rr >[B rs , B rs /] and is therefore antisymmetric in ss', and setting 

K ss > = £ ( B ~}B; k l + B'lB' 1 - B^B'l - B-^B~]^V k Vi 
k,i 

we get 

E^ifc 1 '- 8 ^ 1 ]^ = ^2^BT r 1 BT r 1 [B rs ,B rs/ }K ss > 

k,l rss' 

= - £ ^ B ]r B ir l \- B rs, B rs i]K ss i 
rss' 

= £ g t" 8 ^ 1 ' B jr] i B rs , B rs '} K ss t 
rss' 
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The last two terms in eq.(4.40) are simpler, we get 



k,l rsk 



The quantities Hi will commute if 



ir ' jr 

This is true as shown in the next Lemma. 



[B- 1 ,B- 1 ] = 0, Vi,j,r 



(4.41) 



The condition eq.(4.41) says that the elements on the same column of -B -1 commute 
among themselves. In a sense this is a condition dual to the one on B. It is true 
semiclassically because 

{Br 1 , Bp-} = B i? B ja'{ B ab, Ba'v}BfrByl = ^ B' 1 Bp{B ab , B ab ,}Bp B^ = 

a,a' ,b,b' a,b,b' 

where in the last step we use the antisymmetry of the Poisson bracket. We show that it 
is also true quantum mechanically 

Lemma 1 Let B be a matrix whose elements commute if they do not belong to the same 
line 

[B ik ,B fl } = Q iii^j 

Then the left inverse B^ 1 of B is defined without ambiguity and moreover elements on 
a same column of B^ 1 commute 



Proof. We want to show that 



denote by ffi the vector with components B k i, k / r. Then we have (with j > i) 



[B-\Bj r l ]=Q 



- A rj B i7 } 



A r jB-} 



-l)^A^A--/3} r, A-/3f A-^ 1 

-^y+i+g-j+1 R (r) A R (r) 



ff> A A • • • ) A • • • /?« A • • • #> A £ fi ] B, 



kr 



. iy +i+g-j+lp(r) A A . . . p (r) A . . . dr) A . . . Mr) A B ~ 



1)^(3^ A 4 r) A • • • ^ A • • • /?« A • • • fflB? 



(r) 
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where in the third line we used eq.(4.38). In the above manipulations, we never have two 
operators Bij on the same line so we can use the usual properties of the wedge product. 
Moreover it is important that the line r is absent in the definition of (3^ r \ Remark that 
this equation can also be written A r iD~ 1 A r j = A r jD~ 1 A r i which is a Yang-Baxter 
type equation. ■ 

With this Lemma, we have completed the proof of our theorem. It is remarkable 
that, again, only the separated nature of the variables Aj, /Ltj is used in this construction, 
but the precise commutation relations between Aj, fii does not even need to be specified. 
This is the origin of the multi Hamiltonian structure of integrable systems, here extended 
to the quantum domain. 
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Figure 4.3: Red dots are the Bethe roots for the one spin system. Green dots are the 
branch points. The thin black curve is the solution of eq.(4.35). {h = 1/30, s = 1/fr, 
M = 4/h, highest energy state). 
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Chapter 5 

The Heisenberg spin chain. 



5.1 The quantum monodromy matrix 



We have seen that in the case of two dimensional integrable field theories, the analog of 
the Lax matrix is the monodromy matrix 



T L (X) = Pexp 



J 



Udx 



(5.1) 



In order to define this quantity at the quantum level, it is convenient to start from the 
discretized version of the theory. Thus, we consider a lattice with TV sites and lattice 
spacing A = L/N . To each site n of the lattice, we attach the local transport matrix 
LnW (we retain traditional notations). The matrix elements of L n (X) are functions 
of the local quantum fields of the model. Over each site of the lattice, there is a local 
Hilbert space H n on which the field operators act non trivially. The total Hilbert space 
of the discretized system is the tensor product of all these local Hilbert spaces. 



N 



n 



n=l 



We now define the quantum monodromy matrix T\r(A) by the discretized version of 
eq.(5.1) 



N 



T N (X) = l[L n (X). 



n=l 



(5.2) 



This is our basic object of study in this Chapter. 
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5.2 The XXX spin chain. 



We now provide the canonical example of this construction: the XXX spin chain. Con- 
sider a lattice with iV sites. On each lattice site n, we attach a Hilbert space Tt n = C 2 . 
The total Hilbert space is of dimension 2^. 



■c 2 N 



On each site we introduce the spin operators s l n acting non trivially on 7i n only where 
they are represented by Pauli matrices 



s n = 1 <g> 1 <g> ■ • • 1 



h 



<7<g> 1 



1 



where the non trivial innsertion is on the ra-th position. We recall that 



i / 1 \ 2 / -j 
° =f l o ' a =\i 



We will also use = ± is 2 . We define 



a = 



P\2 = -(Id Id + a® a) 



On each lattice site we introduce the local L n -operator 



1 
-1 



(1 
















1 








1 








Vo 








1/ 



L n (X) 



A + isl 



1S T 



is 



+ A - is 



A Id + i a s Tl 



(5.3) 



Define now the following matrices of operators, 

Lin{\) = L n (X) 1, L 2n (A) = 1 L n {\) 

where the tensor product now refers to the auxiliary matrix space. It is a simple exercise 
(multiplication of 4 x 4 matrices) to check that it satisfies the relation 



i?i 2 (A - n)L ln (X)L 2n {fJ.) = L 2n (lJ.)Li n (\)Ri2{X - n) 



(5.4) 



with 



RMX) 



1 

X + ih 



(A id + ihP 12 ) 
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or explicitly 



with 



R12W 



6(A) 



(I 0\ 

c(A) 6(A) 

6(A) c(A) 

\0 1/ 



c(A) 



A 



A + ih A + ih 

Note that the so called ultralocality condition holds 



L ln (X)L 2m (fi) = L 2m (fj,)L ln (X) n^m 



We now construct the monodromy matrix T/v(A) for a lattice of N sites 

N 



T N (X) = J] MA) 



n=l 



As in the classical case, one can go from the local formula eq.(5.4) to a global one for 
the quantum monodromy matrix. 

Proposition 31 If L n (X) is such that (5.4) and the ultralocality condition holds, then 



R l2 (X - /i)Ti(A)T 2 (//) = T 2 ( M )Ti(A) J Ri 2 (A - fi) 



(5.5) 



Proof . We use ultralocality to write 

Ti(X)T 2 (fi) = H[L ln (X)LM] 

Then, using (5.4), we find 



N 



n=l 



N 



i?i 2 (A-/x)Ti(A)r 2 (/x) = \{[L 2n {^)L ln {X)}R 12 {X- M ) 

n=l 

= T 2 (/i)T 1 (A)i? 12 (A-/i) 
In the last step we have used again the ultralocality property. 



Remark. This is a quantum analog of the classical formula 
{71(A), T 2 ( M )} = [r 12 (A, /[ i),T 1 (A)T 2 (^)] 
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to which it reduces in the semi classical limit 



(5.6) 



Remark. Since i?i 2 (A — fi) depends only on one variable, A — fi, the same argument 
works for inhomogenous models, defined by 

T(A; ai, ■ ■ • , a N ) = Li(A - ai)L 2 (A - a 2 ) ■ ■ ■ L N (X - a N ) 

From these commutation relations we can extract a family of commuting Hamiltoni- 
ans. Define 



t(X) = Tr T(A) 



then 



[t(\),t(»)] = 



This is because 

Ti(A)T 2 (/i) = R\i (X - /x)r 2 (/x)Ti(A)i?i2(A - /*) 

then take the trace and the cyclicity property to bring i? _1 (A — /ix) to the right where it 
cancels the R(X — fi) factor. Then use the fact that 

Tn 2 Ti(A)T 2 (^) = t(X)t(fi), TYi 2 T 2 (/x)Ti(A) = t(fi)t(X) 



Local Hamiltonians are obtained by expanding around a point A = Ao such that 
£i2(Ao) = P12, the permutation operator. To see that we get local quantities by expand- 
ing around Ao, we write Li 2 (A) = Ylij e ij ® L %3 \\). Then 



L 



l N~l l N 

N 



The point Aq is such that L ij (Xq) = eji. Replacing L l; >(X)\\ = eji, we get 



t(X)\ Xo 

t-\x)\ Xo 

t'(X)\ Xo 
so that 



E e i 

H,---,i N 

£ e{ Njl 
,—Jn 

E e i 



JU2 



= *JV*JV-1 

'TV 



3n-i3n 



■ ■ ■ e, 



*«-l*n-2 



n-1 



)' ® <+T 



• e 



jN^N-l 

N 



Il, — ,»JV 



^logt(A)| >0 =t- 1 (A)f(A)| Ao =EE e n J ( L n)'®<+! 
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In this expression, we see that only nearest neighbours sites are coupled. More generally, 
a derivative of order p couples p + 1 sites. In our case 



(5.7) 



L ln (X) = (A - f ) Id + ihP ln = (A + f ) R ln (A - f ) 
so that 

Lin (y) = ifiPln 

Expanding around that point we get a local Hamiltonian. 

n 

This is the Heisenberg spin chain Hamiltonian. 

Another remark obvious from eq.(5.7) is that L(\) and i?(A) are essentially the same 
thing and eq.(5.4) is identical to the Yang-Baxter equation for R(X). 

The product of operators is associative. Writing this condition on eq.(5.5) puts a 
constraint on the matrix R(X — fi). 

Proposition 32 A sufficient condition for associativity is that the matrix R(X — fi) 
satisfies the Yang-Baxter equation 



-Rl2(Al — A2)i?13(Al — As)i?23(A2 ~ A3) — i?23(A2 ~ As)i?i3(Ai — As)i?i2(Ai — A2) 



(5.8) 

This is an equation in End(V\ ® V2 <8> V3). The notation R12 means as usual that it acts 
non trivially on the space Vi <8> V2, and is the identity on V3 ... 

Proof . One should check the equation by a direct calculation. However there is a good 
reason for this equation to hold. Consider the product Ti(Ai)T2(A2)l3(A3). There are 
two ways to bring it to the form T3(As)T2(A2)Ti(Ai) with the help of eq.(5.5 ): 

T T T ^ 12^1 3 rj-i rj-i rj-i 

Rl2 -Rl2-Rl3-R23 

/ \ 
TiT 2 T 3 TgTaTi 

R23 R23R13R12 
\ / 

rj~i rj-\ rj-i ^ 2 3 ^1 3 rj-i rj-\ rjn 

Following the upper path of the diagram produces the combination of i?-matrices in 
the left hand side of eq.(5.8), and the lower path produces the right hand side of this 
equation. The Yang-Baxter equation ensures the commutativity of the above diagram, 
which itself reflects the associativity of the product of the T's. ■ 
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Remark. In the classical limit eq.(5.6), the Yang-Baxter equation becomes the 
classical Yang-Baxter equation 

[n2,n 3 ] + [ri2,J"23] + [»"13,?"23] = 



We end this section by listing a few supplementary important properties of the mon- 
odromy matrix T\r(A) of the XXX spin chain. Writing 

m _(A(\) 5(A) \ 

we see that A(X) and D(X) are polynomials of degree N in A, while B(X) and C(A) are 
polynomials of degree N — 1. In fact we have 



T N (X) = A^Id + iX*- 1 a-S + 0(X 
where S is the total spin operator 

n 

Also, we have the important conjugation property (for real A) 

B\X) = -C{X) 



N-2\ 



We note the obvious relation (SU(2) symmetry) 



which implies as well 



[L n (X),^* a + s a n }=0 



[T N (X),^a a + S a }=0 



Finally, we have the relation (quantum determinant) 

a 2 L n (A)a 2 L^(A - ih) = (A(A - ih) + h 2 s(s + l))Id 
which implies the same type of relation on T(A) 



a 2 T(X)a 2 T t (X - ih) = (A(A - ih) + ti 2 s(s + l)) JV Id 



(5.9) 



(5.10) 
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5.3 Algebraic Bethe Ansatz. 

One can use the fundamental commutation relations 

i?i2(A,Ai)Ti(A)r 2 (/x) = T 2 (ii)T 1 (X)R 12 (X,n) 



(5.11) 



to diagonalize t(X) = TrT(A). The method applies to the case of 2 x 2 matrices. Larger 
matrices will be treated in the next Chapter. Let 



T ^-{C(X) D(X)) 



and let us write the i?-matrix in the form 



(\ 

c(A,/x) 6(A,/i) 

6(A,/i) c(A,/i) 

\0 






1/ 



so that the formulae will apply to cases more genaral than the XXX spin chain. Writing 
explicitly eq.(5.11) we obtain the following set of commutation relations: 



[A(A),A(/.)] 
[D(X),D(ri] 
[B(\),B(»)] 
[C(X),C(ri] 
B(X)A^) 
B(fj)D(X) 
C(A)A( M ) 
C(riD(\) 
c(A,/i)[C(A),B(M)] 
c(X,^[A(X),D(^] 








6(A,/x)S(/x)A(A) + c(A, M )A(/x)B(A) 
6(A,Ai)5(A)D(M) + c(A,Ai)I>(A)B( M ) 
c^AJA^CCAJ + ft^AJC^ACA) 
c( M ,A)Z)(A)C(m) + 6(^A)C(A)D(m) 
-6(A,Ai)(A(A)D(/x)-A(/x)£>(A)) 
-6(A, M )(CXA)£( M )-CV)B(A)) 



Although all these relations are useful, we will need here only three of them which we 
rewrite in the convenient form: 



A(X)B^) 



D{X)B{n) 



[B(X),B(li)] = 

" B(»)A(X) - b ^B(X)A(») 



c(fx, A) 
1 



B(»)D(X) 



c(/i, A) 
6(A >A i) 



B(X)D(fi) 



(5.12) 
(5.13) 

(5.14) 



c(A,/i) c(X,fi)' 

The idea of the algebraic Bethe Ansatz is to use the operator B(fj,) as a creation operator 
to generate eigenstates of the operator 



t(X) = A(A) + D(X) 
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We recall that this operator is a generating function for the commuting Hamiltonians. 

We need a reference state to start with. On each site, introduce a local vacuum |0} n 
in the local Hilbert space H n - 

|0)n 1 1 





The action of the local operator L n {\) eq.(5.3 ) on |0) n is given by 



L n {X)\Q) n 



a(A)|0) n (L n (A))i 2 |0) r 
<5(A)|0) n 



The action of the monodromy matrix T(A) on the state 

|0> = |0)i® |0) 2 ® \0)n 

is obtained by taking the product of these local triangular matrices. Therefore, it is also 
triangular and the diagonal elements are just the products of the local diagonal elements 



T(A)|0) 



a(X) N \0) B(X)\0) 
S(X) N \0) 



Thus we have 



A(X)\0) = a(X) N \0) 
D(X)\0) = 6(X) N \0) 
C(A)|0) = 

Notice that |0) is an eigenstate of 9(A) = A(X) + D(X). We shall now generate new 
eigenstates by applying operators B(fii) to |0). 

Proposition 33 The vectors fl(/J,i, ... , \xm) defined by: 

• • ■ ,Mm) = B{^)B^ 2 ) . . . B( m )\0) (5.15) 

are eigenstates of t(X) if the following relations are satisfied : 

S(H)J ~\\c{^M) (5 ' 16) 

The corresponding eigenvalue of 9(A) is given by 

M M 
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Proof . To evaluate the action of the operator A(X) on the vector . . . , /jm), we use 

the commutation relation eq.(5.13). When pushing A(X) to the right of the chain of 
B(fij) operators, we generate 2 M terms, coming naturally in M + 1 groups 



M 

A{X)B{fi 1 )B{^ 2 )...B{ m )\0) = A(A;/ii,...,/i M )n B (w)l ) 

1=1 

M M 

+ ^A J (A;/xi,...,/x M )5(A)n^(w)|0) 

where A and Aj are numerical coefficients. The first term is of the form required for an 
eigenvector, and we call it a wanted term, the other ones are the unwanted terms. 
The wanted term is obtained by using only the first term in eq.(5.13) for the commutation 
of A(X) through the B((ij). Otherwise one of the B operators would carry the argument 
A. Since this is the only way to obtain it, this gives immediately 

M 



A(A; in,..., mm) = a N {X) TT 1 

fJl c(w,AJ 

To obtain the coefficient Aj(X; . . . , hm) of the unwanted term, we first bring the 
operator B(/ij) in the first position in the chain of S's (the B operators commute by 
eq.(5.12)). In the first commutation of A(X), we use the second term of eq.(5.13). The 
operator A carries now the argument fij and must keep it until hitting the reference state 
|0). Hence we commute it using only the first term in eq.(5.13). This gives uniquely 

A, (A; M ) = -^ f N (H) ft ^) 

Similarly, we have 

M 

D(X)B( f i l )...B( m )\0) = A'(X; f i 1 ,..., m )l[B(n l )\0) 



M M 

^A;.(A;/i 1 ,..., / u M ) J B(A)n 5 (w)|0) 

3=1 1=1 



with 



M 



A'(\; fr,..., m ) = 6 N (X)l\ 



* c(A,/ij) 

A^X;^,..., m ) = - b -^\s N ( H )Y[-^— 
J c(X,Hj) A A c(/Xj,w) 
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The vector £l(fxi, . . . , hm) will be an eigenvector of 9(A) if all the unwanted terms van- 
ish. This gives the conditions Aj + A'j = 0. Since the function 6(A) /c( A) is an odd 
function of A, these conditions reduce to eq.(5.16). The corresponding eigenvalue is sim- 
ply A(A; Mi, ... , n M ) + A'(A; m,..., fi M )- ■ 



Let us write these equations explicitly in the case of the XXX spin chain. The Bethe 
equations take the form 



- v, 



N 



n 



f*j - (JLk-ifi 



(5.18) 



and the corresponding eigenvalue is 



t(X;{ H }) = A + 



ih 



N 



ih 
2~ 



N 



n 



A — fJ-k + 
A - 



(5.19) 

A special property of the XXX spin chain is that the Bethe eigenvectors are all highest 
weights vectors for the total spin operator S. To prove it we start from eq.(5.9) which 
implies 

[S 3 ,B(X)] = -hB(X), [S+,B(X)] = ^(A(X) - D(X)) 
The first relation implies immediately 



S 3 tt(m, ■ ■ -,hm) = ( y - M ) Qf// 1 . ■ ■ ■ .//_ w ) 



To prove that 



S + tt(/j,i, ■ ■ -,hm) = 



we use the second relation to get 

s + n(n 1: ■■■,hm) = \Y, b ^) ■ ■ ■ b(h-i)(Ah) - d(h)) b (h+i) ■ ■ ■ B(»m)\0) 

j 

when we push the oprators A and D to the right, we get a sum of terms of the form 
S + Q(m, ■■■ , hm) = \ ^2 r i(^i' ■ ■ ■ ' ^m)B(^i) ■ ■ ■ B(fij) ■ ■ ■ B(fj, M )\0) 
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Let us compute Ti(fii, • • • , /j,m)- The only way to get this term is to start from 

(A( f i 1 )-D( f i 1 ))B(n 2 )---B( m )\0) 

and push A(jjL\) and D{ni) to the right using only the first terms in eq.(5.13,5.14). Hence 
we find 

r l(w , • • • , , M ) = a"^) n - g 3^ 

This vanishes if the Bethe equations are satisfied. Clearly the same is true for the other 
T j (fx 1 , ■ ■ ■ , fi M ) by symmetry 



5.3.1 Baxter equation. 

Let us introduce the polynomial 



M 



Q(X) = H (A - nr, 



m=l 

Then the Bethe equations eq.(5.18) ca be rewritten as 

ih\ N ( ih\ N 

I 1 1, + 2") <3(Atfe -^)+(^--J Q(nk + ih) = 

This means that the polynomial of degree N + M 

A + y) Q(\-ih)+(\- l -£\ Q(X + ih) 
is divisible by Q(X). Hence there exists a polynomial t(X) of degree N such that 



A + y) Q(A - ih) + ( X - y') Q(X + ///)- /(A)Q(A) 



A? 



This is Baxter's equation. The polynomial i(A) is the same as in eq.(5.19) because that 
equation can be rewritten as 



t(X;{ H }) 



ih\ N Q(X-ih) f ih\ N Q(X + ih) 
+ 2j Q(A) + V ~2J Q(X) 



hence the coefficients of this polynomial are just the eigenvalues of the set of commuting 
Hamiltonians. 

Just as in the Gaudin model, it is interesting to introduce the Riccati version of this 
equation. We set (do not confuse this S with the total spin !!) 



S(X) 



Q(X - ih) 

Q(X) 



127 



Then Baxter equation becomes 



A + 



ih 



N 



ih 



N 



5(A)+(A-y) S~ 1 {X + ih) = t(X) 



This equation determines both 5(A) and t(X). To find the equation for t(X), we expand 
around A = —ih/2 getting 

(e - ih) N S-\e + ih/2) = t(e - ih/2) - e N S{e - ih/2) 

Similarly, expanding around A = ih/2 we get 

(e + ih) N S{e + ih/2) = t{e + ih/2) - e N S-\e + 3ih/2) 

Multiplying the two, we find 




(5.20) 

This is a system of N equations for the N + 1 coefficients of i(A) which determines it 
completely if we remember that t(X) = 2X N +0(X N ~ r ). In fact the sub leading coefficient 
is also easy to compute 

/ N(N — 1) 



t(X) = 2X 



N 



h l 



M(M - N - 1) ) A^" 2 + 0{X N ~ :i ) 



,AT-3\ 



For chains with a small number of sites N, we find by solving directly eq.(5.20) (we set 
^=1): 



N 


t(X) 


M 


S :i 


mult 


dim 


2 


2X*-\ 





1 


3 




2 


2A 2 + | 


1 





1 


2 2 = 3 + 1 


3 


2A 3 - |A 





3 
2 


4 




3 


2A 3 + |A - & 


1 


1 

2 


2 




3 


2A 3 + |A + ^ 


1 


1 

2 


2 


2 3 = 4 + 2 x 2 


4 


2A 4 - 3A 2 + | 





2 


5 




4 


2A 4 + A 2 -| 


1 


1 


3 




4 


2A 4 + A 2 - 2A + | 


1 


1 


3 




4 


2A 4 + A 2 + 2A + I 


1 


1 


3 




4 


2A 4 + 3A 2 - | 


2 





1 




4 


2A 4 + 3A 2 + f 


2 





1 


2 4 = 5 + 3x3 + 2x1 
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This construction can be generalized to the case of a spin-s chain. The Baxter 
equation then reads 



(A + ihs) N Q(X - ih) + (A - ihs) N Q(X + ih) = t(X)Q(X) 



The Riccati equations reads 



(A + ihs) N S(X) + (A - ihsf S-\X + ih) = t(X) 



For s = 1 for instance we expand around A = ih, X = 0, A = — ih to get 

(e + 2ih) N S(e + ih) = t(e + ih) + 0{e N ) 

{e + ih) N S{e) + {e-ih) N S- 1 {e + ih) = t(e) 

(e-2ih) N S- 1 (e) = t(e - ih) + 0(e N ) 

from which we deduce (s = 1) 



(5.21) 



t(e + ih)t(e)t(e - ih) = (e - ih) N {e + 2ih) N t(e - ih) + (e + ih)* (e - 2^) Jv t(e + ih) + 0(e JV ) 



Clearly, for a spin-s, s > 0, the degree of the equation is 2s + 1. If however s < the 
equations generically do not lead to a finite degree equation. 

5.3.2 Separated variables. 

Suppose for a while that the spin in eq.(5.3) is classical. Then the spectral curve reads 

det(T(A) -n) = n 2 - t(X)n + (A 2 + s 2 ) N = (5.22) 

where t(X) = A(X) + D(X) = 2X N + 0(X N ~ 2 ). This is a hyperelliptic curve. To put it in 
canonical form, we set y = ji — \t{X) so that 



1 



= ^(X) - (A + sT ~ (t 2 - Ns')X^- 2 + ■■■ 

from what we see that the genus is N — 2. Remark that the dynamical moduli Hi appear 
linearly in t(X). If we vary them, we get 



5fi 



dX 



8t(X)dX 



St(X) 



dX = holomorphic 



H ^-(A^s 2 )^" 1 2y 

Hence the natural (and correct !) Poisson bracket for the separated variables is 

{Xk,^k>} = $kk'Vk (5.23) 



129 



Setting otherwise 

the spectral curve becomes 



H = (X + is) N S 



(A + is) N S + (A - is) N S- L = t(X) 



We see that the Riccati equation eq.(5.21) can be viewed as a deformation of the spectral 
curve. The fact that this deformation involves difference operators instead of differential 
operators can be understood because the natural quantization of eq.(5.23) is that fj, k is 
a shift operator 



Mfc^fc = (Afe - ih)nk 



(5.24) 



We now prove this formula in the quantum XXX spin chain. Following Sklyanin, we 
introduce the quantum separated variables as the zeroes of B(X). 



N-l 

B(X) = iS- H(X-Xk) 

fc=i 



Because [B(\),B(X')] = 0, we have 



[5-,A fc ]=0, [A fc ,Afc/] = 



For an operator X(X) = Yln-^-nX 71 , we define 

M[X(X k )} = J2^Xn 

n 

where the operator A^ is ordered on the left. Substituting jjl —>■ Xk according to this rule 
in eq.(5.13) which we write in the form 

(A - n)A(n)B(X) = (X-n + ih)B(X)A(n) - ihB{^)A(X) 

we get 

AT[(X - X k )A(Xk)]B(X) = AT[(X -X k + ih)B(X)A(X k )} - iW[B(X k )]A(X) 

Now J\f[B(X k )] = B(X k ) = because all coefficients of B(X) commute with Afc. Next, 
obviously M[(X — X k )A(X k )] = (A — X k )ftf[A(Xk)]- Finally again because B(X) commutes 
with X k , we have Af[(X- X k + ih)B(X)A(X k )] = (A- X k + ih)B(X)M[A(X k )} Hence, setting 

= M[A(X k )} 
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we get 



and therefore 



(A - A fc ) B(\) = (A - X k + ih)B(X) 



(A - A fc ) S- H(\ - A,) = (A - A fc + ih)S~ ]J(X - A>£ 



(+) 



By eq.(5.9), we have 



which implies also 



[S-,A(\)] = -B(X) 



(+) 
k 



Simplifying on the left by (A — X k )S~, we get 

4 +) I[( A " A i) = ( A " A * + ^) JJ(A - A>| 

Mi +) Afc = (Afc-^)4 +) 
at least when fi k + ^ acts on symmetric functions of the Afc's. By exactly the same argu- 

JSr[D(A fc )] 

/ii _) Afc = (Afc + ^)/x^ _) 



this means 

at least wh< 
ment, defining 

we show that 



From their definitions, we have 

^+)_ (A A^(Afc) + J D(Afc))] + Mfe - ) =0 
To finish the identification of /ij^ with fi ±:L in eq.(5.22), we show that 



(5.25) 



/4 +) /4 } = (Afc — ihXk + K z s(s + 1)) 



A? 



(5.26) 



We start with the identity eq.(5.10) which implies in particular 

A(X)D(X - ih) - B{X)C{X - ih) = (A(A - ih) + h 2 s{s + 1)) N 
Now, we have 

Af[A(X k )D(X k -ih)} = ^(Afc - ih) n J\f[A(X k )]D n 

n 

= J>fc - ihr^D n = X>< +) AJD B = Mfc +) A/p(Afc)] = M W^~) 
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More simply, we have J\f[B(\k)C(\k — ih)] = 0. This completes the proof of eq.(5.26). 
We can solve eq.(5.26) by redefining 

= (A* + itu>) N S k , 4 _) = ( X k ~ ihs) N S^ 
where Sk is the shift operator for the variable A& 

Sk^k = (Afe - ifo)Sk 

Then eq.(5.25) becomes 



(A fc + iters* + (A fc - fc)"^ 1 = AA(t(A fc )) 



As in the Gaudin model, the Bethe states can be expressed in terms of the separated 
variables. The formula is 



n(n 1 ,---, m ) = (iS-) M l[Q(\ k )\o) 



At this stage, what seems to be still missing is the scalar product in this representation. 
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Chapter 6 

Nested Bethe Ansatz. 



There exists a non trivial generalization of Bethe equations to the case of SU (n) spins. 
We present here the algebraic version of this construction. 



6.1 The i?-matrix of the Affine sl n+ \ algebra. 

Our starting point will be the i?-matrix associated to the quantum loop algebra U q {sl n+ i® 
C(A, A -1 ). It can be written as 



R{\) 



XR — X~ X R 
Xq-X-iq- 1 



with 



R = ^2 en ® ejj + q ^ en ® en + (q - q x ) ^ ejj <8> ejj 
R = ^ en ® ejj + q" 1 en <g> en — (q- g^ 1 ) ^ ejj <g> e 



Here denotes the matrix with elements [eij] a f3 = bia^jp- The matrices R and R are 
the i?-matrices of U q (sl n +i). 

The matrix R(X) satisfies the Yang-Baxter equation. 

Let us write explicitly the matrix elements of R(X). They read 



R ab( X ) - 



A — A" 



Xq-X^q 



1„-1 « 



5? 8? 



1 



A- A" 



Ag- A-V 



l/i-i 



q-q 



Xq — X l q 



—X< a - h \l-5 ab ) 
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6.2 sl n+ i generalization of the XXZ model. 

We use R to define the matrix of statistical weights of a vertex model, and we construct 
as usual the transfer matrix 

T 1N (X) = R U (X)R 12 (X)---R 1N (X) 

This matrix satisfies 

R 12 (X/n)T 1N (X)T 2N (fi) = T 2N (^T 1N (X)R 12 (X/n) (6.1) 
As usual, the traces TJv(A) = triTuv(A) commute 

[T N (X),T N (»)} = 

Since R(X)\\ = i = P, the permutation operator, we can construct local commuting quan- 
tities 



Proposition 34 Let 



R= q ~ q' 1 x dlogT N (X) 



dX 
we have 



A=l z 



N 



n=l 



E + {l± f 1 E « +1 - E - 



Proof . We write i?i 2 (A) = ^ <g> L^'(A) with 

then remark that L l i(X)\\=i = eji, so that the general construction applies: 

T N \X)U(X)\ X=1 = EE e «'( L «)'® e n+l 

n ijk 



This Hamiltonian acts on the Hilbert space 

N 

nN = Y[®hj; h j = C n+1 
For n = 1, we recover the XXZ model. 
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6.3 Commutation relations. 



We can arrange the matrix elements of R(X) as follows 



R(X) 



where 



^ A ^' = Xq- X^q- 1 ^ ^ X ^ ij = Xq 1 -X- 1 q- l X ^ ^ X ^ ij = \ q - PyT j 'j 

Since b, b, c are all proportionnal to the identity matrix, we will treat them as ordinary 
C-numbers in what follows. 

To take advantage of the block structure of the i?-matrix, we also decompose the 
matrix T(A) in a similar way 



11 lj Jl n 
11 /1 \ 
\i c(A) 6(A) 
il 6(A) c(A) 
w V RW J 



q-q 



A — A" 



(A(X) B(X) A 



where 



A(A) = T n (A); ^(A) = T H (A); Q(A) = T a (A); A; (A) = T^(A) 

Thus -B(A) is a line vector and C(A) is a column vector of dimension n. D(X) is a n 2 x n 2 
matrix. The fundamental equation eq.(6.1) yields the following relations 



D 1 (X)B 2 (ti) 



' B(»)A(X)- b ^-B(X)A( fl ) 



c(X,n) 

B 1 (X)B 2 (^) = B 2 (^B 1 (X)R ( S(X,^) 
We will use the second relation eq.(6.3) in the form 



c(/x,A) c(n,xy 

-^—B 2 ( f i)D 1 (X)R^(X, fJ ,) - h S^B l {X)D 2 { l x) 



c(X,n)' 



D 1 (X)B 2 ( fi ) 



-B^D^R^iX^) 



b(X,fi) 



c(X,fi) ^' ^ ' 12 ^ c(A, M ) 
where P12 is the permutation operator Pjjj? = 5^6 ; j. 
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B 2 {X)Di(n)Pi 2 



(6.2) 

(6.3) 
(6.4) 



6.4 Reference state. 



To apply the algebraic Bethe Ansatz, the first step is to find a reference state to start 
with. In our case this is provided by- 



Proposition 35 Let |1) be the vector 



ID = ID 



IN, 



where ID 





w 



The action of operatorial entries ofT(X) on |1) is given by 

A(X)\1) = a N (X)\l) 
D(X)\1) = 5 N (X)Id\l) 
C(X)\1) = 



where 



Proof. We have 



a(A) = 1; 5(X) = c(A) = 



A - A 



-i 



Xq-X^q 



l/T-l 



4 1 (A)|l>n 
^(A)|Dn 





a(A)|l)„ 

*(A)|l>n 



mX)\l) n = 0, i^j 
The result is obtained by multiplying these triangular matrices at each site of the lattice. 



The idea of the algebraic Bethe Ansatz is to look for eigenstates of the form 

ii—ip 

Since B is seen as a line vector, we consider X as a column vector. We can rewrite our 
state in a tensor notation 

y x ({ri) = B 1 (vi)---B p (v p )\l)X 

We want to find the action of 

T(A) = A{X) + Tr D(X) 
on this vector using eqs. (6.2,6. 3,6.4). 
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Proposition 36 The action ofT(X) = A(X) + Tr D(X) on the vector is of the 

form 



T(A)* X (W) = B 1 (vi)---B p ( f i p )\l)Y 



(6.5) 



+ ^Bi(\)B i+1 (fH+i) ■ ■ ■ B p (fi p ) ■ ■ ■ B^i^VjYi (6.6) 



The first term is called the wanted term, Yq is given by 



(6.7) 



The other terms are called the unwanted terms. Setting them to zero yields the equations 




(6.8) 



The matrix T(A, {/x}) appearing in these equations is defined as 



T(A, M) = TYotfSWi) • • • < } (A//i P ) 



JTiis is t/ie transfer matrix of a spin model where the spin takes n values instead o/n + 1. 

Proof . Wanted terms. Let us start with A{\). The wanted term is obtained by commut- 
ing A{\) through the .B's using only the first term in the commutation relation eq.(6.2). 
We get 

^(A)*x({/X})|wanted = <* N (\) f[ *x(W) 

Let us now evaluate Tr D(A)*x({m}) Iwanted- This term is obtained by commuting D(\) 
with the B's using only the first term in eq.(6.3). We have 

A)(A)*x(M) = fl -j^—M^) • • • BMD (\)\1)B$(\, mi) ■ ■ ■ <?(A, Mp )X 

P 1 

= ^ II "7x \ Bl (Vi) • • • ^(^)|1)< } (A, A*i) ■ ■ ■ < } (A, M P )X 
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and taking the trace, we get 

p i 

Tr £(A)*x(M)|wantcd = 5 N (X) [] -jr-^rBM ■ ■ ■ B P (^)\1)T(X; {fi})X 

Adding these two contributions, we obtain exactly the wanted term and the expression 
forlo- 

Unwanted terms. The unwanted term Y\ has a contribution proportional to a N (fii) 
coming from the commutation of A(X), and a contribution proportionnal to 6 N (fii) 
coming from the commutation of Tr D(X). 

To calculate the first contribution, we commute first A(X) with B([xi) using the 
second term in eq.(6.2), and then we commute A(fii) with the other i?'s using only the 
first term in eq.(6.2). We get 

A(A)*x(M)|u„wa. = f[ -T^—MX)B 2 {H2) ■ ■ ■ BM\\)X 

Let us now evaluate the term proportional to 5 N (fii) It is obtained by commuting 
first D (A) and B{^l\) using the second term in eq.(6.3) and then pushing D(fii) through 
the right using only the first term in eq.(6.3). We get 

D (A)*x(M) = -^^[[^-tBiIA^^)-^)!!) 

^02 (^i>M2) • • • (Ml. A*p)^bi-^ 
taking the trace Tro and identifying Poi with /2q^(//i, /xi), we get 

trD(A)*x(M)Uwa. = -^(Mi)^4 f[ -T^-^^ilA)^^) • • • B p (^,)|l)T( Ml ; W)X 
Putting all this together we obtain 

* = -^m 6 ^ n * - ^m*tv4 n -^^ r ^ 

c(Mi,A) f^c(A*i,Mi) c(A,/ii) f^c(/ii,A*i) 

It remains to examine what happens for an unwanted term where A takes the place 
of fj,i,i > 2. It is enough to assume that it is fx 2 - To compute it, we first push in 
the last position using eq.(6.4). Thus, we write 

= 1] Mto) ■ ■ ■ B^BMiyiiMfauto) ■ ■ ■ RMfaupjX 

Using the relation TroPoiM)2 = ^fi2 we can rewrite this formula as 
*xm) = ]jB 2 ^ 2 )---B p ^ p )B 1 ( f i 1 )\l)T^ 1 ;{ f i})X 
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Performing exactly the same analysis as before, the unwanted term takes the form 

B 2 (X)B 3 (fi 3 )---B p ( f i p )B 1 ( f i 1 )\l)Y 2 

where 



Y 2 



and since T(n2\ {h}) commutes with T(fj,i; {/J,}) we can rewrite Y2 as 
y 2 = -T(mi;M) 



AT, ^%^A) tt 1 



C(M2,A) c(^i,/i 2 ) 



II ~T~~ — ^(^2; W)x 



Notice that when we set the unwanted terms to zero, the A dependence drops out because 

c(fH,X) c(X,ih) 



6.5 Bet he equations. 

The vector X must be a simultaneous eigenvector of the matrices T(X, {//}) for arbitrary 
values of A. These equations are compatible since 

[T(A,M),T(A',M)]=0 

We are back to the original problem of diagonalizing a tranfer matrix for a model 
where the spin takes n values instead of n + 1. Therefore, the solution will be obtained 
by repeating n-times the procedure. 

At each step we introduce a set of /i's 

= {/4 m) ,..., M M} m = l,---,n 
To the set {// m )} is associated the transfer matrix of an inhomogeneous vertex model 
TM(A, {„<■»>}) = Tro#(A/ M S^)4?(A/^) • • • R^JX/^) 
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In the matrix the spin takes the values m + l,---,n+l. In particular, (A) = 1 
is a C-number. 

The eigenvalues A( m )(A, {/i (m) }) of T (m \\ {^ rn) }) depend on all the subsequent 
k = m + 1, ■ • • n. 

Proposition 37 The sets of numbers {^ m ^},m = 1, • • • ,n are determined by Bethe's 
equations 

,(1), 11 J 11 , (1) „(!), 



Pfc + l Pfc-l ! Pfc „(',,( fe ) ,,( fc )^ 

TT c (M (fc+1) ^ (fc) ) TT = TT j — ' 

11 ; ii (fc) (fe _ 1} 11 , (*) (*). 

j=l j=l C W > J j^i C (A > fJ'j ) 



"fr 1 1 fic(Mf ,/4 w) ) 

^ rf// (n) // (n_1) l ^rf// (n) // (n) l 

j=l , V-j ) j^i , ) 

The eigenvalues are given by solving the set of recursion relations 

Pfe+i 

A«(A, {^)}) = II (fc+1) (6.9) 
i=i c(/ii %A) 

Pfc Pfc+i 1 

+n^."f)n- r 4wK A(l+1) ( A '{" ( ' +1, » 

j=i i=i cy^iHi ) 



Pn 1 Pn-1 Pn 1 

A«»-)(A, = n T ^rr;+n^fr i, )n7AoT 

i= i c(//> %A) j=1 i=1 c(A,/i) y ) 

t/iese equations completely determine the 's once the {^}'s are known. 

Proof . Eqs.(6.9) are just coming from the expression of Yq eq.(6.7) in the k th step and 
remembering that a = 1 and S N (X) should be replaced by Hj=i C (A smce we are 
dealing with an inhomogeneous model. The expression for A^™ -1 ** takes into account the 
fact that in the last step, R^ n \\) = 1 is a C-number. 

The {/i( fc )}'s are solutions of Bethe's equations eq.(6.8), which in the k th step read 

Pfe + l -, Pk Pk+l 1 

n - n*r o "f > n . w ^. ^n^ 

rfi c(n) ',fi\ ') j=1 m c{n\ ') 
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Now, since c(A, A) = 0, from eq.(6.9), we have 



Pk+i 



A(*)( A 4 fc ) > { M w»=n 



so that Bethe's equations become 

Pk + l j Pfe Pfc + l j Pfe + 2 

n ; (fc+i) (fc+iK = n c ^ ) n n , ^+2) 
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